Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,118,831 C→A 5.4% intergenic (+71/+14) zapB → / ← rraA FtsZ stabilizer, septal ring assembly factor, cell division stimulator/ribonuclease E (RNase E) inhibitor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,118,8310CA 5.4% 324.2 / 4.7 117intergenic (+71/+14)zapB/rraAFtsZ stabilizer, septal ring assembly factor, cell division stimulator/ribonuclease E (RNase E) inhibitor protein
Reads supporting (aligned to +/- strand):  ref base C (42/68);  new base A (4/3);  total (46/71)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

ACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCA  >  NC_000913/4118684‑4118977
                                                                                                                                                   |                                                                                                                                                  
acaacCATCTGAAAGCACAGCAGAACGGCTGTCAGGACCGGCTTCAGTCCCTGCTGGGTCGCATGGAAGATGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGGTGAAAAAAGGCACCTTGCg                                                                                                                                                   <  2:274189/149‑1 (MQ=255)
    ccATCTGAAATAACAGCAGAACTGCTGGCAGGAACGTCTGCAGGCCCTGCTTGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGc                                                                                                                                               <  1:496134/149‑1 (MQ=255)
       tCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCtt                                                                                                                                            >  2:558860/1‑149 (MQ=255)
        cTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGCGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCttt                                                                                                                                           <  1:796176/149‑1 (MQ=255)
         tGAAAGAAGAGCAGAACGGCTTGCAGGATCGACTGCAGGCCCTGCTGGGTCGCACGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTCCTGTGTAGGCCAAGAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCtttt                                                                                                                                          <  2:868344/149‑2 (MQ=255)
         tGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTc                                                                                                                                          <  2:210598/149‑1 (MQ=255)
          gAAAGTACTTCAGAACTTCTGGCAGGAACGTCTTAAGGCCATGATGGGTCTCTTGGAAGAGGTCTGATTACGCTGCATCTTCAGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCt                                                                                                                                         <  2:1420441/149‑1 (MQ=255)
            aacaaCAGCAGAACGGCTGGCACGAACGTCTGCAGGCCCTGATTGGTAGCATGTAATATTCCAGATTACGGTTCATCTTCAGCGTATTTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGTCACCTTGCGTTGCCTTTCTTa                                                                                                                                       <  2:1400093/146‑1 (MQ=255)
            aaGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTa                                                                                                                                       <  1:143431/149‑1 (MQ=255)
             aGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCGCAAGGTGCCttttttct                                                                                                                                      <  1:670370/149‑3 (MQ=255)
                 cagcagAACGGCTGGCAGGAACGTCTTCAGGCCCTGCTGGGGCGCATGGAAGAGGTCTGATTACGCGTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAtt                                                                                                                                  <  2:1227342/149‑1 (MQ=255)
                   gcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTGTTTAATGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTGTTATCATTCa                                                                                                                                >  2:166880/1‑149 (MQ=255)
                   gcagAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCa                                                                                                                                >  2:1130620/1‑149 (MQ=255)
                   gcagAACGGCGGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGTCACCTTGCGGTGCCTTTCTTATCATTCa                                                                                                                                <  1:486221/149‑1 (MQ=255)
                       aaCGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCAACTTGCGGTGCCTTACTTATCATTCAataa                                                                                                                            >  1:696562/1‑148 (MQ=255)
                          ggCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCa                                                                                                                         >  2:834510/1‑149 (MQ=255)
                             tGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCg                                                                                                                      <  2:445882/149‑1 (MQ=255)
                                cAGGAACGTCTGCAGTCCCTGCTGGGTCGAATGGTAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAt                                                                                                                   <  1:1376184/149‑1 (MQ=255)
                                cAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAt                                                                                                                   <  1:464000/149‑1 (MQ=255)
                                   gAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtctt                                                                                                                <  1:252625/149‑1 (MQ=255)
                                   gAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCAAGCGGAtctt                                                                                                                >  1:405424/1‑149 (MQ=255)
                                    aaCGTATGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGAAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtcttc                                                                                                               <  1:1078062/149‑1 (MQ=255)
                                     aCGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAtcttct                                                                                                              <  2:888438/149‑1 (MQ=255)
                                      cGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTg                                                                                                             <  2:250824/149‑1 (MQ=255)
                                       gTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGa                                                                                                            >  1:14949/1‑149 (MQ=255)
                                        tCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCGAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGAActgctgta                                                                                                           >  1:1444572/1‑147 (MQ=255)
                                        tCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGaa                                                                                                           >  1:1124640/1‑149 (MQ=255)
                                        gcggCAGGCCCTGCGGGGTCGCATGAAGACGTCTGATTACGCGTCAGCATCCGCGTAATTTCTTTCCTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGaaa                                                                                                          <  2:1433701/146‑1 (MQ=255)
                                         cTGCAGGCCCTGCTGGGTCGCATGGAAGCGGGCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGaaa                                                                                                          <  2:1061185/149‑1 (MQ=255)
                                         cTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGaaa                                                                                                          <  1:1011346/149‑1 (MQ=255)
                                          cgCAGCCCCTGCTGGGTCGCATGGAAGAGCTCTGCTTCCGCGTCATCTTCCGCGGAATTGCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAACGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAg                                                                                                         <  2:327261/148‑1 (MQ=255)
                                           gCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTGTCATTCAATATCCAGCGGATCTTCTGAAAGa                                                                                                        <  1:1294075/149‑1 (MQ=255)
                                                cccTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGaataat                                                                                                   >  1:473125/1‑149 (MQ=255)
                                                 ccTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATc                                                                                                  >  1:245507/1‑149 (MQ=255)
                                                  ctgctgGGTCGCATGGAAGATGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                                 <  2:595504/149‑1 (MQ=255)
                                                  ctgctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATcc                                                                                                 >  1:616338/1‑149 (MQ=255)
                                                  ctgctgGGTAGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATACAGCGGATCTTCTGAAATAATAATcc                                                                                                 >  2:541817/1‑149 (MQ=255)
                                                     ctgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCgg                                                                                              <  2:816768/149‑1 (MQ=255)
                                                      tgGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGt                                                                                             <  1:1566657/149‑1 (MQ=255)
                                                          tCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTg                                                                                         >  2:925198/1‑149 (MQ=255)
                                                           cGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGt                                                                                        <  2:1509742/149‑1 (MQ=255)
                                                            gCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTc                                                                                       <  2:750016/149‑1 (MQ=255)
                                                             cATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAATCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCg                                                                                      >  2:1178852/1‑149 (MQ=255)
                                                             cATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAAATTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCg                                                                                      >  2:664128/1‑149 (MQ=255)
                                                              aTGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCgg                                                                                     <  1:568285/149‑1 (MQ=255)
                                                                ggAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCa                                                                                   <  1:1360844/149‑1 (MQ=255)
                                                                 gAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCAt                                                                                  >  2:1342439/1‑149 (MQ=255)
                                                                   agagGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATaa                                                                                <  1:1028912/149‑1 (MQ=255)
                                                                     agGTCTGATTACGCTTCATCTGCCGCGTAATTTCATTACTGTGTAGGCCAGCAGGTAACGCAGAAAAAAGGCACCTAGAGGTGCCTTTCGTATCATTCAATATCCACCGGATCTTCTGAAAGAATAATCCCGGTATTCTCGGCATAAAg                                                                              >  1:342696/1‑149 (MQ=255)
                                                                      ggTCTGATTACGCCTCATCTTCCGCGTAATTTCTTAACTGTGTAGGCCAACAGGTAACGCAGAAAAAACGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGa                                                                             <  2:289102/149‑1 (MQ=255)
                                                                        tCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATg                                                                           <  2:1442200/149‑1 (MQ=255)
                                                                        tCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATg                                                                           <  2:445368/149‑1 (MQ=255)
                                                                         cTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCCGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATgg                                                                          <  1:1177626/149‑1 (MQ=255)
                                                                         cTGATTACGATTCATCTTCCGCGTAATTTCTTTACTGTGTATGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATgg                                                                          <  1:840315/149‑1 (MQ=255)
                                                                          tGATTACGCTTCATCTTCCGCGTAATTTCTATACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGt                                                                         <  1:103906/149‑1 (MQ=255)
                                                                           gATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCAAGCGGATCTTCTGAAAGATTAATCCCGGTATTGTCGGCATAAAGATGGTc                                                                        >  2:603819/1‑149 (MQ=255)
                                                                            aTTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTGCTGAAAGAATAATCCCGGTATTGTCGGGATAAAGATGGTCg                                                                       >  2:1483155/1‑149 (MQ=255)
                                                                               aCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCg                                                                    <  1:570121/149‑1 (MQ=255)
                                                                               aCGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCg                                                                    <  1:1349515/149‑1 (MQ=255)
                                                                                 gCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGAGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAAAAATCCCGGTATTGTCGGCAAAAAGATGGGCGCAGga                                                                  >  1:855726/1‑149 (MQ=255)
                                                                                  cTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCGTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGgag                                                                 >  2:1483574/1‑149 (MQ=255)
                                                                                  cTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGAAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGTTGGTCGCCGgag                                                                 <  1:400105/149‑1 (MQ=255)
                                                                                    tCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGagaa                                                               <  1:834510/149‑1 (MQ=255)
                                                                                      aTCTTAGGCGTAATTTATTTACTGTGTAGGCCAACAGGTAACGCAGAAAAACTGAACCTTGCGGTGCCTTTCTTATCATTCTATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGagaaga                                                             <  1:919515/149‑1 (MQ=255)
                                                                                        cTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGaagaag                                                           <  2:471047/149‑1 (MQ=255)
                                                                                          tCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGGAGAAAAAAGGCACCTTGCGGTTCCTTTCTTATCATTCAATATCCATCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGt                                                         <  2:1436202/149‑1 (MQ=255)
                                                                                          tCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGt                                                         <  2:598595/149‑1 (MQ=255)
                                                                                            cgcgTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGa                                                       <  2:1149101/149‑1 (MQ=255)
                                                                                             gcgTAATTTCTTTACTGTGTAGGCCAACAGGAAAAGGAGAAAAAAGGCACCTTGCGGTGCCTATCTTATCATTCAAAATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGac                                                      >  2:76284/1‑149 (MQ=255)
                                                                                                tAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGGCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATAACGGTATTGAAGTCATAAAGATGGTCGCCGGAGAAGAAGGTGACAcc                                                   >  1:1178272/1‑149 (MQ=255)
                                                                                                tAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACAcc                                                   >  1:1195300/1‑149 (MQ=255)
                                                                                                  aTTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGc                                                 <  2:701408/149‑1 (MQ=255)
                                                                                                      cTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCaaa                                             >  2:407043/1‑149 (MQ=255)
                                                                                                       tttACTGGGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCATTATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGGGACACCGCCaaaa                                            <  1:1344658/149‑1 (MQ=255)
                                                                                                        ttACTGTGTAGGCCAACAGGTAATGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAAt                                           >  2:12694/1‑149 (MQ=255)
                                                                                                        ttACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAAt                                           >  2:582991/1‑149 (MQ=255)
                                                                                                          agtgtgTAGGCCTACAGGTATCTAATAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGCATCTTCTGAAAGAATAATCCCTGAATTGGCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTg                                         <  2:1298163/147‑1 (MQ=255)
                                                                                                          aCTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTg                                         <  1:1483574/149‑1 (MQ=255)
                                                                                                           cTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGa                                        <  1:926045/149‑1 (MQ=255)
                                                                                                           cTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGa                                        <  1:57994/149‑1 (MQ=255)
                                                                                                               gtAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgcg                                    <  1:53847/149‑1 (MQ=255)
                                                                                                                   gCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCATCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAATATGGTCGCCGGAGAAGAACGTGACACCGCCAAAATTGACGCGGACa                                >  2:931350/1‑149 (MQ=255)
                                                                                                                     cAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCAATCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATc                              <  2:412285/149‑1 (MQ=255)
                                                                                                                      aaCAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCg                             <  2:320868/149‑1 (MQ=255)
                                                                                                                      aaCAGGTAACGCAGAAAAAAGGCACCTAGCGGTGCCTTTCATATCATTCAATATCCAGCGGATCTTCTGAAAGCATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGCAGCAGGTGACCCCGCCAAAATTGACGCGGACATCg                             >  2:491591/1‑149 (MQ=255)
                                                                                                                      aaCAGGTAACGCAGAAAAAACGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTAATGTCGGAATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCAGACATCg                             >  1:504398/1‑149 (MQ=255)
                                                                                                                           gTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTc                        >  2:369909/1‑149 (MQ=255)
                                                                                                                            tAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCg                       <  1:1114580/149‑1 (MQ=255)
                                                                                                                                gCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGCCACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                   >  1:778388/1‑149 (MQ=255)
                                                                                                                                gCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                   >  2:1002160/1‑149 (MQ=255)
                                                                                                                                 cAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAAt                  >  1:218897/1‑149 (MQ=255)
                                                                                                                                  aGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATg                 >  1:1238230/1‑149 (MQ=255)
                                                                                                                                    aagaaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc               >  1:528100/4‑149 (MQ=255)
                                                                                                                                    aagaaaGGCACCGCAAGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc               >  2:1277383/4‑149 (MQ=255)
                                                                                                                                    aaaaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc               >  1:1363322/1‑149 (MQ=255)
                                                                                                                                     aaaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCt              >  1:27686/1‑149 (MQ=255)
                                                                                                                                         aGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGc          <  2:1383820/149‑1 (MQ=255)
                                                                                                                                              ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCTGGATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTCACACCGCCAAAATTGACGCGGACATCGCGATCGCCAAAGCCTGCGCCAGCg     >  2:880584/1‑149 (MQ=255)
                                                                                                                                              ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCg     <  2:1529845/149‑1 (MQ=255)
                                                                                                                                                 tGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCa  <  1:1501525/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACAACCATCTGAAAGAACAGCAGAACGGCTGGCAGGAACGTCTGCAGGCCCTGCTGGGTCGCATGGAAGAGGTCTGATTACGCTTCATCTTCCGCGTAATTTCTTTACTGTGTAGGCCAACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCA  >  NC_000913/4118684‑4118977

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: