Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,469,041 A→G 5.6% S271S (AGT→AGC chiA ← periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,469,0410AG 5.6% 326.9 / 3.0 107S271S (AGT→AGCchiAperiplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base A (40/59);  new base G (3/3);  total (45/62)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

TCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGACT  >  NC_000913/3468892‑3469189
                                                                                                                                                     |                                                                                                                                                    
tCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTa                                                                                                                                                      >  1:793488/1‑150 (MQ=255)
 ccGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTTGCTGATTTCTGTTGCCGTTGCCGCGCCAACATAACGCCATGGGTTAc                                                                                                                                                     >  2:1379206/1‑150 (MQ=255)
           ggTCACCGGTGTTGCCGGAGTTGGATCAACCGGATTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGtt                                                                                                                                           <  1:383126/150‑1 (MQ=255)
            gTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGttt                                                                                                                                          >  1:103512/1‑150 (MQ=255)
             tCACCGGTGTCTCCGGAGTCGGAGCCACCGGAGTCACAGCACCGCCATTATTATACTGCTTAACGGATCAGGAGCCATGGTTCCCAAACTGGCTGATTTGTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGGTAGTGGCATCGtttt                                                                                                                                         <  2:1638412/150‑1 (MQ=255)
              cACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttc                                                                                                                                        >  2:1831221/1‑150 (MQ=255)
              cACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttc                                                                                                                                        >  1:724334/1‑150 (MQ=255)
                cctgtgtTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTGGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGTCATCGTtttctt                                                                                                                                      <  1:1331249/147‑1 (MQ=255)
                  ggTGTTTCCGGTGTTGTATCAACCGGAGTCACAGCTCCGCCATTATTATCCGGCTTAACGGAACAGGAGCTAGGGTTACCATTCTGGCTGATTTCTGTTGCCGTTTCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctttc                                                                                                                                    <  1:1524573/150‑1 (MQ=255)
                  ggTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTtttctttc                                                                                                                                    <  1:1168432/150‑1 (MQ=255)
                  ggTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGGtttctttc                                                                                                                                    >  2:1210062/1‑150 (MQ=255)
                   gtgtTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCg                                                                                                                                   <  2:1979173/150‑1 (MQ=255)
                    tgtTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTCTCTTTCGc                                                                                                                                  <  2:25028/150‑1 (MQ=255)
                    tgtTTCCGGAGTTGGAGCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGc                                                                                                                                  <  2:1190799/150‑1 (MQ=255)
                      tttCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCAt                                                                                                                                >  1:1011582/1‑150 (MQ=255)
                       ttCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATc                                                                                                                               >  1:754284/1‑150 (MQ=255)
                       ttCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCACGCACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATc                                                                                                                               >  1:1602125/1‑150 (MQ=255)
                       tgccGGTGTTGGATCCTCCGGAGGCACCTCAGCGCCATTATTCTGCGGCGAAACCGAACAGGACCCTGGTTTACCCTACTGGCTGATTTCTGTTGCCGTTGCGGCGCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATc                                                                                                                               <  2:1511669/148‑1 (MQ=12)
                        tCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTTTCCGGCTTAACGGTACAGGAGCCAGGGTTACCATACTGGCTGATGTCTGTTGCCGTTGCCGCGCGCACATTACGCCATGGGTTACTGGCATCGTTTTCTTTCGCAtac                                                                                                                              >  2:1791931/1‑148 (MQ=255)
                         ccGGAGTTGGATCATCCGGAGTCACAGCTCCGCCATTATAATCCGGCTTAACGGAACAGGAGCCAGGGGTACCATACTGGCTGATTTCTGTTGCCGTTGCAGCGCGACCATAACGCCATGGGTTTCTGGCATCGTTGAGTTTCGCATAGc                                                                                                                             >  1:1061253/1‑150 (MQ=255)
                         ccGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGAGTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATAGTTTTCTTTCGCATCGc                                                                                                                             >  2:1234180/1‑150 (MQ=255)
                          cGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGcc                                                                                                                            >  1:765112/1‑150 (MQ=255)
                                tGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCa                                                                                                                      >  2:176234/1‑150 (MQ=255)
                                tGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCa                                                                                                                      >  1:1657164/1‑150 (MQ=255)
                                 ggATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCATTGGATTTCGCATCCCCCGGCCag                                                                                                                     >  2:370243/1‑149 (MQ=37)
                                 ggATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCaa                                                                                                                     >  1:249407/1‑150 (MQ=255)
                                  gATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCTGGCATCATTTGGTTTCGCATCTCCAGGTCatt                                                                                                                    >  1:773857/1‑148 (MQ=38)
                                  gATCAACCGGAGTCACAGCACCGCCATTATTATACGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTTCTGGCATCATTGATTTTCGCATCCCCTGGCCAAt                                                                                                                    >  1:1585013/1‑150 (MQ=37)
                                   atcatcCGTAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCTTGGTTACCATACTGGCTGATTTATGTTGACGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAAtt                                                                                                                   <  2:388214/150‑1 (MQ=255)
                                    tCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAAttt                                                                                                                  >  1:1279782/1‑150 (MQ=255)
                                      aaCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATACCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGc                                                                                                                >  1:847856/1‑150 (MQ=255)
                                      aaCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGc                                                                                                                <  1:1437416/150‑1 (MQ=255)
                                      aaCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGc                                                                                                                <  1:1934764/150‑1 (MQ=255)
                                        ccGGAGTCACATCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCTCGAACATAACGCCATGGGTTGCTGGCATCATTGCTTTTCGCATCACCCGGACAGttagaag                                                                                                              <  1:1367659/150‑6 (MQ=18)
                                        ccGGAGTCACACCACCGCCATTATTATCAGCCTTAACGGAACAGGAGCCAGGCATACAATCCTGGCAGAATTCTGTTGCCGTTGCCGCGCGAACATAACGCCATCGGTTACTGGCATGCTTATCATACGCATCGCCTGGGCAATTTGCcg                                                                                                              >  2:404056/1‑150 (MQ=255)
                                          ggAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCcgcg                                                                                                            >  2:665415/1‑150 (MQ=255)
                                            aGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGGTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGcc                                                                                                          <  1:535653/150‑1 (MQ=255)
                                               cacaGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAc                                                                                                       <  1:1860696/150‑1 (MQ=255)
                                               cacaGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAc                                                                                                       >  1:1785186/1‑150 (MQ=255)
                                               cacaGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAc                                                                                                       >  2:1594150/1‑150 (MQ=255)
                                                acaGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTAcc                                                                                                      <  1:1153831/150‑1 (MQ=255)
                                                   gCACCGCCAGTATTATCCGGCTTAACGGACCAGGAGCCAGGGTTACCTTACTGGCTGATTTCTGTTGCCGTTGCAGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTGCGCATCGCATCGGCAATTTGCCGCGCCTACccac                                                                                                   >  2:308366/1‑150 (MQ=255)
                                                      ccGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAg                                                                                                >  1:953726/1‑150 (MQ=255)
                                                        gcgattattATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCCCCAGGc                                                                                              >  2:1447974/4‑150 (MQ=255)
                                                            ttattaTCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTc                                                                                          >  2:1831837/1‑150 (MQ=255)
                                                             tattttTCGTCTAAACTGAACAGGAGCCAGGGTTTCCATTCTGGCTGATTTCTGTTGCCGTTGCCGCGCGTCCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                                                                         <  1:1510544/150‑1 (MQ=255)
                                                             tattctCCGGCTTAACGGAACAGCAGCCAGGGTGACCATACTGGCGGATTTCTGTTGCCGTTGCCGCGCGAACATCACGCCATGGGTTACTGGCAGCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                                                                         <  2:518000/150‑1 (MQ=255)
                                                             tattatCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                                                                         <  2:1611183/150‑1 (MQ=255)
                                                             tattatCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGGTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                                                                         <  2:1513938/150‑1 (MQ=255)
                                                             tattatCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCGGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCt                                                                                         <  2:1964201/150‑1 (MQ=255)
                                                                 atCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTtata                                                                                     <  2:180278/150‑1 (MQ=255)
                                                                 atCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTtata                                                                                     <  2:531464/150‑1 (MQ=255)
                                                                   ccGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTg                                                                                                                 >  2:1939452/1‑120 (MQ=255)
                                                                   ccGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTg                                                                                                                 <  1:1939453/120‑1 (MQ=255)
                                                                     gtcTTAACGGAACAGGAGCCAGGGTGACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                                                                 <  1:56940/148‑1 (MQ=255)
                                                                     ggCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                                                                 <  1:729745/150‑1 (MQ=255)
                                                                     ggCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                                                                 <  1:497647/150‑1 (MQ=255)
                                                                     ggCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATg                                                                                 <  1:2004624/150‑1 (MQ=255)
                                                                        ttAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgcct                                                                              <  2:349234/150‑2 (MQ=255)
                                                                          aaCGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATgccgcc                                                                            <  2:1895223/150‑1 (MQ=255)
                                                                                  aGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGa                                                                    <  2:1490706/150‑1 (MQ=255)
                                                                                   ggAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGAt                                                                   >  1:1974028/1‑150 (MQ=255)
                                                                                     aGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGTCGTTCTTATAAATGCCGCCGTTATAGATaa                                                                 >  2:1262186/1‑150 (MQ=255)
                                                                                      gCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAAc                                                                <  1:843689/150‑1 (MQ=255)
                                                                                         aGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTa                                                             <  2:793487/150‑1 (MQ=255)
                                                                                           ggTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTAcc                                                           >  1:1672416/1‑150 (MQ=255)
                                                                                            gTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCa                                                          <  2:754283/150‑1 (MQ=255)
                                                                                              tACCATACTGTCTGATTTCTGTTGACGTTGCCGCGCGCACATAACGCCATTGTTTACTGGCATCGTTTTCTTTCGCATCGCCTTGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATg                                                        <  1:999213/150‑1 (MQ=255)
                                                                                               ttcATACTTGCTGAATTCGGGTGCAGTCGCCTCGCGCACCTGACGGGATGGGATACTTGCTGCGTTCCGATCCGCACCCCCTGGGCAATTTTCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt                                                       <  1:1197774/148‑1 (MQ=255)
                                                                                               aCCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt                                                       <  1:1831222/150‑1 (MQ=255)
                                                                                               aCCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCACCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGt                                                       <  1:1292499/150‑1 (MQ=255)
                                                                                                  aTACTGGCTGATGTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGTGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTaga                                                    <  1:1253939/150‑1 (MQ=255)
                                                                                                   tACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTGCt                                                                                                                                                    <  2:137005/53‑1 (MQ=255)
                                                                                                    aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaag                                                  <  1:1676069/150‑1 (MQ=255)
                                                                                                    aCTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTagaag                                                  <  1:1343933/150‑1 (MQ=255)
                                                                                                      tGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa                                                <  2:1663180/150‑1 (MQ=255)
                                                                                                      tGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAgaagaa                                                <  1:176234/150‑1 (MQ=255)
                                                                                                       ggCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATCTATAACggcg                                                            >  2:1358718/1‑134 (MQ=255)
                                                                                                       ggCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATCTATAACggcg                                                            <  1:1358719/137‑4 (MQ=255)
                                                                                                        gCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAAcc                                              <  2:847855/150‑1 (MQ=255)
                                                                                                            aTTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTcc                                          <  2:1796073/150‑1 (MQ=255)
                                                                                                             tttCTGTTGCCGTTGCCGCGCGAACATAATGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCt                                         <  1:465223/150‑1 (MQ=255)
                                                                                                                          tgccgCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGc                            <  1:1594151/150‑1 (MQ=255)
                                                                                                                                 cgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCTCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACa                     <  2:1861317/150‑1 (MQ=255)
                                                                                                                                 cgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACa                     >  2:163123/1‑150 (MQ=255)
                                                                                                                                 cgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACa                     <  1:1398086/150‑1 (MQ=255)
                                                                                                                                 cgAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACa                     <  1:1831838/150‑1 (MQ=255)
                                                                                                                                    aCATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCa                  <  2:1785185/150‑1 (MQ=255)
                                                                                                                                     cATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCAt                 >  2:260550/1‑150 (MQ=255)
                                                                                                                                      aTAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATa                >  2:511769/1‑150 (MQ=255)
                                                                                                                                       tAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAg               >  2:890415/1‑150 (MQ=255)
                                                                                                                                        aaCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGt              <  1:599929/150‑1 (MQ=255)
                                                                                                                                         aCGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGACGCGCCTACCCACCAGGCATTGTTATATACCGCGCCTTTATAGACGACATACCAGGTATTACAGCCATCCTGAGTATTCCAGGCTACTACAGAATAGtt             >  2:931642/1‑150 (MQ=255)
                                                                                                                                           gCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATACCGTACCAAGTAGCAGAACCTTCCTGACCATTCCTGGCTACTACAGCATAGTTGt           >  1:91485/1‑150 (MQ=255)
                                                                                                                                               tGGGTTACTGGCATCGTTATCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTAAACATAACGTAACATGTAGAAGAACCTTCCTGACATTTCCAGGCTACTACAGCATAGTTGTCGtt       >  2:1354266/1‑150 (MQ=255)
                                                                                                                                                gggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGttt      <  2:615734/150‑1 (MQ=255)
                                                                                                                                                 ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTg     <  2:724333/150‑1 (MQ=255)
                                                                                                                                                 ggTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGAGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTg     <  2:1279781/150‑1 (MQ=255)
                                                                                                                                                  gTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGa    >  2:629935/1‑150 (MQ=255)
                                                                                                                                                    tACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGACt  <  2:1974027/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTTGCCGCGCGAACATAACGCCATGGGTTACTGGCATCGTTTTCTTTCGCATCGCCTGGGCAATTTGCCGCGCCTACCCACCAGGCGTTCTTATAAATGCCGCCGTTATAGATAACGTACCATGTAGAAGAACCTTCCTGACCTTTCCAGGCTACTACAGCATAGTTGTCGTTTGACT  >  NC_000913/3468892‑3469189

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: