Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 19.9% 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,9850TC19.9% 423.0 / 93.3 182intergenic (‑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 (63/82);  new base C (34/2);  total (98/84)
Fisher's exact test for biased strand distribution p-value = 9.69e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

Only 100 of 187 total aligned reads displayed.

TGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705889‑3706084
                                                                                                |                                                                                                          
tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCAGATGCGGCGTTAACGCCTTATc                                                                                                         <  1:1338895/100‑1 (MQ=255)
 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        >  2:3239895/1‑100 (MQ=255)
 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        <  2:2684741/100‑1 (MQ=255)
     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  2:3080071/100‑1 (MQ=255)
     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  1:142860/100‑1 (MQ=255)
      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  1:1020673/100‑1 (MQ=255)
      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  2:2946884/100‑1 (MQ=255)
      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                   <  1:1487344/100‑3 (MQ=37)
         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                               <  1:984224/100‑1 (MQ=255)
              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                          >  2:2689433/1‑100 (MQ=255)
               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                         <  2:1273863/100‑1 (MQ=255)
                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                      <  2:1839423/100‑1 (MQ=255)
                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                      >  1:2278993/1‑100 (MQ=255)
                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                      <  1:2258861/100‑1 (MQ=255)
                     gagGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAATGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                                   >  1:2007920/1‑100 (MQ=255)
                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                                <  1:3087489/100‑1 (MQ=255)
                         cacaGCGCCACCGGGTATACAGATTGTTGACATAGTGCGGTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                               <  1:2346787/100‑1 (MQ=255)
                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                               <  2:1583771/100‑1 (MQ=255)
                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                               <  1:3320623/100‑1 (MQ=255)
                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                               >  1:3197023/1‑100 (MQ=255)
                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                             >  1:956647/1‑100 (MQ=255)
                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                         >  2:139067/1‑100 (MQ=255)
                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                        >  1:644144/1‑100 (MQ=255)
                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                      <  2:1691793/100‑1 (MQ=255)
                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                      <  2:1203671/100‑1 (MQ=255)
                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                      <  2:1153347/100‑1 (MQ=255)
                                       gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑At                                                             <  2:3143498/100‑1 (MQ=255)
                                        tatataGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAATTCGTGCAAATTCAATACATTGC‑‑A‑GA‑At‑t                                                           <  2:1898055/100‑1 (MQ=255)
                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTc                                                       <  2:535277/100‑1 (MQ=255)
                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTc                                                       <  1:1219074/100‑1 (MQ=255)
                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGTCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCg                                                      <  1:1827005/100‑1 (MQ=255)
                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCg                                                      <  2:196868/100‑1 (MQ=255)
                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCg                                                      <  2:1716538/100‑1 (MQ=255)
                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCg                                                      >  2:581642/1‑100 (MQ=255)
                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTa                                                    <  1:2238269/100‑1 (MQ=255)
                                               aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAg                                                   <  2:1527489/100‑1 (MQ=255)
                                                ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAgg                                                  <  2:1027915/100‑1 (MQ=255)
                                                gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATAGC‑‑A‑GA‑AT‑T‑TTCGTAgg                                                  >  1:1956472/7‑100 (MQ=18)
                                                 cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGc                                                 >  1:3043652/6‑100 (MQ=18)
                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGcc                                                >  1:2320545/1‑100 (MQ=255)
                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGcc                                                >  1:1108108/1‑100 (MQ=255)
                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCt                                               <  2:1467663/100‑1 (MQ=255)
                                                    tgACAGAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AGAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTg                                              <  1:1509216/100‑1 (MQ=255)
                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGa                                             >  1:2210422/1‑100 (MQ=255)
                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGAt                                            <  2:160962/100‑1 (MQ=255)
                                                         aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAg                                         >  2:2904942/1‑100 (MQ=255)
                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCg                                       <  2:1547524/100‑1 (MQ=255)
                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCgg                                      <  1:501890/100‑1 (MQ=255)
                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCgg                                      >  1:1679104/1‑100 (MQ=255)
                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAgcg                                  <  1:1908111/100‑1 (MQ=255)
                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAgcg                                  <  1:581642/100‑1 (MQ=255)
                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCa                                >  1:467457/1‑100 (MQ=255)
                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCa                                >  1:101451/1‑100 (MQ=255)
                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑A‑AT‑CT‑C‑CGTGTAGGCCTGATAAGCGTAGCGCAt                               >  1:2916067/8‑100 (MQ=18)
                                                                   cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑‑A‑TG‑AG‑C‑CACGTAGGCCTGATAAGCGTAGCGCAt                               >  2:2597289/8‑100 (MQ=17)
                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATc                              >  1:395965/1‑100 (MQ=255)
                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  2:1108111/100‑1 (MQ=255)
                                                                     tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGC‑‑A‑GG‑AT‑A‑CGTGTAGGCCTGATAAGCGTAGCGCATCa                             >  2:3311802/6‑100 (MQ=18)
                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg                            >  1:2736833/1‑100 (MQ=255)
                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAACACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg                            >  1:2145762/1‑100 (MQ=255)
                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCATGCAAATTCAATATATTGC‑‑A‑GA‑GT‑T‑CACGTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:815794/5‑100 (MQ=17)
                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑‑AAACCGCGCAAATTCAATATATTGCGGA‑GA‑AA‑A‑T‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  1:2491171/5‑100 (MQ=18)
                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCC‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:2566539/4‑100 (MQ=18)
                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCC‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2545684/3‑100 (MQ=18)
                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑‑AGGG‑AC‑T‑AT‑GTAGGCCTGATAAGCATAGCGCATCAGGc                          >  1:2687291/3‑100 (MQ=18)
                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGC‑‑A‑GAGAT‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2334551/2‑100 (MQ=17)
                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:2164984/1‑99 (MQ=18)
                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGC‑‑A‑GA‑AA‑C‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:182183/1‑99 (MQ=18)
                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑ATCC‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  1:1802484/1‑100 (MQ=17)
                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGA‑‑T‑CG‑AT‑C‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  2:1289661/1‑99 (MQ=21)
                                                                           ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑A‑GG‑AC‑C‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  1:2438299/9‑98 (MQ=18)
                                                                           ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  1:1042575/1‑100 (MQ=255)
                                                                             ngataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTGC‑‑A‑CA‑GC‑T‑AGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagatt                     >  1:124068/7‑96 (MQ=37)
                                                                             ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑‑GA‑AT‑CCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                    >  1:1090501/7‑98 (MQ=21)
                                                                             ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGC‑‑A‑GG‑AA‑A‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                     >  2:2270858/7‑100 (MQ=18)
                                                                              gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                    <  2:382968/100‑1 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttttn                   <  1:869271/100‑2 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:3018944/100‑1 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:139066/100‑1 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:2007934/100‑1 (MQ=255)
                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  2:2145778/100‑1 (MQ=255)
                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                <  1:3083147/100‑1 (MQ=255)
                                                                                      gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            >  2:1038959/1‑100 (MQ=255)
                                                                                       ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:1630361/100‑1 (MQ=255)
                                                                                        tAACGCCTTATCCGGCCTAC‑AACATCGTGCACATTCATTACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          <  1:2253522/100‑1 (MQ=255)
                                                                                         aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAACTTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCATAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  1:2800636/1‑99 (MQ=255)
                                                                                          aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATTCATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:1531702/100‑1 (MQ=255)
                                                                                          aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:700750/100‑1 (MQ=255)
                                                                                          aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:3057161/100‑1 (MQ=255)
                                                                                           cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       >  2:1990075/1‑100 (MQ=255)
                                                                                           cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAGCa       >  2:511778/1‑100 (MQ=255)
                                                                                            gCCTTATCCGGCCTAC‑AAAATCGTGTAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  2:1634866/100‑1 (MQ=255)
                                                                                            gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  1:2376126/100‑1 (MQ=255)
                                                                                             ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  2:766982/1‑100 (MQ=255)
                                                                                              cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc    >  2:2402876/1‑100 (MQ=255)
                                                                                               ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  2:2114964/100‑1 (MQ=255)
                                                                                               ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  2:3262999/100‑1 (MQ=255)
                                                                                                tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  2:1939554/1‑100 (MQ=255)
                                                                                                tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  1:2429667/1‑100 (MQ=255)
                                                                                                tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   >  2:656396/1‑99 (MQ=255)
                                                                                                |                                                                                                          
TGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑‑A‑GA‑AT‑T‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705889‑3706084

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