Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 775,620 C→G 5.0% L156L (CTC→CTG tolQ → membrane spanning protein in TolA‑TolQ‑TolR complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913775,6200CG 5.0% 315.6 / 6.7 100L156L (CTC→CTGtolQmembrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base C (47/47);  new base G (2/3);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

Only 100 of 101 total aligned reads displayed.

GGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTA  >  NC_000913/775472‑775769
                                                                                                                                                    |                                                                                                                                                     
gggCGTCGCGTGCTATGCGTATCTCCATGACCCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCGGGGGGATCATGCACGCCTTTATCGCCCTCg                                                                                                                                                      >  1:130223/1‑150 (MQ=255)
    gTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGc                                                                                                                                                  <  1:989828/150‑1 (MQ=255)
    gTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGc                                                                                                                                                  <  1:1393923/150‑1 (MQ=255)
       gcgTGCTATGCGTATCTCCATGAACCGGGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGt                                                                                                                                               <  1:576897/150‑1 (MQ=255)
        ggTGCTATGCGTATCTGCATGAACCGTGATCTTGAAAATCTGGAAACGCAGATTCCGTTCCGCGGTACGGGTGGCTCCCTCACCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTa                                                                                                                                              <  1:897335/149‑1 (MQ=255)
         ttGCTATGCGTATCTCCCTGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTaa                                                                                                                                             <  1:1674887/149‑1 (MQ=255)
           gCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGGATATTGGTCTGGTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTaaaa                                                                                                                                           >  1:1603733/1‑150 (MQ=255)
             tATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAACCGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACa                                                                                                                                         <  2:592572/150‑1 (MQ=255)
             tATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACa                                                                                                                                         <  2:546992/150‑1 (MQ=255)
              aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACaa                                                                                                                                        >  2:129354/1‑150 (MQ=255)
              aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTACCGCCCCGAGGGCGATaaaggcg                                                                                                                                        <  1:1410975/150‑5 (MQ=255)
              aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTCTGGTACGGTCTGGGGGATCATGCACGCCTTTACCGCCCCGAGGGCGATaaaggcg                                                                                                                                        <  1:1449618/150‑5 (MQ=255)
               tGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAg                                                                                                                                       >  2:395692/1‑150 (MQ=255)
               tGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAg                                                                                                                                       <  1:195809/150‑1 (MQ=255)
                gCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGc                                                                                                                                      <  2:1649605/150‑1 (MQ=255)
                gCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGc                                                                                                                                      <  2:1364818/150‑1 (MQ=255)
                  gTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGc                                                                                                                                                  <  1:947065/136‑1 (MQ=255)
                   tATCTCCATGATCCGTGATCTTGAAAATCTGGAAACGCACATTCCGGTCCTGGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGGTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAac                                                                                                                                   <  2:662946/150‑1 (MQ=255)
                   tATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAac                                                                                                                                   <  1:1480412/150‑1 (MQ=255)
                   tATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGGCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAac                                                                                                                                   <  2:1325437/150‑1 (MQ=255)
                    aTCTCCATGAACCGTGAAATTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTCGGACGGGCTGGGGGATCATGCACGCCTTGATCGCCGTCGGGGGGGTAAGACAAggacca                                                                                                                                  >  1:289406/1‑145 (MQ=255)
                       tcCATGAACCGTGAACTTGAAAATCTGGTAACGCACATTCCGGTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGGTTGCTACGGGCTGGGGGATCATGCACGCCTTTATCGCCCTAGGGGCGGTAAAACAATCATCACtc                                                                                                                               >  1:250324/1‑149 (MQ=255)
                       tcCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTg                                                                                                                               >  1:319662/1‑150 (MQ=255)
                        ccATGAACCGTGAACTTGAAAATCTGGAAACGCACCTTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACTGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGc                                                                                                                              <  1:1542114/150‑1 (MQ=255)
                              aCCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATgg                                                                                                                        >  1:1528764/1‑150 (MQ=255)
                               ccGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGt                                                                                                                       >  1:991382/1‑150 (MQ=255)
                                 gTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTg                                                                                                                     <  2:1133357/150‑1 (MQ=255)
                                 gTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTg                                                                                                                     <  2:319662/150‑1 (MQ=255)
                                    aaCTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTgcgc                                                                                                                  >  1:1038032/1‑150 (MQ=255)
                                      cTTGAAAATCTGGAACCGCACATTCCGTTCCTAGGGACGGTTGGCTCCCTCAGCCCTTATATTGGTCTGTTTGGTACGGTATGGGGGATCGAGCACGACTTTTTCGCACGCGGGGCGGGAACACAAGCAACACGGCAAATGGTTgcgcac                                                                                                                >  2:1141754/1‑148 (MQ=255)
                                       ttGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCg                                                                                                               <  1:1016380/150‑1 (MQ=255)
                                               cTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCa                                                                                                       <  1:259515/150‑1 (MQ=255)
                                                     aCGCACATTCCGTTCCTCGGGACGGTTGGCTCCATCAGCCAGTATATTGGTCTGTTTGGTACGGGCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCg                                                                                                 <  1:1157967/150‑1 (MQ=255)
                                                      cGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGt                                                                                                >  2:664893/1‑150 (MQ=255)
                                                      cGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGt                                                                                                >  1:1533550/1‑150 (MQ=255)
                                                       gCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGtt                                                                                               <  2:425989/150‑1 (MQ=255)
                                                        cacaTTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATAGTGGTCTGTTTGGTACTGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTg                                                                                              <  2:1608953/150‑1 (MQ=255)
                                                         acaTTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGa                                                                                             >  1:871766/1‑150 (MQ=255)
                                                             tCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGc                                                                                         <  1:996879/150‑1 (MQ=255)
                                                               cGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTCTCGCAGAAGCGTTGATTGCgc                                                                                       >  2:1522027/1‑149 (MQ=255)
                                                                gTTCCTCGGTACGGGTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGAc                                                                                      <  1:689170/150‑1 (MQ=255)
                                                                gTTCCGCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGAc                                                                                      <  1:986649/150‑1 (MQ=255)
                                                                  tcatCGGGACGGCTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATTCACTCCGTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTg                                                                                    <  1:1483803/147‑1 (MQ=255)
                                                                   ccTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGc                                                                                   >  1:808403/1‑150 (MQ=255)
                                                                   ccTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCGTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGc                                                                                   <  2:1717666/150‑1 (MQ=255)
                                                                      cGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAAt                                                                                >  1:1257754/1‑150 (MQ=255)
                                                                         tACGGGTGGCTCCCTCAGCCCGTATATTGGTCTGTTTGGTACGGTCGGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTgg                                                                             <  1:692330/150‑1 (MQ=255)
                                                                              ttGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGTGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGt                                                                        <  2:935123/150‑1 (MQ=255)
                                                                               gggCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGtt                                                                       <  2:383141/149‑1 (MQ=255)
                                                                                 gCTCCATCAGCCCGTTTTTTGGTGTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAATACTAGCAACTCTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTg                                                                     <  1:1588508/150‑1 (MQ=255)
                                                                                   tCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATAGCCCTCGGGGCGGTAAAACAAGCACCACTGCAAATGGTTGCGCCCGGTATCGCAGCAGCGGTGAATGCGACTGCAATTGGTATGTTAGcc                                                                   >  1:135951/1‑150 (MQ=255)
                                                                                   gccATCAGCCCGTATATTTGTCTGTTTGGTACGGTCGGGGGGATCATGCACGCCTTTATCGCCCGCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGcc                                                                   <  2:773452/149‑1 (MQ=255)
                                                                                           cccGTATATTGGTCTGGTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCgggatagc                                                           <  1:918659/150‑8 (MQ=255)
                                                                                                  atTGGTGTGTTTGGTACGGTCTGGGGGATCATGCACGCGTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGGTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGtt                                                    <  2:176707/150‑1 (MQ=255)
                                                                                                  atTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGtt                                                    <  2:794491/150‑1 (MQ=255)
                                                                                                    tGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTAt                                                  >  2:607175/1‑150 (MQ=255)
                                                                                                     ggTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCATATGGTTGCGCCGGGTATCGCTGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTTTg                                                 >  1:1496876/1‑150 (MQ=255)
                                                                                                      gtttGTTTGGTACGGGCTGGGGGAGCATGCACGCATTTCTCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATgg                                                <  2:1083138/147‑1 (MQ=255)
                                                                                                      gTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGTAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATgg                                                >  1:998523/1‑150 (MQ=255)
                                                                                                      gTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATgg                                                <  2:371568/150‑1 (MQ=255)
                                                                                                        cTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCATATGGTTTAGCCCTGTATCGCAGAAGAGTTGATTGCGACTGAAATTGGTCTGTTTGCCGCTAGCCCGGCAGTTATGGcc                                              >  1:1062549/1‑150 (MQ=255)
                                                                                                          gTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATGCCGGCAGTTATGGCCTa                                            >  2:739627/1‑150 (MQ=255)
                                                                                                          gTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTa                                            >  2:739628/1‑150 (MQ=255)
                                                                                                          gTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTa                                            >  2:1250849/1‑150 (MQ=255)
                                                                                                          gTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTa                                            >  1:1067422/1‑150 (MQ=255)
                                                                                                              ggTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAAc                                        >  2:7008/1‑150 (MQ=255)
                                                                                                                tACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGCCTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCg                                      >  1:1556965/1‑150 (MQ=255)
                                                                                                                tACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCg                                      >  1:63849/1‑150 (MQ=255)
                                                                                                                  cGGGCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACGGCAAATGGGTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGcc                                    <  2:1045665/150‑1 (MQ=255)
                                                                                                                   ggTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCt                                   <  1:395692/150‑1 (MQ=255)
                                                                                                                     tCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCa                                 <  2:1528764/150‑1 (MQ=255)
                                                                                                                     gcTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAATCAACACTGCAAATGGTTGCTCCCGGTATCGCAGATGCGTTGATTGCGAATGCAATTGGTCTTTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCa                                 <  1:441381/149‑1 (MQ=255)
                                                                                                                      cTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCaa                                >  1:552502/1‑150 (MQ=255)
                                                                                                                        gggggATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAAcc                              >  2:476162/1‑150 (MQ=255)
                                                                                                                        gggggATCATGAACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAAcc                              >  2:1637655/1‑150 (MQ=255)
                                                                                                                           ggATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCATCACTGCAAATGGTTGCGCCCGGTATCTCAGAAGCGTTGATTGGGACTGCATTTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTTCAACAGCCTCAAACAgg                           >  1:321934/1‑149 (MQ=255)
                                                                                                                           ggATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAgc                           >  2:1716602/1‑150 (MQ=255)
                                                                                                                           ggATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAgc                           >  1:1215663/1‑150 (MQ=255)
                                                                                                                               cATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGt                       <  2:1178573/150‑1 (MQ=255)
                                                                                                                                  tgcttgttttACCGCCCCGAGGGCGATAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaa                    >  1:660641/8‑150 (MQ=255)
                                                                                                                                  gCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaa                    >  2:1031173/1‑150 (MQ=255)
                                                                                                                                   cACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCACCCAGCGCGTaaac                   >  2:1536892/1‑150 (MQ=255)
                                                                                                                                   cACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaac                   >  2:213228/1‑150 (MQ=255)
                                                                                                                                    cttgttttACCGCCCCGAGGGCGATGAAGCAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTa                                                               <  1:1533152/100‑1 (MQ=255)
                                                                                                                                    cttgttttACCGCCCCGAGGGCGATGAAGCAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTa                                                               >  2:1533152/6‑105 (MQ=255)
                                                                                                                                    aCGACTTTATAGCACTCGGGGCGGTAAAACAAGCACCACTGCATATGTTTGCGCCAGGTATCGCAGAAGCGTTGATTGAGACTGCAATTGGTCTGTTTGACGATATCACGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaaca                  <  1:318848/150‑1 (MQ=255)
                                                                                                                                      gCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaacaaa                >  2:1261674/1‑150 (MQ=255)
                                                                                                                                       ccTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaacaaac               >  1:1722946/1‑150 (MQ=255)
                                                                                                                                       ccTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaacaaac               >  1:1137397/1‑150 (MQ=255)
                                                                                                                                        cTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTCTCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACt              <  1:430878/150‑1 (MQ=255)
                                                                                                                                        cTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACt              >  2:1039740/1‑150 (MQ=255)
                                                                                                                                         tttATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTg             <  1:1004687/150‑1 (MQ=255)
                                                                                                                                          ttATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTgg            >  1:1051857/1‑150 (MQ=255)
                                                                                                                                          atATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCAGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTCTCCCGGCAGTTTAGGCCTACAACCGACTCATCCAGCGCGTAAACAAACTgg            >  1:1521621/2‑150 (MQ=255)
                                                                                                                                           tATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGTCTGCAATTGGTCTGTTTGCCGCTATCACGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGa           >  1:1616712/1‑150 (MQ=255)
                                                                                                                                             tCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAAc         >  1:1521608/1‑150 (MQ=255)
                                                                                                                                              cGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACt        <  2:1611968/150‑1 (MQ=255)
                                                                                                                                               gCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGACGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTg       >  2:848279/1‑150 (MQ=255)
                                                                                                                                                cccTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGa      <  2:552502/150‑1 (MQ=255)
                                                                                                                                                    cGGGGCGGTAAAACACGCAACCCTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTCCTATTGGTCTGATTGCCCCTATCCCGGCAGTTATGGCCTACAACCGCGGCGACCAGCGCGTCAACAAACTGGAACGGaagtc  >  2:886865/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTA  >  NC_000913/775472‑775769

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