Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,776,132 A→C 5.6% S72R (AGC→CGC)  yfjX → CP4‑57 prophage; putative antirestriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,776,1320AC 5.6% 328.9 / 5.2 121S72R (AGC→CGC) yfjXCP4‑57 prophage; putative antirestriction protein
Reads supporting (aligned to +/- strand):  ref base A (59/53);  new base C (2/6);  total (61/60)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

CGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTC  >  NC_000913/2775983‑2776281
                                                                                                                                                     |                                                                                                                                                     
cGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACCCGCTCa                                                                                                                                                       >  2:962538/1‑150 (MQ=255)
          cGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCgg                                                                                                                                             <  2:241356/150‑1 (MQ=255)
            aTGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAg                                                                                                                                           <  2:1424027/150‑1 (MQ=255)
              gAACAGCATTTACGCTGCTGGCCGCCGCACTTTGGCACCAACCCGCAGTGGATAACGCTGGAACAGCGCATTTTCGCCGGGATGGACCGCTTCTGTGAGGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCCGCCATTGCTGAGCg                                                                                                                                         <  2:370138/150‑1 (MQ=25)
              gAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCGTTTTCGGCTGGATGGACCGGTTCTGTGATGAGTACTGCGGGGGTTTGTGGTGCTTTTACGCGCTCAGCAATGGCGGAGCg                                                                                                                                         >  2:192015/1‑150 (MQ=255)
              gAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCg                                                                                                                                         >  1:1619192/1‑150 (MQ=255)
               aaCAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGGATCTGGTCCTTTTTCACGCTCAGCAATGGCGGGGCgg                                                                                                                                        >  1:128558/1‑149 (MQ=255)
               aaCAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACAGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGt                                                                                                                                        <  1:1086702/150‑1 (MQ=255)
                aCAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATGACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCAGGTACAGTTACACGCTCAGCAATGGCGGAGCggg                                                                                                                                       >  2:399049/1‑148 (MQ=37)
                aCAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCCGCCATTGCTGAGCGtg                                                                                                                                       <  1:1765400/150‑2 (MQ=255)
                 cAGCATATCCGCTTCTGGCCGGAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCCGCCATTGCTGAGCGtgt                                                                                                                                      <  2:796926/150‑3 (MQ=37)
                 cAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGttt                                                                                                                                      <  1:2040386/150‑1 (MQ=255)
                 cAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGttt                                                                                                                                      <  1:1431553/150‑1 (MQ=255)
                  aGCCTATCCGCTTCTGGCCGCAGCTCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCCGCCATTGCTGAGCGTGTa                                                                                                                                     <  2:112480/150‑1 (MQ=37)
                   gCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTAt                                                                                                                                    >  1:1336514/1‑150 (MQ=255)
                   gCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTAt                                                                                                                                    >  2:647256/1‑150 (MQ=255)
                   ccATATCCGCTTCT‑GCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATg                                                                                                                                   <  1:1724742/149‑1 (MQ=255)
                         ccGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCAGTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTcc                                                                                                                              <  1:1550074/150‑1 (MQ=255)
                          cGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCt                                                                                                                             >  2:2039300/1‑150 (MQ=255)
                          cGCTTCAGGGCGCACCACTTTGGCAGCTTGCCGCAGTGGTTATCGCTGGATCCGGTCATTTTCGCCTGGTTGGAGCGCTTCTGTGATGAGTACTGCGGGGGTTTCTGGTGCTTTTACACGCTCTGCAATGGCGGTGCGTTTATGGGTCCt                                                                                                                             <  2:766621/150‑1 (MQ=255)
                           gCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTg                                                                                                                            >  2:1602297/1‑150 (MQ=255)
                           gCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCAATTTCGACTGGATGGACCGCTTCTGTGATGAGTAATGCGGGGGTAACTGGTCCTTTTACACGCTCAGCAATTGCGGAGCGTTTATGGctgctg                                                                                                                            >  2:1342699/1‑150 (MQ=37)
                                tGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGGGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCt                                                                                                                       >  1:1230692/1‑150 (MQ=255)
                                    cGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACg                                                                                                                   <  1:1759098/150‑1 (MQ=255)
                                    cGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCGTCTGTGATGAGTACTGCGGGGGGATCTGGGCCTGTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACg                                                                                                                   >  1:1405904/1‑150 (MQ=255)
                                         caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGAc                                                                                                              <  2:606929/150‑1 (MQ=255)
                                         caCTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGAc                                                                                                              <  1:250940/150‑1 (MQ=255)
                                         caCTATGGCACCCTCCCGCAGGGGATAACGCTGGAACCGCGCATTTTCGCCTGGAGGGCCCGCTTCTGTTATGAGTACTGCGGGGGTATCTGGTCCTTTAACACTCTCAGCAATGGCGGAGCGTTTATGGGTCCTGATTCTGACGGTGAc                                                                                                              <  2:112937/150‑1 (MQ=255)
                                             ttGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATa                                                                                                          <  2:1533874/150‑1 (MQ=255)
                                               ggCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATaaa                                                                                                        <  2:563040/150‑1 (MQ=255)
                                                 cACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATg                                                                                                      <  2:1840296/150‑1 (MQ=255)
                                                 cACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATg                                                                                                      <  2:1155319/150‑1 (MQ=255)
                                                  aCCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATgg                                                                                                     <  2:1506595/150‑1 (MQ=255)
                                                  aCCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGGGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATgg                                                                                                     >  2:412057/1‑150 (MQ=255)
                                                   ccATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACCCGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGc                                                                                                    >  2:66246/1‑150 (MQ=255)
                                                   ccATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTACTGTGATGAGTACTGCGCGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGc                                                                                                    >  2:547080/1‑150 (MQ=255)
                                                     aTCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCAt                                                                                                  >  1:1561811/1‑150 (MQ=255)
                                                      tCCCGCAGTGTATAACGCGGGAACCGCGCTTTTTCGCCTGTATGTTCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTAACACGCTCAGCAATGGCGTAGCGTTTATGGCTACTGATGCTGACGGTGACGATAAATGGCAtt                                                                                                 <  2:127691/150‑1 (MQ=255)
                                                      tCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGGTTATGGCTCCTGCTGCTGACGGTGACGATAAATGGCAtt                                                                                                 >  2:1302416/1‑150 (MQ=255)
                                                       cccGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGTTGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCAttt                                                                                                >  2:1440506/1‑150 (MQ=255)
                                                              tGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGCTGACGATAAATGGCATTTGCTCAAc                                                                                         >  2:486067/1‑150 (MQ=255)
                                                               ggATAACACTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACg                                                                                        <  1:389397/150‑1 (MQ=255)
                                                                taTAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCATTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACgg                                                                                       <  2:1598141/149‑1 (MQ=255)
                                                                  tAACGCTCGAACCGAGCATTTTCGCCTGGATGGACCGCTTCTGTGTTGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCATCGGCa                                                                                     <  2:224188/150‑1 (MQ=255)
                                                                        tGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATg                                                                               <  2:1734848/150‑1 (MQ=255)
                                                                         ggAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATgg                                                                              >  1:1395514/1‑150 (MQ=255)
                                                                         ggAACCGCGCATTTTCGCATGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATCGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATgg                                                                              >  1:1149347/1‑150 (MQ=255)
                                                                           aaCCGCGCATTTACGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGTGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCa                                                                            <  1:1052814/150‑1 (MQ=255)
                                                                               gcgcATTTTCGCCTGGATGGACCGCTTCTGTTATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCGCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATgg                                                                        <  2:516017/150‑1 (MQ=255)
                                                                                 gcATTTTCGCCTGGATTGAACGCTTCTGTGATGAGTACTGCGGTGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACTGTGACGATAAATGGCATTTTCTAAACGGCATGAATGGCAATGGTg                                                                      <  2:920943/150‑1 (MQ=255)
                                                                                    ttttCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCgg                                                                   >  2:925842/1‑150 (MQ=255)
                                                                                    ttttCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCgg                                                                   <  1:1076150/150‑1 (MQ=255)
                                                                                        cGCCTGGATGGACCGCTTCTGTGATGAGTACTGAGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAAt                                                               >  2:639829/1‑150 (MQ=255)
                                                                                         gCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATg                                                              <  2:2063141/150‑1 (MQ=255)
                                                                                         gCCTGGATGGACCGCTTCTGTGAGGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATg                                                              <  1:589736/150‑1 (MQ=255)
                                                                                          ccTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGa                                                             >  1:1296441/1‑150 (MQ=255)
                                                                                          ccTGGATGGACCGCTTCTGTGATGAGTACTGAGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGa                                                             <  2:42934/150‑1 (MQ=255)
                                                                                           gtggatggaCCGCTGCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAg                                                            <  2:1280773/149‑1 (MQ=255)
                                                                                             ggatggaCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAgcg                                                          >  2:6919/1‑150 (MQ=255)
                                                                                              gatggaCCGGTTCTGTGTTGAGTACGGCGTGGGTATCTTGTCCTTTTACACGCGCAGCAATGGCGGAGCGTTTATGTCTCCTGATGCTGACGGTGACGTTAAATGGCATTTGCTCAACGGCTTGAATGGCAATGGTGCGGAAATGAgcgc                                                         <  2:1518545/150‑1 (MQ=255)
                                                                                                 ggaCCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGa                                                      >  1:1665753/1‑150 (MQ=255)
                                                                                                    ccGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGc                                                   >  1:1617575/1‑150 (MQ=255)
                                                                                                          cTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAgg                                             <  1:157994/150‑1 (MQ=255)
                                                                                                          cTGTGAGGAGGACTCCGGGGGTATCTGGTCCTTTAACCCGCTCCGCAATGGCGGAGCGTTTATGGCTCCTGATCCTGACGGGGACGATAAATGGCATTTGCGCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAgg                                             <  1:147344/150‑1 (MQ=255)
                                                                                                            gtgATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTa                                           >  2:920367/1‑150 (MQ=255)
                                                                                                            gtgATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTa                                           >  1:1997068/1‑150 (MQ=255)
                                                                                                            gtgATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTa                                           >  1:548806/1‑150 (MQ=255)
                                                                                                              gatgaGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATc                                         >  1:1736256/1‑150 (MQ=255)
                                                                                                               atgaGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAATTGAGCGCGGAAGCCGCAGGTATCg                                        >  1:357284/1‑150 (MQ=255)
                                                                                                                tgaGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGTATCCTGAGGCTGACGGTGACGATAAAGGTAATTTGCTCAACGGAATGCATGGCAATGGTGAGGAAATTAGCTAGGAAGCCGCAGTTATCGc                                       >  1:1281382/1‑150 (MQ=255)
                                                                                                                tgaGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGc                                       >  1:1809976/1‑150 (MQ=255)
                                                                                                                  aGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTg                                     <  2:2101391/150‑1 (MQ=255)
                                                                                                                  aGTACTGCGGGGGTATCTGCTCCTTTTACACGCTCCGCAATGGCGGAGCGTTAATGGCTACTGATGCTGACGGTGACGATAAATGGCATATGCGCTCAGGCATGAATGGCAATGGTGCGGAAATTATCTCGGCAGACGCAGGTATCGCTg                                     >  2:347612/1‑150 (MQ=255)
                                                                                                                      cTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTg                                 <  2:18391/150‑1 (MQ=255)
                                                                                                                      cTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTg                                 <  2:1794702/150‑1 (MQ=255)
                                                                                                                      cTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTg                                 <  2:1693685/150‑1 (MQ=255)
                                                                                                                        gCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGcc                               >  2:429037/1‑150 (MQ=255)
                                                                                                                          gggggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCGCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTg                             <  2:1010324/150‑1 (MQ=255)
                                                                                                                             ggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAtt                          <  2:1913186/150‑1 (MQ=255)
                                                                                                                             ggTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAtt                          <  1:66246/150‑1 (MQ=255)
                                                                                                                              gTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAttt                         >  2:699967/1‑149 (MQ=255)
                                                                                                                                aTCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGaa                       >  1:1519219/1‑150 (MQ=255)
                                                                                                                                aTCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGaa                       >  1:951821/1‑150 (MQ=255)
                                                                                                                                 tCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAACTGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCATGTATCGCTGGCTGCCTGATTGAat                      <  2:128558/150‑1 (MQ=255)
                                                                                                                                 tCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAat                      <  1:1433306/150‑1 (MQ=255)
                                                                                                                                  cTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                     >  1:804247/1‑150 (MQ=255)
                                                                                                                                    ggTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                   >  1:261868/1‑150 (MQ=255)
                                                                                                                                      tCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAgt                 <  1:518326/150‑2 (MQ=255)
                                                                                                                                      tCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGc                 >  1:190885/1‑150 (MQ=255)
                                                                                                                                       ccTTTTACACGCTCAGCACTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGcc                <  2:1348424/150‑1 (MQ=255)
                                                                                                                                       ccTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGcc                <  2:1522646/150‑1 (MQ=255)
                                                                                                                                             acacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACg          <  2:619402/150‑1 (MQ=255)
                                                                                                                                             acacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTCACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACg          >  2:159699/1‑150 (MQ=255)
                                                                                                                                             acacGCTCAGCAATGGCGGAGCGTTTATATCTCCTGATGCTGCCTATCACTAAAAAATACATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGAATCGCTGTCTGCCTGATTGAATATAGCCATCCCg          >  2:1269095/1‑150 (MQ=255)
                                                                                                                                              cacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgc         >  1:1668875/1‑150 (MQ=255)
                                                                                                                                              cacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgc         >  1:328305/1‑150 (MQ=255)
                                                                                                                                                 gCTCAGCAATGGCGGAGCGTATATGGCTCCTGAAGCACACGGTGACGATAAATGGCAGTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGACGCCGAAGGTATCGCTGAGATCCTGATTGAATATAGCAATCACgcctg      >  1:833252/1‑150 (MQ=255)
                                                                                                                                                  cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGATTATAGCCATCACgcctgc     >  1:198369/1‑150 (MQ=255)
                                                                                                                                                  cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgc     >  1:891678/1‑150 (MQ=255)
                                                                                                                                                  cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgc     >  2:1829121/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCAGGTACCCGATGAACAGCATATCCGCTTCTGGCCGCAGCACTTTGGCACCATCCCGCAGTGGATAACGCTGGAACCGCGCATTTTCGCCTGGATGGACCGCTTCTGTGATGAGTACTGCGGGGGTATCTGGTCCTTTTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTC  >  NC_000913/2775983‑2776281

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: ‑