Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_013198 2,160,010 G→A A459T (GCG→ACG)  LGG_02095 type II restriction enzyme methylase subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0131982,160,0100G→A100.0% 926.5 348A459T (GCG→ACG) LGG_02095type II restriction enzyme methylase subunit
Reads supporting (aligned to +/- strand):  new base (148/136):  ref base (28/34):  total (176/172)

Only 100 of 382 total aligned reads displayed.

CGGCTGGTCGATTTTGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAGTTGA  >  NC_013198/2159914‑2160110
                                                                                                |                                                                                                      
  gCTGGTCGATTTTGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACt                                                                                                <  18:299117/101‑1
  gCTGGTCGATTTTGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACt                                                                                                <  3:299117/101‑1
    tGGTCGATTTTGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACttt                                                                                              <  18:1054385/101‑1
    tGGTCGATTTTGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACttt                                                                                              <  12:924354/101‑1
    tGGTCGATTTTGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACttt                                                                                              <  3:1054385/101‑1
      gTCGATTTTGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTga                                                                                            >  19:2738/1‑101
      gTCGATTTTGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTga                                                                                            >  4:2738/1‑101
        cGATTTTGATTCACAGGCACATAAGGCAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTgaga                                                                                          <  20:649909/101‑1
          aTTGTGATTCACAGTCACAGAAGGAAAAAGCTTATCTTCCTTACGTTAAAGATCCAAATCTTAAGTATGGTTCTCGTACTAAACAAAC..GAACTTTGagaag                                                                                        <  18:567240/101‑1
          aTTGTGATTCACAGTCACAGAAGGAAAAAGCTTATCTTCCTTACGTTAAAGATCCAAATCTTAAGTATGGTTCTCGTACTAAACAAAC..GAACTTTGagaag                                                                                        <  3:567240/101‑1
            tttGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGAt                                                                                      <  26:1164660/101‑1
            tttGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGAt                                                                                      <  9:164206/101‑1
              tGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATAc                                                                                    >  11:361639/1‑101
              tGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATAc                                                                                    <  11:202927/101‑1
              tGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATAc                                                                                    <  27:789601/101‑1
              tGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATAc                                                                                    <  17:534964/101‑1
               gATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATAcc                                                                                   >  27:274585/1‑101
               gATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATAcc                                                                                   >  26:533796/1‑101
                aTTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAACGATCCAAGTCTTAAGTATGTTTATCGTACTCAACAAAC..AACATTTGAGTAGAGACCt                                                                                  >  6:487316/1‑101
                 ttCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTg                                                                                 >  28:341518/1‑101
                   cacaGGCACATAAGTAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGa                                                                               <  26:204378/101‑1
                   cacaGGCACATAAGAAAAAATCTGATCGTGCGTCCCTTAAAGATCCAAAGCTTAAGGATGTTGAGCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGa                                                                               <  2:104023/101‑1
                      aGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGt                                                                            <  10:547597/101‑1
                      aGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCAAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGt                                                                            >  14:190313/1‑101
                       ggCACATAAGGAAAAGGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTc                                                                           <  27:1204676/101‑1
                       ggCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGGTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTc                                                                           >  2:216858/1‑101
                          acaTAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGa                                                                        >  21:378981/1‑101
                           caTAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGAt                                                                       >  16:376995/1‑101
                           caTAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAAGCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGAt                                                                       >  20:371602/1‑101
                            aTAAGGAAAAAGCTTATCTTGCTGGCGTTAAAGATCCAAATCTTAAGTATGTTTATCTTACTAGAAAAGA..GAACTTTGAGAAGATACCTGGAAGTCGGAtt                                                                      >  24:16585/1‑100
                            aTAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATa                                                                      >  23:1225480/1‑101
                              aaGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGc                                                                    >  13:583657/1‑101
                              aaGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGc                                                                    >  22:386985/1‑101
                              aaGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGc                                                                    >  20:1144832/1‑101
                                ggAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCtt                                                                  <  2:591809/101‑1
                                ggAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCtt                                                                  >  5:745962/1‑101
                                 gAAAAAGCTTATCTCGCGGGCGGTAAAGCTCCAAATCTTAAGTATGTTCATGCCATTCAACAATC..GAATTTAGCAAAGATACCTGGCAGTCCTATAGCttt                                                                 >  23:1141581/1‑100
                                  aaaaaGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTAt                                                                >  24:1245887/1‑101
                                  aaaaaGCTTACATTGCTACCTTTAAAGATCCAAATCTTAAGTATGTTTCTCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGGCCGATAGCTTAt                                                                <  1:499301/101‑1
                                   aaaaGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTCAGTATGTTTATCGTAATAAACAAAC..GAACTTGGAGAAGATACCTGGAAGTCCGATAGCTTAtt                                                               >  21:533575/1‑101
                                    aaaGCTTATCTTGCTGCCGTTAAAGATCTAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTg                                                              >  10:1092510/1‑101
                                    aaaGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTg                                                              >  5:399581/1‑101
                                     aaGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATCCCTGGAAGTCCGATAGCTTATTgg                                                             <  5:946683/101‑1
                                     aaGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTgg                                                             >  28:47810/1‑101
                                           aTCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGt                                                       >  12:1197591/1‑101
                                             cTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTaa                                                     >  8:1201772/1‑101
                                             cTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTaa                                                     >  28:132623/1‑101
                                             cTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTaa                                                     <  27:104078/101‑1
                                                  tgcCGTTAAAGATCCAAATCTTCAGTATGTTTATCGTACTAAACAAAC..GACCTTTGAGGAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAAcccc                                                <  14:942915/101‑1
                                                   gcCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCt                                               <  11:56128/101‑1
                                                    ccGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCtt                                              <  8:949739/101‑1
                                                     cGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTa                                             >  29:618605/1‑101
                                                      gTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCTATAGCTTATTGGGCGAGTAAACCCCTTAt                                            >  22:1059699/1‑101
                                                      gTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTAt                                            <  13:243178/101‑1
                                                      gTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTAt                                            >  12:7028/1‑101
                                                      gTTAAAGACACAAGTCGTAAGTATGTATACTGTACTAA.CA.GCATGAATTTTGAAAAGATCCGTGGAAGTCCGATAGCTTATTGGGCGAGTAAAACCCTTAt                                            <  11:817825/101‑1
                                                       ttAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GACCTTTCAGAAGCTACCTGGAAGTCAGACATCTTATTTGGCGAGTAAAGCACTTAtt                                           >  21:1032666/1‑101
                                                          aaGATCCAAATCTTAGGTAGTTCGATCGTACTAATAAAAC..GAACTTTGAGAAAAGACCTGGGAGCCCGATAGCTTGTTGGGCGAGTAAACCCCTTATTTcc                                        <  13:996131/101‑1
                                                           aGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCg                                       <  1:1223950/101‑1
                                                            gATCCAAATCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGa                                      >  17:1167011/1‑101
                                                                 aaaTCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTg                                 >  14:11040/1‑101
                                                                  aaTCTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGCTAGCTGATGGGTCGAGTAAACCGCTTCTTACCAGTTTTGa                                >  5:234822/1‑101
                                                                  aaTCTTAAGTATGTTTATCGTAAAAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGTTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGa                                <  10:985140/101‑1
                                                                     cTTAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAAt                             <  8:1137789/101‑1
                                                                     cTTAAGTATGTTTATCGTACTAAAAAAAA..GAACTTTGAGAAGATACCTGGAAGTCCGATGGTTTATTGTGAGAGCAAACAACTCCATTCAGATTTTGTAAt                             >  25:718041/1‑101
                                                                      ttAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTTGAAGTCCTATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAAtt                            <  8:539536/101‑1
                                                                      ttAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAAtt                            <  10:842917/101‑1
                                                                       tAAGTATGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTg                           <  12:258146/101‑1
                                                                        aaGTATGTTTATCGTACTAAACAAAC..CAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGGAAACCCCCTATTTCGGCTCTTGCAATTgc                          >  23:344253/1‑100
                                                                            aTGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGATGTCCGATAGCTTATTGGGTGATTAAATGCCTTCTTTCCGCTTTTGACATTGGCATa                      >  16:1130266/1‑101
                                                                             tGTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATAc                     <  6:746534/101‑1
                                                                              gTTTATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATAcc                    >  10:180058/1‑101
                                                                               tttATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCa                   >  16:881899/1‑101
                                                                               tttATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCa                   >  8:301558/1‑101
                                                                               tttATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTACGATTTTGAAATTGGCATACCa                   <  3:846180/101‑1
                                                                               tttATCGTACTAAACAAAC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTACGATTTTGAAATTGGCATACCa                   <  18:846180/101‑1
                                                                                 tATCGTACTAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACt                 <  29:1019981/101‑1
                                                                                 tATCGTACTAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATCGGCATACCACt                 <  29:1203096/101‑1
                                                                                  aTCGTACTAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTc                <  13:448387/101‑1
                                                                                  aTCGTACTAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTc                >  7:278756/1‑101
                                                                                      tACTAAACAAGC..GAACTTTGAGAAGATACCTCGAATGCCGATAGCTTATGGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAg            <  25:113105/101‑1
                                                                                        cTAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGAc          >  3:154808/1‑101
                                                                                        cTAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGAc          >  18:154808/1‑101
                                                                                         tAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGCTAGCTTATTGGGCGAGTAAAACCCTTATTTGCCATTTTGAAATTGGCATACAACTCAAAGACt         >  28:632362/1‑101
                                                                                         tAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACt         <  26:658830/101‑1
                                                                                          aaaCAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACtt        >  23:829509/1‑101
                                                                                          aaaCAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACtt        <  19:984282/101‑1
                                                                                          aaaCAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACtt        <  4:984282/101‑1
                                                                                           aaCAAGC..TAACTTTGAGAAGAGACCTGTAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGACTTTGAAATTGGCATACCACTCAAAGACTTa       <  24:551600/101‑1
                                                                                            aCAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAg      <  19:898000/101‑1
                                                                                            aCAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAg      <  4:898000/101‑1
                                                                                             cAAGC..CAACTTCGAGAAGATACATGTAAGTGCGATAGCTTATTGGGCAAGTAAACCCCTGATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAGt     <  3:946830/101‑1
                                                                                             cAAGC..CAACTTCGAGAAGATACATGTAAGTGCGATAGCTTATTGGGCAAGTAAACCCCTGATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAGt     <  18:946830/101‑1
                                                                                              aaGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGGTTATTGGGCGAGTAAACCCCTTATTTCCAACGTTGAAATTGGCAGACGACTCACAGACTTCGtt    >  20:777097/1‑101
                                                                                              aaGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGCTTTCGAAATTGGCATACCACTCAACGACTCAGtt    >  9:530878/1‑101
                                                                                              aaGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAGtt    <  19:674463/101‑1
                                                                                              aaGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAGtt    <  6:244404/101‑1
                                                                                              aaGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAGtt    >  25:859/1‑101
                                                                                              aaGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAGtt    <  4:674463/101‑1
                                                                                               aGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAGTTg   >  25:449683/1‑101
                                                                                                |                                                                                                      
CGGCTGGTCGATTTTGATTCACAGGCACATAAGGAAAAAGCTTATCTTGCTGCCGTTAAAGATCCAAATCTTAAGTATGTTTATCGTACTAAACAAGC..GAACTTTGAGAAGATACCTGGAAGTCCGATAGCTTATTGGGCGAGTAAACCCCTTATTTCCGATTTTGAAATTGGCATACCACTCAAAGACTTAGTTGA  >  NC_013198/2159914‑2160110

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 28 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG