Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,877,629 C→T 10.9% intergenic (‑97/‑86) yoaF ← / → dgcP DUF333 domain‑containing lipoprotein YoaF/diguanylate cyclase DgcP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,877,6290CT10.9% 192.8 / 2.4 72intergenic (‑97/‑86)yoaF/dgcPDUF333 domain‑containing lipoprotein YoaF/diguanylate cyclase DgcP
Reads supporting (aligned to +/- strand):  ref base C (38/24);  new base T (2/7);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 3.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGCTGACCCCTTGCGCTGTCTGCACTGGAATGAGCGTCCCGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCC  >  NC_000913/1877355‑1877907
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                      
tgttgCTGACCCCTTGCGCTGGCTGCACTGGAATTGGCATCCCGCCCTTCTGCTGGCAATAGACCGCCGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGttt                                                                                                                                                                                                                                                                                  <  2:101367/281‑1 (MQ=255)
              tGCGCTGTCTGCACTGGAATGAGCGTCCCGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTtgatg                                                                                                                                                                                                                                                                    >  1:338339/1‑281 (MQ=255)
                       tggagtGGAATGGTGGTCCCGCCGTTCTGCTGGCAAGAGACCGACGCCGGATTGACCATAACAATCGGCGGCGGTGTAGGTGCTTCTGGCTGAGAAGGAGTTGAACAACCAGCGAGGACCAGCAAGCAAGGCACAACAAAACTGATAAGTTTCATTTATTGATCTCACATATTTATCCAAGGTTTGAGTATCGCGGTATCGTTTTGTTTTGCGGAACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTGCAGCTTTGTTGATGAGGTTTTAt                                                                                                                                                                                                                                                           <  2:106182/276‑1 (MQ=255)
                               aaTGAGCGTCCCGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCACGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTAAAAACATATTACGTCTTGGTTCATCTTTGTTGATGATGTTTTATTATGCCTg                                                                                                                                                                                                                                                   >  2:195973/1‑281 (MQ=255)
                               aaTGAGCGTCCCGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTg                                                                                                                                                                                                                                                   >  2:257104/1‑281 (MQ=255)
                               aaTGAGCGTCCCGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTATTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTg                                                                                                                                                                                                                                                   >  1:333797/1‑281 (MQ=255)
                                  gaccGTCCCGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCaa                                                                                                                                                                                                                                                <  2:81190/278‑1 (MQ=255)
                                  gAGCGTCCCGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCaa                                                                                                                                                                                                                                                <  2:686439/281‑1 (MQ=255)
                                   aGCGTCCCGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGGATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTAGTTTCATCTTTGTTGTTGATGTTTTATCATGCCTGCaaa                                                                                                                                                                                                                                               >  1:196386/1‑281 (MQ=255)
                                        cccGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCCCTCAAAACATATTACGTCTTGTTTCATCTTTTTTGATGATGTTTTATCATGCCTGCAAAGATTa                                                                                                                                                                                                                                          >  1:421855/1‑281 (MQ=255)
                                        cccGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTTTTGTTTCATCTTTGTTTATGATGTTTTATCATTCCTGCAAAGATTa                                                                                                                                                                                                                                          >  2:970353/1‑281 (MQ=255)
                                              cTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAaataat                                                                                                                                                                                                                                    <  1:257104/281‑1 (MQ=255)
                                                           tAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACccg                                                                                                                                                                                                                       <  2:338339/281‑1 (MQ=255)
                                                             gACCCACGCCTGTTTTGCACTACCAATCTGGGGGGGTTAAGGGGTTTCTGTCTGGAAAGGGGTTAAACAACCAGCCAGGACCAGCAAGAAAGGCAAAACAAAACTAATAATTTTCATTTATTGTTCCTGCATATTTATCCAAGTTTGGAGTATCGCGGTATCTTTTTGTTTTGGAGCACTATTTTTATTACATTCCCTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATTTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACccgcc                                                                                                                                                                                                                     <  2:421855/281‑1 (MQ=255)
                                                                  aCGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCCAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATc                                                                                                                                                                                                                >  1:710966/1‑281 (MQ=255)
                                                                        gagtggCCATCCCAATCTGCGGGGGTTTAGGTGCTTCTGGCTGAGAAGAGTTTGAACAACCCACCAGGCCCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTCTCCAAGATTAGAGTATCGCGTTATCTTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAg                                                                                                                                                                                                          <  2:972630/276‑1 (MQ=255)
                                                                              ccATACCAATCTGCGGCGGTTTGGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAATGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCCCTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTc                                                                                                                                                                                                    <  2:503382/281‑1 (MQ=255)
                                                                              ccATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGGGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTc                                                                                                                                                                                                    <  1:1048957/281‑1 (MQ=255)
                                                                                  aCCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATc                                                                                                                                                                                                                                                                              <  1:804151/203‑1 (MQ=255)
                                                                                  aCCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATc                                                                                                                                                                                                                                                                              >  2:804151/1‑203 (MQ=255)
                                                                                       tCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGtcaga                                                                                                                                                                                           >  2:333240/1‑281 (MQ=255)
                                                                                                gTTTAGTGGCTTCGGGCTGAGAAGGGTTTGAAAAACCAGCCAGGACCAGCAAGCAAGGAAGAACAAAACTGATAATTTACATTTATTGATCTCACATTGTTATCCAAGATTAGATAATCGGGATATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATTTTTTTTGAGGATGTTTTATCATGCCTGCAAAGATTAAATAATAAGCATTTAACCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGAtt                                                                                                                                                                                    <  2:283783/279‑1 (MQ=255)
                                                                                                                   agaAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATcc                                                                                                                                                               >  2:131489/1‑281 (MQ=255)
                                                                                                                        gggTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATGTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGAACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCAGGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTg                                                                                                                                                           <  1:409509/280‑1 (MQ=255)
                                                                                                                         ggTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGTGTCTCGCATTCCCTTGcc                                                                                                                                                         >  1:1057285/1‑281 (MQ=255)
                                                                                                                                 aacCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGaa                                                                                                                                                   <  1:131489/279‑1 (MQ=255)
                                                                                                                                  acCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTTATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTACTACATTCACTCAAAACATAATACGTCTTGTTCCATCTTTGTTGATGATGTTTTATTATGCCTGCAAAGATTAAATTATCAGCATTTACCCGCCGTATCCAGGAGTTGTTCCGTTGCAGATCAGATTACCGCCCGCGTCTCGCAATCCCTTCCCAAAGAACAg                                                                                                                                                >  1:175138/1‑281 (MQ=255)
                                                                                                                                  acCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAg                                                                                                                                                >  1:362098/1‑281 (MQ=255)
                                                                                                                                                gcaagcaaGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCt                                                                                                                                  >  1:335335/1‑281 (MQ=255)
                                                                                                                                                        ggCAGAAAAAAACTGATAAATTTCATTTATTGAGCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTTTTGTTTCATCTTTGTTGATGAGGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGAc                                                                                                                          <  1:321275/281‑1 (MQ=255)
                                                                                                                                                                                 tttATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGt                                                                                                 >  1:823772/1‑281 (MQ=255)
                                                                                                                                                                                    attaatccccaatatTTATCCAAGATTAGAGTATCGCGGTATCGTTTTTTTTTGCAGCACTATTTTTATTACTTTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTCCCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCAc                                                                                              <  1:970353/270‑1 (MQ=255)
                                                                                                                                                                                      tgttttcACATATTTATCAAAGTTTAGAGTATCGCGTTATCGTTTTTTTTTCCAGCCCTATTTTTATTACATTCCCCCAAAACATATTACGTCTTTTTTCATCTTTGTTTATGATGTTTTATCATGCCGGCAAAGATTAAATAATCAGAATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTg                                                                                            <  2:823772/276‑1 (MQ=255)
                                                                                                                                                                                        ttttcAAATATTTTTCCAATATTAAGTTAGTGCGATTTCTTTTTTTTTTTCAGCATTATTTTTATTTCATTCACTCAAAACAAATTTCGTTTTGTTTAATCTTTGTTGATGATGGTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGata                                                                                         <  2:1057285/278‑1 (MQ=255)
                                                                                                                                                                                             taattttcacacAATATTATCGTATCACGTTTTGGTTTTGTTTTGAAGCACTTTTTTTATTACATTCCCTCAAAACATTTTTCGTTTTTTTTCTTCTTTGTTGATGATGTTTTATCATGCCTTCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGa                                                                                     <  2:362098/270‑1 (MQ=255)
                                                                                                                                                                                             acaTATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTTTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATAAGATTATCTCCCGCGTCCCGCAATCCCTTGCAAAAGAACAGCCACTGGAAAGTATGGTCCGACAGCTTTCGGAAATGCTGGAAATGGTACATGCTATGGa                                                                                     >  2:1006871/1‑281 (MQ=255)
                                                                                                                                                                                             acaTATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGa                                                                                     >  1:991261/1‑281 (MQ=255)
                                                                                                                                                                                              caTATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAAc                                                                                                                                                  <  2:885447/219‑1 (MQ=255)
                                                                                                                                                                                              caTATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAAc                                                                                                                                                  >  1:885447/1‑219 (MQ=255)
                                                                                                                                                                                                       ccAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAacctac                                                                           >  1:887191/1‑281 (MQ=255)
                                                                                                                                                                                                                    aTCGGGGTATCGTTTTGTCTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTTTTGTTTCATCTTTGTTGATAATGTTTTATCATGCCTGCAAAGATAAAATATTCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTgg                                                              <  1:505026/281‑1 (MQ=255)
                                                                                                                                                                                                                     tCGCGGTATCGTTTTGTTTTGCCGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGAACAGATTATCTCCCGCGTCTCGCAATCCCTTGCCAACGAACAGTCACTGGAAAGTCTGGTCCGACCGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGAACAACGTgg                                                              >  1:1038882/1‑280 (MQ=255)
                                                                                                                                                                                                                     tCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCCAAGTGGa                                                             >  1:997743/1‑281 (MQ=255)
                                                                                                                                                                                                                     tCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTgg                                                              >  1:171861/1‑280 (MQ=255)
                                                                                                                                                                                                                     tCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTCTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCCCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGCCAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGCATCAACCTACTTGAACAAAGTGGa                                                             >  2:662290/1‑281 (MQ=255)
                                                                                                                                                                                                                      cccgGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGACGATGTTTTATCATGCCTGCCAAGATTAAATAACCAGCATTTACCTGCCGTATTCTGGCGTTGCTCTGTGCTAGACCAGATTATCCACCGCGTCTCTCAATCTCTTGCCAAAGATCAGTCACTGGGAGGTCCGGCCCGAAAGCTTCTTGAGACGCTGGAAATGGGCACTGCTATAGAATCAAACTACACGACCAAAggggtg                                                            >  2:1032599/3‑279 (MQ=255)
                                                                                                                                                                                                                             tCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCa                                                                                               >  2:196097/1‑239 (MQ=255)
                                                                                                                                                                                                                             tCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCa                                                                                               <  1:196097/239‑1 (MQ=255)
                                                                                                                                                                                                                                ttttgttttgCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTACGACAGCTTCTGGAGAGTATGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGCGGAACTCGAAgcg                                                   >  1:408911/1‑280 (MQ=255)
                                                                                                                                                                                                                                  tttttttctatatcTATTTTTTTTACTTTCACTCAAAACAAATTCGGTTTTTTTTCATCTTTGTTGGTGGTGTTTTTTCATGCCTGCAAAGATTAACTAATCAGCATTTACCCGCCTTATCCGGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAAGCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAgcgcgc                                                <  2:991261/268‑1 (MQ=255)
                                                                                                                                                                                                                                       tttctcaccatttttttttTCATTCACTCTAAAATTCTTCCGTTTTTTTTTGTCTTTTTTTTTGTTTTTTTTCCTTGCCTGCAAAGATTAAATATTCATCATTTCCCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgca                                           <  2:196386/270‑1 (MQ=255)
                                                                                                                                                                                                                                       acatactgctttttttATTACTTTCCCTCAAAACATTTTTCGTCTTTTTTCATCTTTTTTGATGATGTTTTATCATGCCTGCAAAAATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCCGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTgca                                           <  1:903319/270‑1 (MQ=255)
                                                                                                                                                                                                                                                  tttttATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGAGCCTGCAGCATATTATGt                                >  1:360215/1‑281 (MQ=255)
                                                                                                                                                                                                                                                      tattaCATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGc                            <  1:333240/281‑1 (MQ=255)
                                                                                                                                                                                                                                                       attaCATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCTCGCCTGCAGCATATTATGTTTGc                            >  2:498032/1‑280 (MQ=255)
                                                                                                                                                                                                                                                       attaCATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGc                            >  1:160557/1‑280 (MQ=255)
                                                                                                                                                                                                                                                          aaaTTCACTCAAAACAGTTTACTCTTTGTTTCATCTTTGTTTATGAGGTTTTATCATTCCTACAAAGATTAAATAATACACATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGt                        <  2:887191/279‑1 (MQ=255)
                                                                                                                                                                                                                                                          aCATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATTCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTAACCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCACTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGAGGCTGGAAATGGTCACTGATAAGGAATCAACCTACCTGACAAAAGTGGAGGTCGAAGCGCCACTGCACCATATTAAGTCTgtctct                        >  2:903319/1‑276 (MQ=255)
                                                                                                                                                                                                                                                            aTTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCTTGTCAGATCAGATTATCCACCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGCTGCTGGAAATGGTCACTGATCTGGAATCAACCTACCTGACCAAAATGGACGCCGACGACCCCCCGCCGCATATCTAGTTTCCCCGTaa                      >  2:761233/1‑281 (MQ=255)
                                                                                                                                                                                                                                                            aTTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATg                                                                                       >  1:910535/1‑216 (MQ=255)
                                                                                                                                                                                                                                                              tCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCAt                                                                                                                                                                                                                                                        >  2:588029/1‑53 (MQ=255)
                                                                                                                                                                                                                                                              tCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCAt                                                                                                                                                                                                                                                        <  1:588029/53‑1 (MQ=255)
                                                                                                                                                                                                                                                              tCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGGCCCTGCAGCATATTATTTTTGCCCGTAAcc                    >  2:164634/1‑280 (MQ=255)
                                                                                                                                                                                                                                                                  tCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCaa                                                                                                                                                                                                                                                >  2:128699/1‑57 (MQ=255)
                                                                                                                                                                                                                                                                  tCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCaa                                                                                                                                                                                                                                                <  1:128699/57‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tatTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGAGGATGTAGAAGCGCGCCTGCAGCATATTAGGTTTGCCCGTAACAGCCAGAAAATGt        >  1:847348/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                          tatTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCACTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGAGCCTTCAGCATATTCTTTTTGCCCGGAACAGCCAGAAAATGt        >  2:145874/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                            ttACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAg                                                                                                                                                                                                                                              <  2:789956/49‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ttACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAg                                                                                                                                                                                                                                              >  1:789956/1‑49 (MQ=255)
                                                                                                                                                                                                                                                                            ttACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTAc      >  2:68995/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                             tACGTCTTTTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACa     <  1:68995/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                                gTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAACGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTCCCTGACCAAAGTGGAGGTCGAAGCGCGCCTGCAGCCTATTCGGTTTGCCCGTAACAGCCAGAAAATGTAAATcc  >  2:226077/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                      
TGTTGCTGACCCCTTGCGCTGTCTGCACTGGAATGAGCGTCCCGCCCTTCTGCTGGCAATAGACCGACGCCGGATTTGCCATACCAATCTGCGGCGGTTTAGGTGCTTCTGGCTGAGAAGGGGTTGAACAACCAGCCAGGACCAGCAAGCAAGGCAGAACAAAACTGATAATTTTCATTTATTGATCTCACATATTTATCCAAGATTAGAGTATCGCGGTATCGTTTTGTTTTGCAGCACTATTTTTATTACATTCACTCAAAACATATTACGTCTTGTTTCATCTTTGTTGATGATGTTTTATCATGCCTGCAAAGATTAAATAATCAGCATTTACCCGCCGTATCCTGGAGTTGTTCCGTGTCAGATCAGATTATCGCCCGCGTCTCGCAATCCCTTGCCAAAGAACAGTCACTGGAAAGTCTGGTCCGACAGCTTCTGGAGATGCTGGAAATGGTCACTGATATGGAATCAACCTACCTGACCAAAGTGGATGTCGAAGCGCGCCTGCAGCATATTATGTTTGCCCGTAACAGCCAGAAAATGTACATCC  >  NC_000913/1877355‑1877907

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