Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 33.7% intergenic (‑195/+113) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC33.7% 255.1 / 176.3 181intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (59/61);  new base C (57/4);  total (116/65)
Fisher's exact test for biased strand distribution p-value = 4.76e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-02

Only 100 of 183 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084
                                                                                                   |                                                                                                         
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACACAGTGCGCTCTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                          >  1:2775941/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                          >  1:1070749/1‑100 (MQ=255)
  ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                         >  2:525132/1‑100 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                    >  2:2292391/1‑100 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  2:1348758/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  2:575896/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  2:1673983/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                   <  1:2577210/100‑2 (MQ=37)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                  <  2:1617370/100‑3 (MQ=37)
          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                 <  2:3350030/100‑1 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                              >  2:1976710/1‑100 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaa                                                                                               <  1:973800/100‑3 (MQ=25)
             cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                             >  2:433335/1‑100 (MQ=255)
              gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                            >  2:453470/1‑100 (MQ=255)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                          <  1:481235/100‑1 (MQ=255)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                         >  2:2647779/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTACCCGGCCTAC‑AAAATCGTg                                                                                      <  1:490811/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                     <  1:1576694/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                     <  2:303189/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                     <  1:1246125/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                    <  1:2412074/100‑1 (MQ=255)
                         agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                                 <  2:1376918/100‑1 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                              >  1:1514883/1‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                              <  2:1070752/100‑1 (MQ=255)
                             acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                                             >  2:652931/1‑100 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                            <  1:2673634/100‑1 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                        >  2:2882506/1‑100 (MQ=255)
                                  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                        <  2:2850252/100‑1 (MQ=255)
                                  gcCACCGGGTATACAGATTCTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                        <  1:360147/100‑1 (MQ=255)
                                  cccACCGGGTATACAGATTGATGACAAAATGCGCCATGTTTATGACAGAAACACCGTTAACGACTTGTCCGTCCTAC‑AACATAATGCAAATTCAAGACa                                                                         >  1:2873693/2‑99 (MQ=255)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCGTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                       >  2:2377219/1‑100 (MQ=255)
                                      ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                    <  1:1618408/100‑1 (MQ=255)
                                        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                                 >  1:417066/1‑100 (MQ=255)
                                             taCAGATTGTTCACAATGTTCGCTTTGTTTGTTCCGGATGCGGCGTTTAGCCCTTATCCGGTCTAC‑AGTATCGTTCAATTTCATTTCATTGAA‑GA‑Gt‑tt                                                           >  1:2597208/1‑99 (MQ=12)
                                              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTc                                                         <  2:1808234/100‑1 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑Gt                                                     <  2:92022/100‑1 (MQ=255)
                                                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTa                                                    <  1:1000320/100‑1 (MQ=255)
                                                  aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAg                                                   >  2:2161586/1‑100 (MQ=255)
                                                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTa                                                    >  1:1692304/1‑98 (MQ=255)
                                                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAgg                                                  <  2:1223669/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGcc                                                >  1:1379029/1‑100 (MQ=255)
                                                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCt                                               <  2:3246252/100‑1 (MQ=255)
                                                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCt                                               <  1:2990399/100‑1 (MQ=255)
                                                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCt                                               >  1:3114353/1‑98 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATa                                           <  1:73175/100‑1 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATa                                           <  1:1185480/100‑1 (MQ=255)
                                                            aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAg                                         >  2:2525288/1‑100 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                                      <  1:1464915/100‑1 (MQ=255)
                                                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                                     <  1:2525268/100‑1 (MQ=255)
                                                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                    >  2:1512176/1‑100 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTCC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                   <  1:2167561/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                   <  2:1545345/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                   <  1:156416/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                   <  1:859607/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  2:1265712/100‑1 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  2:2254074/100‑1 (MQ=255)
                                                                     ctttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAgcgcagcgca                                >  2:3390066/9‑100 (MQ=21)
                                                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                               >  1:2121174/1‑100 (MQ=255)
                                                                         gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:1291721/5‑100 (MQ=18)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:3358721/5‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:3149123/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:2382679/3‑100 (MQ=21)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:3132533/3‑100 (MQ=21)
                                                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                          <  2:209180/100‑1 (MQ=255)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:1085467/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:2398409/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCCTTGCGAATTCAATATATTGTATGG‑AA‑ATTC‑‑‑‑AGGCCTGATAAGCGTAGCGCATCAGGctt                        >  2:1840620/2‑98 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2699243/2‑100 (MQ=18)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:879412/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:247189/1‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:1782752/1‑99 (MQ=17)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:3018624/1‑100 (MQ=21)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  2:546887/9‑98 (MQ=18)
                                                                              ccggataaggcgTTCACGCCTCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:2072701/9‑98 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                       >  2:224570/9‑98 (MQ=25)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  2:781015/9‑98 (MQ=21)
                                                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatt                      >  2:2618651/8‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1897287/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1667946/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:671829/7‑100 (MQ=14)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CCAC‑‑GTAGGCCTTATAAGCGTAGCGCATCAGGCAgttt                     >  1:2621897/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattc                     >  1:2426462/7‑99 (MQ=21)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                    >  1:3025458/7‑98 (MQ=21)
                                                                                ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCA‑GG‑AA‑CGAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                     >  2:569447/7‑96 (MQ=25)
                                                                                 gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  2:3318869/6‑100 (MQ=18)
                                                                                 gataaggcgTTCACGCCTTATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GA‑GA‑TGGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttttt                    >  2:994446/6‑100 (MQ=14)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:2500960/100‑1 (MQ=255)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:3231696/100‑1 (MQ=255)
                                                                                   taaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAAGCAgttt                     >  2:779446/4‑97 (MQ=21)
                                                                                     cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  2:852676/1‑100 (MQ=255)
                                                                                      ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               >  2:2436216/1‑100 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:318497/100‑1 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  2:1308958/100‑1 (MQ=255)
                                                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTCTTTGTCa          >  1:826406/1‑100 (MQ=255)
                                                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  1:1605822/1‑100 (MQ=255)
                                                                                            aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTTGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt         >  1:364156/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:1691979/100‑1 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  2:1653727/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:3286066/100‑1 (MQ=255)
                                                                                                ccTTATCCGGCCTAC‑AAAATCGTACAAATTCAATACATTGCA‑GA‑CT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATGTTTGCCTTTGCCCTCAGt     >  1:3326085/1‑100 (MQ=21)
                                                                                                   |                                                                                                         
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑