Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,127,577 C→T T106T (ACG→ACA yiiX ← putative lipid binding hydrolase, DUF830 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,127,5770CT98.2% 170.9 / NA 55T106T (ACG→ACAyiiXputative lipid binding hydrolase, DUF830 family protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/30);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCCGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGC  >  NC_000913/4127284‑4127865
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                
tATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCCCGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGACCAGTACTGGCGATCGTCGCCCCATGAGAAGCTGAAATCATACGGTTTACCGAGTTAACGTTTTGCTGTTTgcg                                                                                                                                                                                                                                                                                           >  2:398129/1‑301 (MQ=255)
   tACCAGCAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGGTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGcc                                                                                                                                                                                                                                                                                         <  2:566291/300‑1 (MQ=255)
        cGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTg                                                                                                                                                                                                                                                                                             <  1:264848/291‑1 (MQ=255)
          aaaaCAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTTGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGCACAGTCCTGGCGATCGTCGCCCCATGCAACGCCGAAATCCTCCGGTTTTCCGACAAAACGTTTTACTGCTTGCGCCGGGTTTTg                                                                                                                                                                                                                                                                                 >  1:559105/1‑301 (MQ=255)
           aaaCCGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAACCGGCTGGTGGCTACACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGTGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGc                                                                                                                                                                                                                                                                                <  1:208019/301‑1 (MQ=255)
                gCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCCTTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTCGCGCTCGTCGCCCCATGGCCATCTGAAATCATCCGGGTTACCGAGATAACGTTTTCCTGTTTGCGCCAGTTTTTGCtgttg                                                                                                                                                                                                                                                                           >  2:556658/1‑301 (MQ=255)
                 ccgCCATTCTTTGGCTACCGTGGTTAGTTGTCGCGCATCTACAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCtgttgt                                                                                                                                                                                                                                                                          <  1:13116/299‑1 (MQ=255)
                  gtcctttctttgcgagcagtggtgGGTTGGGGCGCATCGAAAACGGCTTTGGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTgttgtt                                                                                                                                                                                                                                                                         <  2:373016/283‑1 (MQ=255)
                             ggCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCCCTTCCGAACAGTCCCGGCGATCGTCTCCCCATGCGAAGCTGAAATCCTACGGTTTACCGAGATAACCTTTTGCTGTTTGCGCCAGTTTTTGCCGTTGTTCAACACTCAGt                                                                                                                                                                                                                                                              >  2:540927/1‑301 (MQ=255)
                             ggCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGCTAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGt                                                                                                                                                                                                                                                              >  1:87154/1‑301 (MQ=255)
                                gACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCCTACGGTTTACCGTGATAACGTTTTGCTGTTTGCTCCGGTTTTTGCTGTTGTTCAACACTCAGTACg                                                                                                                                                                                                                                                           >  2:238747/1‑301 (MQ=255)
                                gACCGTGGTGAGTTGGGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTc                                                                                                                                                                                                                                                                 <  1:238811/295‑1 (MQ=255)
                                  ccGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTCCTGGTAAACCTTCCACCCCACTTCCGCACAGTACTGGCGACCGTCGCTCCATGAGAACCTGAAATCATCCGGTTTACCCAGCACACATTTTCCTGTTTGTTCCAGGTTTTGCTTTTGTTTCACCCTCaacccccc                                                                                                                                                                                                                                                         >  2:337322/1‑294 (MQ=255)
                                          aattGGGGCCCCTCGAAACCTGCTGGTTGCTAGCCGCCCGTTTCCGCAAGCGGGAAGTTGTTGCCAGAACGTCCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGACCTCTTTCAACTTCTGCTGCTCGCCCACTGGCATTCCCAGCGCGTTCTGATAAGCCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAAcgc                                                                                                                                                                                                                                                 <  2:229299/299‑1 (MQ=255)
                                                   gcATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACCTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGCTCGTCGCTCCATGAGAAGCTGACACCATACGGTTTACCGAGATTACGTTTTGCCGTTTGCGCCCGTCTTTGCTGTTGGTCAACACCCAGTCCGCCTTCTCCGCGGCGaccc                                                                                                                                                                                                                                         >  2:227600/1‑297 (MQ=255)
                                                   gcATCGAAAACGGCTGGTAGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGaacaa                                                                                                                                                                                                                                        <  1:36868/301‑1 (MQ=255)
                                                    cATCGAAAACGGCTTGTGGCGAGACCGCCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGaacaac                                                                                                                                                                                                                                       <  1:398222/301‑1 (MQ=255)
                                                          gaaaggtttGTGGCGTAACCACCGTTGCCTAGAGCGGAATATTCTTGCCGTAACGTTCTTTGCGCTATGCCGGAACTTGTGGGCTGCGCAGGTCGAACTCTTTCAACTTCTGCGGCTCGCCCACGCGCAGTCCCAGCGCGTTCTGGTCAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTGGTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTTAACGCGGCGAACAACGTAttt                                                                                                                                                                                                                                 <  2:558973/294‑1 (MQ=255)
                                                                cTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGTTTTACAGAGATAAAGTTTTGATGTTTGCTCCAGTTTTTGCTGTTGTTCCAACCTCATTCCGCCTTCCACGCGGCGAACCAAGTATTTCCCCtt                                                                                                                                                                                                                           >  1:13179/1‑301 (MQ=255)
                                                                  ttgtgCGAGCCCACCTTTTCCTCCAGCGGTATATTCTTGCCGTTACGTTCTTTGAGCTTTGCCTGAACTTGGGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCttt                                                                                                                                                                                                                          <  2:87139/296‑1 (MQ=255)
                                                                     tgggagccacccgtttcgtgagcggggatatCCTTGCCGGAACGGTGTTTGAGGTTTGCCTGAACCTGTGGGTTGCTCAGGTCGGACTCTTTCAACTTCTGCTGCTCGCCCACGAGCATTGCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCa                                                                                                                                                                                                                       <  1:227664/273‑1 (MQ=255)
                                                                            ccaccGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAac                                                                                                                                                                                                                                                                     >  1:561492/1‑247 (MQ=255)
                                                                                      cTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGCTCCACTGc                                                                                                                                                                                                     >  1:473408/1‑301 (MQ=255)
                                                                                              gAATATTATTGCCGTAACGTTCTTTGCGTTTTGCCGGAACTTGTGGGCTGCTCAGGTCGACCTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTCAACCTTCCACACCACTTCCGAACAGTACTGGCTATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCgg                                                                                                                                                                                             <  1:423493/301‑1 (MQ=255)
                                                                                                            tAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGCTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTAAACGCGGCGAACCACGTATTTGCCCTTTTCACCATGGGCg                                                                                                                                                                                                              >  1:315609/1‑270 (MQ=255)
                                                                                                             atcGTTCTGTGGGCTTGTCCGTAACTGGTGGGCTCCTCAGGTCAAACTCCGTACACCTCTGCTGCTCGGCCACGAGCATTCCCCGCGCGCTCTGGTAAACCTGCCACACCACTTCCGAAACGTACTGGCGATCGTGTCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCg                                                                                                                                                                               <  2:8377/298‑1 (MQ=255)
                                                                                                                         gtgctggatgcaCTTTTGGGCTACGCAAGTCGTACTCTTTCAACTTACTCTGCTCGCACCCGCAGATTCACTGAGCGTTCTTGTACTCCTTCCACCCCCCTTCCGAACAGTACTGGAGATCGTCGCTCCATGAGAATCTGAAATCTTACGGATTACCGGAATAACGTTTTGCTGTTTGAGCCAGGTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTc                                                                                                                                                                  <  2:13179/290‑1 (MQ=255)
                                                                                                                             cgCCCGAGCTTGGTGGCGGCGCAGGTCGTACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGCTAAACCTTCCACACCACTTCCGAACAGTACTGGCGAGCGTCGCTCCATGAGCAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCaaa                                                                                                                                                               <  2:275911/299‑1 (MQ=255)
                                                                                                                                              tgctCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTGTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCCCTGCGTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGTCTTTtggtt                                                                                                                                             >  2:251519/1‑301 (MQ=255)
                                                                                                                                               gctCAGGTCGACCTCTTTCAACTTCTGCTGCCCGCCCCCGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCCCACCACTTCCGAACAGGACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTttgtt                                                                                                                                             <  2:473294/300‑1 (MQ=255)
                                                                                                                                                ctCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGACGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGTCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGGTTTTCACAGGGCCGACTGCGTCAAAAACGCAGGGCTTTTTGTTGc                                                                                                                                           >  2:351203/1‑301 (MQ=255)
                                                                                                                                                         aaCTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTAACGAGATACCGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCACTCCGCCTTCACCGCGGCGAACAACGTATTTGCCCTTTTAACCATGGGCGCACCAATGCGTGAGCGGGGTGTATTGCACCGGGCCGACTGCTTACAAAACGCAGGGCGTTTTGTTGCGttacaac                                                                                                                                   >  1:300965/1‑294 (MQ=255)
                                                                                                                                                                   aCTTCTTCTGCTCGCCCCCGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGc                                                                                                                                <  1:454107/293‑1 (MQ=255)
                                                                                                                                                                   aCTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCCTACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTTTATTTCCCCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTGTTTGTGTATCACCAGCagctcggc                                                                                                                        >  2:453996/1‑294 (MQ=255)
                                                                                                                                                                             ctCCCCCCCGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGAACGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCt                                                                                                                                                                    <  1:146714/247‑1 (MQ=255)
                                                                                                                                                                              tCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGTTCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTACGCCTTCAACGCCGCGAACAACGTATTTGCCCTTTTCACCATGCGCGTTCACCTGCTTGAGCGGGGTGTCTTTCACCGGGCCGACTGCTTCCAAAACGTAGGGATGTTTTTTGCGTTTCCCCAGCATACCGGTGTGGTtatact                                                                                                             >  1:341585/1‑299 (MQ=255)
                                                                                                                                                                              tCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAg                                                                                                                                                        >  1:353821/1‑258 (MQ=255)
                                                                                                                                                                               cgacaacgtacattaccagggagtcagagTAAACGTTCCACACCAATACCCACCACTACTAGGAATCGATGCTTCATGTGAAGCTGAAGACATACGTTTATCCGGAATAACGGTTTGCTGGTTGCGCAAGTGTTGGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACACCGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATc                                                                                                            <  2:52646/273‑1 (MQ=255)
                                                                                                                                                                                   cacaCGCGTTCCCAGCGCGTTCTGGTAACCCTTCCACACCACTTCCGCACAGTACTGGCGATTGTCGCTCCATTAGACGCTGAGATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCTCCAGTTTTTGCTGTTGTTCAACACTCATTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTa                                                                                                                                        <  1:321986/269‑1 (MQ=255)
                                                                                                                                                                                   cACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGcc                                                                                                                                                                                                                              <  1:551886/183‑1 (MQ=255)
                                                                                                                                                                                   cACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGcc                                                                                                                                                                                                                              >  2:551756/1‑183 (MQ=255)
                                                                                                                                                                                     cgcgCATTCCCAGCGCTTTCTGGTAAAGCTTCCCACCCCCTACCGAACAGAACTGGCGATCGTCGCTCCATGAGAAGCTGAACTCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTAAACACTCAGTCCGCCTTCAACGCGGCGAACAGCGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATg                                                                                                      <  2:34059/301‑1 (MQ=255)
                                                                                                                                                                                        gcATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCCCTGCTTGAGCGGGGTGTATGTCACCGGGCCGCCTGCTTCAAAAACGTAGGGCGTTGTGTTGCGTATCACTCGAATACCGGGGTggttttctacgctattggg                                                                                                   >  2:526093/1‑284 (MQ=255)
                                                                                                                                                                                            tCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCACAAACGTAGGGCTTTTTGTTGCGTAGCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGcc                                                                                               >  2:110787/1‑301 (MQ=255)
                                                                                                                                                                                                                                          tGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAg                                                                                                                                                        >  2:5012/1‑198 (MQ=255)
                                                                                                                                                                                                                                          tGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAg                                                                                                                                                        <  1:5012/198‑1 (MQ=255)
                                                                                                                                                                                                                                                     gCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCAACAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGCGCGAGCGATCTGAAAgaagggtttgcaggat                                      >  2:557661/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                    gAAATCATACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAgg                         <  1:251589/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gccgccgtgtgctcttccgctctcGTTTTGCTTTTTGCGCCAGTTTTTGCTGTTGTTAAACACTCCGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCaa                      <  1:220846/277‑1 (MQ=255)
                                                                                                                                                                                                                                                                       aacattccgTTTTCCGGGAGAACCGTTTGCTGTTTGCGCCAGTTTTTTCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGGGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAg                    <  1:198507/293‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aTACGGTTTACCGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTAAACGCGGCGAACAACGTATTTGCCCTTTTCCCCATGGGCGATCCACTGCTTGAGCGGGGTGGATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCa                 <  1:556789/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    ccGAGATAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGTTGCCAGGCAAAGGCAGGTAGAGaa        >  1:446155/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                      gagaTTATGTTTTGCTTTTTGTGCCAGTTTTTGCCTTTGTTAAACACTCATTCCGCCTTAAACGCGGCGAACAACTTATTTGCCCTTTTCACCATGGGCTATCCACTGCTTGAGCGGGGTGTATTTCCCCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAAcc     <  1:351292/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                         aTAACGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCTCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGAAGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACcagc  <  1:557792/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             cGTTTTGCTGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCGCCAGCATCCCGGTGTGGCTATTATCGGTATGGGTCGCCAGTTGGAGCGCTTTTCTTTGCGATGAGGGAGAGCTCTGAAAGATGATGTGGCTGGTTTGGGGCTGCCAGGCag                      >  2:220787/1‑276 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                
TATTACCACGAAAACAGCGGCCATTCTTTGGCGACCGTGGTGAGTTGTGGCGCATCGAAAACGGCTTGTGGCGAGACCACCGTTTCCTCAAGCGGAATATTCTTGCCGTAACGTTCTTTGAGCTTTGCCTGAACTTGTGGGCTGCTCAGGTCGAACTCTTTCAACTTCTGCTGCTCGCCCACGCGCATTCCCAGCGCGTTCTGGTAAACCTTCCACACCACTTCCGAACAGTACTGGCGATCGTCGCTCCATGAGAAGCTGAAATCATACGGTTTACCGAGATAACGTTTTGCCGTTTGCGCCAGTTTTTGCTGTTGTTCAACACTCAGTCCGCCTTCAACGCGGCGAACAACGTATTTGCCCTTTTCACCATGGGCGATCCACTGCTTGAGCGGGGTGTATTTCACCGGGCCGACTGCTTCAAAAACGTAGGGCTTTTTGTTGCGTATCACCAGCATACCGGTGTGGCTATAATCGGTATGGGTCGCCAGTTGGATCGCTTTACTTTGCGATGAGCGAGAGATCTGAAAGATGATGTCGCCGGTTTGTGGCTGCCAGGCAAAGGCAGGTACAGAAACCAGC  >  NC_000913/4127284‑4127865

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