Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,159 T→A 6.8% intergenic (‑142/+29) yedQ ← / ← yodC putative membrane‑anchored diguanylate cyclase/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,1590TA 6.8% 295.9 / 10.7 88intergenic (‑142/+29)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (41/41);  new base A (3/3);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-02

CACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCA  >  NC_000913/2027888‑2028430
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                               
cACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                                                                                                                                              >  2:1378/1‑275 (MQ=255)
  cGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATcc                                                                                                                                                                                                                                                                            >  2:163436/1‑275 (MQ=255)
     cAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGc                                                                                                                                                                                                                                                                         <  4:53463/275‑1 (MQ=255)
       tgTGAGCAAAGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACTACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                                                                                                                                                                                                       <  1:17677/274‑1 (MQ=255)
             cAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAACCGGCCt                                                                                                                                                                                                                                                                       >  2:436724/1‑269 (MQ=255)
               aGGTGGAAAAAAGCAATACCACGATAAAGCAGACATTAACGACATCACCAGGCTTCAGGCGGCGTGCGACTTTTTTCAACCAGCTCTGGTTCTCCATTTTTGTCTCGTGCCGCACCCTGACGACTCCCCACCTGCTTTTTTATGCTTCTGGCCTAACCCCGCCCTATCCTCCCGCAAAACAGCGTAAATTTTTCAACGGTGCCATGATTAGACGCTACTATCACTCCTTCAAGGCATTAACACACTCAGCCAGAGGTACTCCGCCATGAAGTTaa                                                                                                                                                                                                                                                               >  4:98560/1‑275 (MQ=255)
                 gTGGAAAAAAGCAATACCACGTAAAACCAGAGATTAATGCCATGACAAGGCCCCAGGCGGCAGGCAAGTTTTTTCAACCAACTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAAAAGCGATAATTTTATAATGGGGGAATGATAAGAAGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                                                                                                                                                                                                   <  4:436723/269‑1 (MQ=255)
                          aGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGCTGTTAATGGTACCCtt                                                                                                                                                                                                                                                     >  2:419798/1‑274 (MQ=255)
                          aGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCtt                                                                                                                                                                                                                                                     >  2:218800/1‑274 (MQ=255)
                             aaTACCAGAATAAAGCAGAGATTAACGACATGCCCAGGCCCCAGGGGGCGTGGGAGTTTTTTCAACCATCTGTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                                                                                                                                                                                                 <  4:1378/275‑1 (MQ=255)
                              aTACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTTTCTCGTTCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCCTAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAGTGGTGGCATGATAAGACGCGAATATCATTCATTGCAGGCATTACCACACCAGGCCGGATGTGATCCGGCCTGCTGTTAATGATACCCTTTccca                                                                                                                                                                                                                                                >  4:430114/1‑275 (MQ=255)
                                            cAGAGATTAACGTCATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATTCCCTTTCACACTTCTTCCGCAGa                                                                                                                                                                                                                                   >  2:122241/1‑274 (MQ=255)
                                                     aCGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccc                                                                                                                                                                                                                         <  4:217241/275‑1 (MQ=255)
                                                                       aGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCCCCCTGACGACTCCCAACCTGCTTTTTTATGCTTCTGGCATAACCTGCCGATAGCGTACCGCAACACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGCCTATCATTCATTGCAGGCATTAGCACACCAGTCCGGATGTTATCCGCCCTGATGTTAATGATACCCTTTCCACCTTCTTCCCCCGCACGACCCCCCTCCCCCggcaacagcca                                                                                                                                                                                                        >  4:374482/1‑265 (MQ=255)
                                                                       aGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCACCTCCCCCGGAACAAGCt                                                                                                                                                                                                         >  2:47896/1‑273 (MQ=255)
                                                                        ggcggcGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                                                                                                                                                                                                        >  4:505290/1‑209 (MQ=255)
                                                                        ggcggcGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACACGCTCgg                                                                                                                                                                                                      >  2:392723/1‑275 (MQ=255)
                                                                            gcgTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                                                                                                                                                                                                    <  2:505291/209‑1 (MQ=255)
                                                                                   aGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCCTTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACCATTCTTCCGCAGAAAGACTCCCCTCCCCCGGAACAAGCCAGGTTTCCTTaaa                                                                                                                                                                                            >  4:409738/1‑274 (MQ=255)
                                                                                       ttttCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTTATCCGGCCTGATGTTAATGATAATTCTCCACACTTCTTCCGTAAAACGACTCTCCTCCCCCGGAACAATCTCGGTTTAATGAAAAGCtt                                                                                                                                                                                       >  2:384270/1‑275 (MQ=255)
                                                                                                                ttttGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCAt                                                                                                                                                               <  4:218799/274‑1 (MQ=255)
                                                                                                                 tttGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATa                                                                                                                                                              >  4:147098/1‑274 (MQ=255)
                                                                                                                    gTCTCGTGATGCACCCCGCCGACTCCCAACCTGCTTTGTTATGATTCGTGCATAACCTGTCGATAGCGTACCGCAAAAAAGCGATAATTTTATAATGGTGGCATGATAAGACGCTACTATCATTCATTGCAGGCATTAGCACCCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCaa                                                                                                                                                          <  1:2045/275‑1 (MQ=255)
                                                                                                                     tctcGTGCTGCAGCCTGACGCATCCCAACCTGCTTTTTGATGTTTCTGGCATAACGTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAAc                                                                                                                                                         <  4:163435/275‑1 (MQ=255)
                                                                                                                         ttgcgtcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAACGTACCGCAAAACAGCGATAATTTTATAATGTTGGCATAGTAAGACGCGACTATCATTCATTGCAGTCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAATCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGaca                                                                                                                                                     <  4:122240/269‑1 (MQ=255)
                                                                                                                          tgctgcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGTATGTGATACGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAANGACTCCCCTCCCCCGGAAAAAGCTAGGTTTACTTAAAAGCTCCCCCCCTTACCCCATAACCGACAtcacaaaccacc                                                                                                                                                    >  3:132915/1‑264 (MQ=255)
                                                                                                                           gctgcACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACt                                                                                                                                                   <  1:91882/275‑1 (MQ=255)
                                                                                                                                 cccTGACGGCTCCCAACTTGCTTTTTTATGTTTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCa                                                                                                                                             <  4:199911/275‑1 (MQ=255)
                                                                                                                                 cccTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTAACGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAGCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCa                                                                                                                                             <  2:421724/275‑1 (MQ=255)
                                                                                                                                    tGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCt                                                                                                                                                                                                                                                                       >  1:138393/1‑150 (MQ=255)
                                                                                                                                    tGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACACTTTCACACTTCTTCCGCATAACGACTCCCCTCCCCCGGAACAAGCTCGGTTCCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATAc                                                                                                                                          >  2:306085/1‑275 (MQ=255)
                                                                                                                                     gacgacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCNGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGACCGACTCCCCTCCNCCGTAACAAGCCCGTTTTCATGCAAAGCTTCCCCCTTTACCCCGTNANCCTCACACACCCAACACCCAACCAAAcc                                                                                                                                         >  3:40659/1‑275 (MQ=255)
                                                                                                                                      ccgacGCCCAACCTTATTTTTTAGAATGTNNNCATAACCTGGGGAATGCGGACCGAAAACCAGCGATAATTTTATAAGGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACCACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCg                                                                                                                                        <  3:146473/274‑1 (MQ=255)
                                                                                                                                      acgacTCCCAACCGGCTTTTTTATGATTCTGGCATAGCCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCg                                                                                                                                        <  1:62497/275‑1 (MQ=255)
                                                                                                                                        gacTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                                                                                                                                                                                                   <  3:138393/150‑1 (MQ=255)
                                                                                                                                         acTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGaca                                                                                                                                                     <  1:40659/259‑1 (MQ=255)
                                                                                                                                           tCCCAACATGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTAACGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGa                                                                                                                                   <  1:132915/275‑1 (MQ=255)
                                                                                                                                            cccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATTTGATCCgg                                                                                                                                                                                                                                                                          <  5:967470/139‑1 (MQ=255)
                                                                                                                                            cccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCgg                                                                                                                                                                                                                                                                          >  6:967470/1‑139 (MQ=255)
                                                                                                                                            cccAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACAACACCCAACCATACCCCTGGAg                                                                                                                                  >  4:340838/1‑275 (MQ=255)
                                                                                                                                              cAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                                                                                                                                              >  8:656765/1‑133 (MQ=255)
                                                                                                                                              cAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                                                                                                                                              <  7:656765/133‑1 (MQ=255)
                                                                                                                                                        ttttATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                                                                                                                                                                                              <  2:394242/123‑1 (MQ=255)
                                                                                                                                                               aTTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCCTAACCGTCATACCAACGACACTCAACCATACCGCCGAACTCACCAGTGACAatcat                                                                                                                 >  4:283056/1‑273 (MQ=255)
                                                                                                                                                                       aTAACCTGGGTATGGCGTACCGCAAAACAGCGAGAATTGTATAATGGTGGCATGATAAGACGCGACTAGCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGgccg                                                                                                       <  4:392722/275‑1 (MQ=255)
                                                                                                                                                                         aaCCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGgccgcc                                                                                                     <  2:409739/275‑1 (MQ=255)
                                                                                                                                                                         aaCCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGgccgcc                                                                                                     <  2:147099/275‑1 (MQ=255)
                                                                                                                                                                                     gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGCCCCCCCTCTTTACCCg                                                                                          >  4:110093/1‑274 (MQ=255)
                                                                                                                                                                                      cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCTCCCTCTTTAACCGt                                                                                         >  2:40940/1‑274 (MQ=255)
                                                                                                                                                                                                aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                                                                                                                                                                                             >  5:288773/1‑148 (MQ=255)
                                                                                                                                                                                                 aCAGCGATAATTTTATAAGTGTGGCATGATAAGACGCGACTATCATTCATGGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTGCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACt                                                                             <  1:65353/275‑1 (MQ=255)
                                                                                                                                                                                                        taatatttaaaTGGTGGCATGATAAGACGCGACTATCATTCCTTGCGGGCATTAGTCCACCAGGCCGGATGTGATCCGGCCTGAGGTTAATGATACCCTGTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGATCAAGCTCGGTTTCATGAAAAGCTACCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATGCCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATa                                                                      <  4:47895/266‑1 (MQ=255)
                                                                                                                                                                                                            tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGcc                                                                                                                                                                                                                                                                        <  7:570885/77‑1 (MQ=39)
                                                                                                                                                                                                            tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGcc                                                                                                                                                                                                                                                                        >  8:570885/1‑77 (MQ=38)
                                                                                                                                                                                                              tataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTATTGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTAATACCAACGACACTCAACCATACCGCTGGATTATCCAGTGACAATCATCCGCGGGCCGCCCCCTTTAACCGTACCTTCCTCCCAAACCATAAAGCTc                                                                >  4:45312/1‑275 (MQ=255)
                                                                                                                                                                                                                           tGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                                                                                                                                                                                             <  6:989661/73‑1 (MQ=255)
                                                                                                                                                                                                                           tGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                                                                                                                                                                                             >  5:989661/1‑73 (MQ=255)
                                                                                                                                                                                                                             aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                                                                                                                                                                                                   >  7:66117/1‑81 (MQ=255)
                                                                                                                                                                                                                             aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTc                                                                                                                                                                                                                                                   <  8:66117/81‑1 (MQ=255)
                                                                                                                                                                                                                                      gACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTa                                                                                                                                                                       <  6:288773/148‑1 (MQ=255)
                                                                                                                                                                                                                                       aCTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGAcc                                                                                                                                                                            >  2:553628/1‑142 (MQ=255)
                                                                                                                                                                                                                                        cTAACATGCAGTTCAGGCATTAGCACACCATGCCGGATGTGATCCGGCCTGATGTTAATGCTACCCTGTCACACGTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGaa                                      <  2:361/275‑1 (MQ=255)
                                                                                                                                                                                                                                           tcattcatTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGt                                                                                                                                                                        <  4:553627/142‑1 (MQ=255)
                                                                                                                                                                                                                                              ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                                                                                                                                                                                                           <  5:663062/56‑1 (MQ=255)
                                                                                                                                                                                                                                              ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                                                                                                                                                                                                           >  6:663062/1‑56 (MQ=255)
                                                                                                                                                                                                                                              ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                                                                                                                                                                                                           >  8:433894/1‑56 (MQ=255)
                                                                                                                                                                                                                                              ttcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGAt                                                                                                                                                                                                                                                           <  7:433894/56‑1 (MQ=255)
                                                                                                                                                                                                                                               tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTCAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                                                                                                                                                                                             <  8:427432/101‑1 (MQ=255)
                                                                                                                                                                                                                                               tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTCAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                                                                                                                                                                                             >  7:427432/1‑101 (MQ=255)
                                                                                                                                                                                                                                                 attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt                                                                                                                                                                                 >  5:457827/1‑127 (MQ=255)
                                                                                                                                                                                                                                                 attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCt                                                                                                                                                                                 <  6:457827/127‑1 (MQ=255)
                                                                                                                                                                                                                                                 attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGtt                               <  2:316050/273‑1 (MQ=255)
                                                                                                                                                                                                                                                     cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTa                                                                                                                                                                                                                                                                >  2:532878/1‑44 (MQ=255)
                                                                                                                                                                                                                                                         cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGa                                                                                                                                                                                                                                                            <  4:532877/44‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                                                                                                                                                                                                   <  5:65486/35‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                                                                                                                                                                                                   >  6:65486/1‑35 (MQ=255)
                                                                                                                                                                                                                                                                  acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTaa                                                                           >  4:79979/1‑212 (MQ=255)
                                                                                                                                                                                                                                                                   tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                                                                                                                                                                                                  <  8:10033/91‑1 (MQ=38)
                                                                                                                                                                                                                                                                   tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                                                                                                                                                                                                  >  7:10033/2‑92 (MQ=38)
                                                                                                                                                                                                                                                                   tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                                                                                                                                             >  7:463616/2‑129 (MQ=255)
                                                                                                                                                                                                                                                                   tcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATAc                                                                                                                                                             <  8:463616/128‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGaa                                                                                                                                                                                                                                  <  7:138048/60‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ccAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGaa                                                                                                                                                                                                                                  >  8:138048/1‑60 (MQ=255)
                                                                                                                                                                                                                                                                      ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                                                                                                                                                                                            >  7:254062/1‑95 (MQ=255)
                                                                                                                                                                                                                                                                      ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaa                                                                                                                                                                                            <  8:254062/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ggCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCAt                                                                       <  2:79980/212‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGc   >  4:151075/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                            aTGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTACCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCTCGCTTGACCCCGTAACCGTCATACCAACGACACCAAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAAACTTAATTTCCCCACTAACCATAAAGCCCATACTCGCACACTTTTTTCGAGTGAAACTTGCTCACCTTTGTTAAAGTTTGCGAATTTTgga  >  4:532743/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                               
CACGCCAGGTGAGCAAGGTGGAAAAAAGCAATACCACGATAAAGCAGAGATTAACGACATGACCAGGCCCCAGGCGGCGTGCGAGTTTTTTCAACCAGCTCTGGTTTTCCATTTTTGTCTCGTGCTGCACCCTGACGACTCCCAACCTGCTTTTTTATGATTCTGGCATAACCTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACTTGTTCACCTTAGTTGAAGATGGCGAATTTTGCA  >  NC_000913/2027888‑2028430

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