Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,630,221 C→T E128K (GAA→AAA)  yhiI ← putative membrane fusion protein YhiI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,630,2210CT98.6% 227.8 / ‑6.6 70E128K (GAA→AAA) yhiIputative membrane fusion protein YhiI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/39);  minor base G (0/1);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTCCGT  >  NC_000913/3629982‑3630467
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
tGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCCGCGGCGGCGCGATCGTCATCCCGCTGTTGCGCAGAAATAGCCCCTCGTTTGGCCAGGGAACGGGAACGCGTATGACGGTTTGCTACGGAGTCCCGTTTTGCCTGGGGtt                                                                                                                                                                                                                                             >  2:223449/1‑251 (MQ=255)
ggCGGCAATGCGCCGTTCAGTGGCTTGGGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGTTCGTCATCCGGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGTGTATGACGTTTGGCTACGGAGTCCAGTTGTGCCTGGCGtt                                                                                                                                                                                                                                             <  2:137647/250‑1 (MQ=255)
   ggCAATGCGCCGTTCCGTGGCTTGGGCCGCTTCGCCGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAATCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGAt                                                                                                                                                                                                                                          <  1:103897/251‑1 (MQ=255)
   ggCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGGTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCCGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAACCCCTCGTTGGGCCAGTGAACGGGAACGCGTCTGACGTTTTGCTACGGTGTCCAGTTTTTCCTGGCGtttat                                                                                                                                                                                                                                          >  2:126993/1‑248 (MQ=255)
          cgcCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTTGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCCGCTGGTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGGGTCCAGTTTTGCCTGGCGTTGATTAAcc                                                                                                                                                                                                                                     >  1:33391/1‑249 (MQ=255)
           gcCGTTCAGTGGCTTGTGCCTCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGc                                                                                                                                                                                                                                  <  1:117401/251‑1 (MQ=255)
               ttCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCCGCGACt                                                                                                                                                                                                                              >  1:75824/1‑251 (MQ=255)
                tctgTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTGAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACt                                                                                                                                                                                                                              <  2:160663/247‑1 (MQ=255)
                tcattGGCTTGTGCCGCTTCGTCGCTGGTTTGCGCCTGAATGATATTGGGGCGTGCCGCTTCTATATCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGAGGCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTGATCCAGCTGTTGCGCAGAAATAGCCCCTCGTGGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACt                                                                                                                                                                                                                              <  1:29561/246‑1 (MQ=255)
                 cAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACtg                                                                                                                                                                                                                             <  2:54416/250‑1 (MQ=255)
                   gTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGc                                                                                                                                                                                                                          <  2:75824/251‑1 (MQ=255)
                     ggCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCg                                                                                                                                                                                                                         <  2:202808/250‑1 (MQ=255)
                         tgtgCCGCTTCTACGCGGTTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTGTAGAAGCCGATGCCTGAGCTTGCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTTTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTGTTGCTACGGAGTCCAGTTTTGCTTGGCGTTGATTAACCAGCGACTGTGCGGCAc                                                                                                                                                                                                                     <  2:232857/250‑1 (MQ=255)
                              cGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGttt                                                                                                                                                                                                               <  1:7067/251‑1 (MQ=255)
                                 ttCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCACTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGc                                                                                                                                                                                                            >  2:234471/1‑251 (MQ=255)
                                       gcgGGTTTGCGCCTGAAGGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGATTTTCGCTTTGTc                                                                                                                                                                                                      <  1:214568/251‑1 (MQ=255)
                                        cgGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATCCCTGAGCTTTCGCCGATTCCAGCGCAGCGCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGCTTCGCTTTGTCg                                                                                                                                                                                                     <  2:268203/251‑1 (MQ=255)
                                                 gccggAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAg                                                                                                                                                                                            <  2:30563/247‑1 (MQ=255)
                                                 gcCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAg                                                                                                                                                                                            <  1:60820/251‑1 (MQ=255)
                                                   cTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGc                                                                                                                                                                                           >  1:215383/1‑250 (MQ=255)
                                                                    gtgcCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGgcagca                                                                                                                                                                         <  1:216225/251‑1 (MQ=255)
                                                                         gCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGgc                                                                                                                                                                    <  1:230771/251‑1 (MQ=255)
                                                                             cTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCtt                                                                                                                                                                >  1:12080/1‑251 (MQ=255)
                                                                               ataGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCtt                                                                                                                                                                <  1:17886/249‑1 (MQ=255)
                                                                               ataGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTtg                                                                                                                                                              >  1:189778/1‑251 (MQ=255)
                                                                                taGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGCCATCCAGCTGTTGCGCCGAACTAGCCCCTCGTTGGGCCAGCGAACGGGAACGCGTATGACGCTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGATCCAGCAACGCCTGCGCGGCAGCAACGGCGCTTtgg                                                                                                                                                             >  1:178394/1‑250 (MQ=255)
                                                                                   ccGCTTTAGAAGCCGATACATGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGcc                                                                                                                                                          >  1:191494/1‑251 (MQ=255)
                                                                                       tttAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCtt                                                                                                                                                      >  1:18025/1‑251 (MQ=255)
                                                                                                       ggagtgttctccGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCCGCTGTTGCGCAGAAATAGCCCCTCGTGTGGCCAGTGAACGGGACCGCGTATGACGTGTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCTCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGAGCTTTGTGCCTCTTTGATTTGCGCTATGGc                                                                                                                                      <  2:246864/241‑1 (MQ=255)
                                                                                                        gAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATgg                                                                                                                                       <  2:218666/249‑1 (MQ=255)
                                                                                                               cGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGACCGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGGCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATTGCTTCCAGt                                                                                                                               <  2:221478/250‑1 (MQ=255)
                                                                                                               cGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTc                                                                                                                              >  1:86916/1‑251 (MQ=255)
                                                                                                               cGCCGATTCCAGCGCAGCGCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTc                                                                                                                              <  2:196025/251‑1 (MQ=255)
                                                                                                                              aGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGATTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCTCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAaca                                                                                                               <  1:35272/251‑1 (MQ=255)
                                                                                                                              aGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAac                                                                                                                >  2:89900/1‑250 (MQ=255)
                                                                                                                                    ggCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGCACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACAcgcg                                                                                                           >  1:118213/1‑249 (MQ=255)
                                                                                                                                      cgcTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGt                                                                                                        >  1:8601/1‑250 (MQ=255)
                                                                                                                                            cAGCGGCGGCGCGATCGTCCTCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATc                                                                                                 >  1:105821/1‑251 (MQ=255)
                                                                                                                                                            gTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCAc                                                                                 >  2:216025/1‑251 (MQ=255)
                                                                                                                                                            gTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCAc                                                                                 <  1:223449/251‑1 (MQ=255)
                                                                                                                                                                 ccAGCTGTTGCGCCGAAGTAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGCGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCACCACGCGAGTATCCATCTTCGCCAGCACTTCACCTTcg                                                                            >  2:221151/1‑251 (MQ=255)
                                                                                                                                                                 ccAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTc                                                                             <  2:215383/250‑1 (MQ=255)
                                                                                                                                                                 ccAGCTGTTGCGCAGAAATACCCCCTCGTTGGGCCAGTGACCGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTACCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTcg                                                                            >  2:256289/1‑251 (MQ=255)
                                                                                                                                                                   aGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCCCCTGcgcg                                                                          <  1:1002/251‑1 (MQ=255)
                                                                                                                                                                    gCTGTTGCGCAGAAATCGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGAttt                                                                                                                                                >  2:127511/1‑180 (MQ=255)
                                                                                                                                                                    gCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCTTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGAttt                                                                                                                                                <  1:127511/180‑1 (MQ=255)
                                                                                                                                                                    gCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTcgcg                                                                          >  1:171281/1‑250 (MQ=255)
                                                                                                                                                                         tGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCCGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCg                                                                                                                                                                                                     >  1:122552/1‑122 (MQ=255)
                                                                                                                                                                         tGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGGTTTGCTACGGAGTCCAGTGTGGCCTGGCGTGGATTAACCAGCGACTGTGCGGCACGAGGTTCGCGTTGTCg                                                                                                                                                                                                     <  2:122552/122‑1 (MQ=255)
                                                                                                                                                                          gcgcAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCGTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCATCACTTCACCTTCGCGAACAAAc                                                                   <  2:33391/251‑1 (MQ=255)
                                                                                                                                                                                           cGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCAcc                                                    >  2:131475/1‑249 (MQ=255)
                                                                                                                                                                                               gggCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATg                                              <  1:95516/251‑1 (MQ=255)
                                                                                                                                                                                                ggCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTGGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCTCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTGCACCAGAATgg                                             <  1:220851/251‑1 (MQ=255)
                                                                                                                                                                                                  ccAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGcgcg                                                                                                                                                                               >  2:164926/1‑119 (MQ=255)
                                                                                                                                                                                                  ccAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGcgcg                                                                                                                                                                               <  1:164926/119‑1 (MQ=255)
                                                                                                                                                                                                   cAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGgtgt                                          <  1:148501/251‑1 (MQ=255)
                                                                                                                                                                                                    aGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTc                                         <  1:216025/251‑1 (MQ=255)
                                                                                                                                                                                                     gTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCg                                        <  1:126993/251‑1 (MQ=255)
                                                                                                                                                                                                       gagcGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCATAGCCTGCGCGGCAGCAGCGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTGCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGAt                                      <  2:105821/248‑1 (MQ=255)
                                                                                                                                                                                                       gAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTGTTGCCTGGCGTTGATTAACCAGCGACTGGGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGAt                                      <  2:205883/251‑1 (MQ=255)
                                                                                                                                                                                                        aaCGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGACTGGCGTTGATTAACCAGCGACTGTGCGGCAAGAGTTTCGCTTTGTCGTTGCTCAAGCAAAGCCTGCGCGGCAGAAACGGCGCTTTGTGACTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGAt                                      >  2:256072/1‑250 (MQ=255)
                                                                                                                                                                                                           gggAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGc                                  >  2:65565/1‑251 (MQ=255)
                                                                                                                                                                                                            ggAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGcc                                 <  2:77830/251‑1 (MQ=255)
                                                                                                                                                                                                                   gTATGACGTTTTGCTACGGAGTCCAGTTTTGCCGGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAAt                          <  1:165876/251‑1 (MQ=255)
                                                                                                                                                                                                                       gACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTg                      >  1:108020/1‑251 (MQ=255)
                                                                                                                                                                                                                               atttcGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAGAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAGCACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGGCGATACGCCCGGCAATTTTGCGGGCAat              <  2:103942/247‑1 (MQ=255)
                                                                                                                                                                                                                                   cGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAat              >  1:103942/1‑247 (MQ=255)
                                                                                                                                                                                                                                    ggAGTCCAGTTTTGCCTGGCGTTGAATAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGGTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCa         <  2:171281/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCACTTTTGCTGGCAATATCCACtt      >  2:17301/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTCCGt  >  1:240280/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      
TGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTCCGT  >  NC_000913/3629982‑3630467

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