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,2210CT96.9% 194.9 / ‑2.3 66E128K (GAA→AAA) yhiIputative membrane fusion protein YhiI
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (28/35);  total (29/37)
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.51e-01

CGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTCCGT  >  NC_000913/3629975‑3630467
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
cGATATCGGCGGAAATGCGCCGTTCAGGGGCTGTGTCCGGTTCGACGCGGGTGTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCGTTAGAAGCCGATACCTGAGCTTTCGCCGAGTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGAGCGTCATCCAGCTGTTGCGCAGAAATAGCCCCGCGTTGGGCCAGTGAACGGGAACGCGTATGACGGTTTGCTACGGAGTCCAGTTTTGcc                                                                                                                                                                                                                                                     <  2:114642/251‑1 (MQ=255)
    atCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTgg                                                                                                                                                                                                                                                  <  1:386886/250‑1 (MQ=255)
     tCTGCGGCAATGCGACGTTCAGTGGCATTTGACGCTTCGCCGCGGGTTTGCGCCTGAATGTTATTTGGGCGTGCCGCCTCTATAGCCGCTTTTGAAGCCCAAAACTGAGCTTTCGCCGATTCCAGCGCAGCTCCGGCGCTCTCCGCGGCGGCGCGATCGTCAACCCGCTGTTGCGCAGAAATAGACCCCCGCTGGGCCAGGGAACGGGAACGCGTATGACGTTTTTCTCAGGAGGCCAGTTTTGCCTgggg                                                                                                                                                                                                                                                >  2:319926/1‑249 (MQ=255)
            cAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCCGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGGACGGGCACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTTCCTGGCGGTGAtt                                                                                                                                                                                                                                          >  1:327288/1‑250 (MQ=255)
                ggcgcctgtgaggggtttggtccGCTTCGACGCGGGTTTGCGTCGGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGAATCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTGGAGTAAcc                                                                                                                                                                                                                                      <  2:475086/230‑1 (MQ=255)
                 cgcCGTTCAGGGGCGTGTGCCGCTTCGACGCGGGTGTGCGCCTGAATGATATTGGGGCGTGCCGTTTCTATAGCTGCTTTAGAAGCCGATACCTGAGCTCTCGCCGATTCCAGCGCAGCGCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGATTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAg                                                                                                                                                                                                                                    <  2:492583/251‑1 (MQ=255)
                         aGTGTCTTGTTCCGCTTCTACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTGAGAAGCCGATACCTGAGCTTTCGCCTATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACtgtg                                                                                                                                                                                                                            <  1:62484/251‑1 (MQ=255)
                            ggCTTGTGCCGTTTCGTCGCGTGTGTGCGCCTGAATGATATTGGTGCGTGCCGATTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTGCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGTAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCgg                                                                                                                                                                                                                         <  1:455654/251‑1 (MQ=255)
                                     ccGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGtt                                                                                                                                                                                                                >  1:10212/1‑251 (MQ=255)
                                              cgcgGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAATCCCATACCTTAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCCGCTGTTGCGCAGACATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCATCGCCTGTGCGGCACGAGTTTCGCTTTGt                                                                                                                                                                                                       >  1:38733/1‑251 (MQ=255)
                                                   gTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTg                                                                                                                                                                                                  >  1:448986/1‑251 (MQ=255)
                                                         gcCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTcc                                                                                                                                                                                              >  1:105511/1‑249 (MQ=255)
                                                             gAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGAATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCaa                                                                                                                                                                                         >  2:485214/1‑250 (MQ=255)
                                                                         ggcgtgcCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGgca                                                                                                                                                                            <  1:132217/250‑1 (MQ=255)
                                                                          gcgggcCGCTTCTACAGCCGCTTTCGAAGCCGATACCTGAGCTTTCGCCGATGCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTTGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGgcag                                                                                                                                                                           <  2:157215/247‑1 (MQ=255)
                                                                           cgtgcCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAACGCCTGCGCGgcagc                                                                                                                                                                          >  2:125982/1‑251 (MQ=255)
                                                                            gtgcCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGGCGTTGCTCCAGCAAAGCCTGCGCGgcagc                                                                                                                                                                          >  1:5012/1‑250 (MQ=255)
                                                                               ccGCTTTTATGGCCGCTCTGGAAGCCGATACCTGACCTTTCGCCGATCCTAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCAACATCCCGCTGTTCCGCAGAACTCGCCCCTCGTTGCGCCAGTGAACGGGAAAGCGGATGCCGTTTTGATACGGCGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGCACCGAGTTTCGCTTTGTCGTTGCTCCGGCACAGCCTCCGCGGCAGCCACg                                                                                                                                                                      >  2:323322/1‑251 (MQ=255)
                                                                                 gCTTCTATAGCCGCTGTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGgc                                                                                                                                                                    <  1:322190/251‑1 (MQ=255)
                                                                                     cTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCACCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGgcgcat                                                                                                                                                                >  1:109430/1‑249 (MQ=255)
                                                                                             gCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCctc                                                                                                                                                        >  1:51180/1‑251 (MQ=255)
                                                                                              cTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCctct                                                                                                                                                       >  1:370485/1‑251 (MQ=255)
                                                                                                 tAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTg                                                                                                                                                    >  2:186424/1‑251 (MQ=255)
                                                                                                    aaGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGAtt                                                                                                                                                 >  1:117404/1‑251 (MQ=255)
                                                                                                      gCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTg                                                                                                                                               <  1:270633/251‑1 (MQ=255)
                                                                                                          atcccTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGAGCTCTCCGCGGCGGCGCGAACGTCAGCCAGCTGGTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGcgcg                                                                                                                                           <  2:200755/248‑1 (MQ=255)
                                                                                                           tACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGcgcg                                                                                                                                           >  1:387702/1‑250 (MQ=255)
                                                                                                                  gCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGCTCGTCATCCAGCTGTTGCGCAGACATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTc                                                                                                                                   >  2:23942/1‑251 (MQ=255)
                                                                                                                  gCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATATCCCCTCGTTGGGCCAGTTAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCATGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTc                                                                                                                                   <  1:305346/251‑1 (MQ=255)
                                                                                                                       cGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGt                                                                                                                               <  1:86695/250‑1 (MQ=255)
                                                                                                                       cGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTc                                                                                                                              <  1:225013/251‑1 (MQ=255)
                                                                                                                        gCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCGCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAggcg                                                                                                                             <  2:370485/251‑4 (MQ=255)
                                                                                                                          cGATTCCAGCGACGCCCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCt                                                                                                                           <  2:169473/251‑1 (MQ=255)
                                                                                                                               ccAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTTGGCCAGTGAACGGGACCTCGTATGACGTTTTGGTACGGAGTCCAGTTTGGTCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTcc                                                                                                                      <  1:98751/251‑1 (MQ=255)
                                                                                                                                 aGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTg                                                                                                                    <  1:104751/251‑1 (MQ=255)
                                                                                                                                          cGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACAcg                                                                                                             <  1:125982/249‑1 (MQ=255)
                                                                                                                                             gcgcTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAgg                                                                                                        >  1:460126/1‑250 (MQ=255)
                                                                                                                                              cgcTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTa                                                                                                       <  1:492775/251‑1 (MQ=255)
                                                                                                                                                    cAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCAt                                                                                                  <  1:61125/250‑1 (MQ=255)
                                                                                                                                                     aGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCt                                                                                                <  1:233228/251‑1 (MQ=255)
                                                                                                                                                       cggcggcgCGAGCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTGTGCCTGGCGTTGATTAACCAGCGACTGGGCGGCACGAGGTTCGCTTTGTCGTTGCTCCAGCAAAGCCGGCGCGGCAGCAACGGCGCGTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTc                                                                                              <  2:217903/251‑1 (MQ=255)
                                                                                                                                                        ggcggcgCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCg                                                                                             >  1:85596/1‑251 (MQ=255)
                                                                                                                                                          cggcgCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGcc                                                                                           <  1:73173/251‑1 (MQ=255)
                                                                                                                                                           ggcgCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGccc                                                                                          >  2:499374/1‑250 (MQ=255)
                                                                                                                                                                       aTCCAGCTGTTGCGCAGAAATAGCCCCCCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTGCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCtt                                                                              <  2:80874/251‑1 (MQ=255)
                                                                                                                                                                        tCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGCGTATCCATCTTCGCCAGCACTTAACCTTc                                                                             >  2:137473/1‑251 (MQ=255)
                                                                                                                                                                         ccAGCTGTTGCGCAGAAATCGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGATTGGCTACGGAGTCCAGTTTGGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTcg                                                                            <  2:112630/251‑1 (MQ=255)
                                                                                                                                                                         ccAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTcg                                                                            <  1:446376/251‑1 (MQ=255)
                                                                                                                                                                               gTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACa                                                                      <  1:101519/251‑1 (MQ=255)
                                                                                                                                                                                    gcAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACtt                                                                 >  1:491548/1‑251 (MQ=255)
                                                                                                                                                                                            aGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCcttctt                                                         >  1:472657/1‑251 (MQ=255)
                                                                                                                                                                                                ccTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAACCTTGCCTTCTTTccc                                                     >  2:142719/1‑249 (MQ=255)
                                                                                                                                                                                                        ggCCAGTGAACGGGAACGCGTATGACTTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATgg                                             <  2:32183/251‑1 (MQ=255)
                                                                                                                                                                                                          ccAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGgtg                                           <  1:373666/251‑1 (MQ=255)
                                                                                                                                                                                                             gTGAACGGGAAAGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGCGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGGTGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGCACAACCTTGCATTCTTTCACCAGCATTGTGTCg                                        >  2:63910/1‑251 (MQ=255)
                                                                                                                                                                                                               gAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGAt                                      <  1:271158/251‑1 (MQ=255)
                                                                                                                                                                                                                 aCGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCGTCTTTCACCAGAATGGTGTCGATAc                                    <  1:354134/251‑1 (MQ=255)
                                                                                                                                                                                                                  cGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACg                                   >  1:363709/1‑251 (MQ=255)
                                                                                                                                                                                                                     gAACGCGTATGTCGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGccc                                <  1:52700/251‑1 (MQ=255)
                                                                                                                                                                                                                              tGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCACtttt                       >  1:234397/1‑251 (MQ=255)
                                                                                                                                                                                                                               gACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTg                      <  2:10212/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tttGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCctt                                                            >  1:201786/1‑208 (MQ=255)
                                                                                                                                                                                                                                    tttGCTACGGAGTCCAGTTTTGCCGGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCctt                                                            <  2:201786/208‑1 (MQ=255)
                                                                                                                                                                                                                                     ttGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGc                 <  2:117404/250‑1 (MQ=255)
                                                                                                                                                                                                                                          aCGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAatat            <  1:186424/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   aGTTTTGCCTGTCGTTGATTAACCAGCAACTGTGCGACACGAGTTTCGCTTTATCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGAGCTTTGTGACTCTTTGATTTGCGCGATGGCTTCCAGTCGATGTTCCTGCAACAAGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGCAACGCCCGGAAATTTTGCTGGCAATATCCACTTCCGt  >  1:312586/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
CGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTCCGT  >  NC_000913/3629975‑3630467

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