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,2210CT94.7% 213.6 / 0.2 76E128K (GAA→AAA) yhiIputative membrane fusion protein YhiI
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (29/42);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTC  >  NC_000913/3629975‑3630464
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGcc                                                                                                                                                                                                                                                 <  1:212778/251‑1 (MQ=255)
  atatCTGCGGCAATGCTCCGTTCAGTGGCTTGTGCCGCTTCGACGCGTGTTTGCTCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATGCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTg                                                                                                                                                                                                                                               <  1:271886/251‑1 (MQ=255)
   gcccTGCGGCAATGCGCCGTTCATTTGCTGGTGCCGCCTCGACGCGGGTTTGCTCCTGAATGTTATTTGTTCTTGCCTCTTCAATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTg                                                                                                                                                                                                                                               <  2:289118/247‑1 (MQ=255)
         cggaaatgcgacgtcaaGTGGCTTGTGCCGCTTCGACGCGGGTTGGCGACTGAAGGATATGGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGa                                                                                                                                                                                                                                        <  2:355026/235‑1 (MQ=255)
                gcgcCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAAcc                                                                                                                                                                                                                                  <  1:230695/250‑1 (MQ=255)
                 ggcctttaaGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGGTTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTGTGCCTGGCGTTGATTAACCAg                                                                                                                                                                                                                                <  2:358057/243‑1 (MQ=255)
                     gtgcagtgggtggtgcgcctgcgaggcgGGTTTGCGCTTGAATTATATGGGTTCGTGCCGCTTCTATAGCCGCTTGAGAAGCCGATACCGGAGCTTTCGGCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCGCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTGTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGAc                                                                                                                                                                                                                            <  2:446105/226‑1 (MQ=255)
                        cAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACtgt                                                                                                                                                                                                                         <  1:46801/251‑1 (MQ=255)
                         aGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACtgtg                                                                                                                                                                                                                        <  1:317683/251‑1 (MQ=255)
                                tgtgCCGCTTCGACGCTGGTTTGCGCCTTAATGATAGTGGTGCTTTCCGCTTCTATAGCCGCTTTAGAAGCCGTGTCCTGAGCTTTCGCCGATTCCAGCGCAGCGCGGGCGCGCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCTAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCa                                                                                                                                                                                                                   <  1:64157/249‑1 (MQ=255)
                                tgtgCCGCTGCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCa                                                                                                                                                                                                                   <  2:362325/249‑1 (MQ=255)
                                     cGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGttt                                                                                                                                                                                                            <  1:360749/251‑1 (MQ=255)
                                          cGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGCTCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTGCGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCtt                                                                                                                                                                                                       >  2:48322/1‑251 (MQ=255)
                                          cGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTCCGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGCGTTTCGCtt                                                                                                                                                                                                       >  2:71063/1‑251 (MQ=255)
                                              gcgGGTTGGCGCCGTAACGATATTGGTGCGTGCGGCTTCGATAGCCGCGTGAGAAGTCGATTCCTGAGCTCTCGCCGATTCCTGCGCAGCTCGGGCGCTCGCATCGGCGGCGCGATCGTCAGCCAGCTGTTGCGTAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTGTTGCGTCGGAGTCCAGTTTTTCCTGGCGTTGATTAATCAGCGACTGTGCGGCACGAGGTTCTCTTTGt                                                                                                                                                                                                    <  2:291289/250‑1 (MQ=255)
                                              ccgGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTc                                                                                                                                                                                                   <  1:222274/250‑1 (MQ=255)
                                               cgGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTGCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCg                                                                                                                                                                                                  <  1:215837/251‑1 (MQ=255)
                                                        cccTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGAGACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTTATTAACCAGCGACTGGGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAg                                                                                                                                                                                         <  1:410808/250‑1 (MQ=255)
                                                            gAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTACAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCGGTGAACGGGAAAGCGTATGACGTTTTGCTACGGCGTCCAGTTTTTCCTGGCGTTGATTAAACAGCGACTGTGCGGCCCGAGTTTCGCTTTGTCGTTGCTCCAGcta                                                                                                                                                                                      >  2:254163/1‑248 (MQ=255)
                                                                ggtttttgtgcgtgcCGTTTCTATACCGGCTTTAGAAGCCGATACCTGAGATTTCGCCGATGCCAGCGCAGATCGGGCTCTCGCAGAGGCGGCGCGATCGTCAACCTGCTGTTGCGCAGAAATAGCCCCTCGGTGGTCCAGTGAACTGGAACGCGTATGACGTTTTGCTACGGCGTCCAGTGTTGCCTGTCGTTGATTAACCAGCGACTGTGCGGCACGATTTTCGCTTTGTCGTTGCTCCAGCAAAGcc                                                                                                                                                                                  <  2:470442/243‑1 (MQ=255)
                                                                       gggcgggcCGCTTCTATGGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCGCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGgc                                                                                                                                                                          <  2:223859/245‑1 (MQ=255)
                                                                             gcggtttcgagggcgtcgttggAAGCCGATACGTGAGCTTGCGCAGATTCCCGCGCAGCTCGGGCGCTCTCAGCGGCGGGGCGAGCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCGGCCTGGCGGTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCCTTGGCGTTTCTCCAGCATAGCCTGCGTGGCAGCAAc                                                                                                                                                                    <  2:319790/230‑1 (MQ=255)
                                                                             gcCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAAc                                                                                                                                                                    <  1:185525/251‑1 (MQ=255)
                                                                             gcCGCTTCTATAGCCGCTTGAGAAGCCGGTACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAAc                                                                                                                                                                    <  1:304907/251‑1 (MQ=255)
                                                                               cGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACgg                                                                                                                                                                  >  2:359033/1‑251 (MQ=255)
                                                                                      ataGCCGCTTTAGAAGCCTATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTtg                                                                                                                                                           >  1:121754/1‑251 (MQ=255)
                                                                                                 aGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTg                                                                                                                                                 <  1:434442/250‑1 (MQ=255)
                                                                                                      ccGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTgc                                                                                                                                           >  1:481411/1‑251 (MQ=255)
                                                                                                               taGCTTTGGCCGATGCCCGTGCCGCTCGGGCGCGCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTGCGCTTTGTCGTTGATCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGc                                                                                                                                   <  2:475708/249‑1 (MQ=255)
                                                                                                               gAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCt                                                                                                                                  >  1:245577/1‑251 (MQ=255)
                                                                                                                aGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCtt                                                                                                                                 <  1:395225/251‑1 (MQ=255)
                                                                                                                   tttCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGACCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGc                                                                                                                                                                                   <  1:484291/198‑1 (MQ=255)
                                                                                                                     tCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTGCGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAg                                                                                                                             <  1:48322/250‑1 (MQ=255)
                                                                                                                           aTTCCAGCGCAGCGCGGGCGCGCGCAGCGGCGGCGCGATCGTCACCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGATTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCTGTCGCTGt                                                                                                                      <  2:176555/251‑1 (MQ=255)
                                                                                                                                     aGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTgacacc                                                                                                            >  2:442695/1‑246 (MQ=255)
                                                                                                                                         cGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACAcgc                                                                                                         <  1:448081/250‑1 (MQ=255)
                                                                                                                                              gcTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTAt                                                                                                   >  1:412137/1‑251 (MQ=255)
                                                                                                                                                 ctcAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTGTTGCGACGGAGTCCAGTGTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTACGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCGTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATc                                                                                                  <  1:206630/249‑1 (MQ=255)
                                                                                                                                                   cAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCa                                                                                                <  2:402705/249‑1 (MQ=255)
                                                                                                                                                   cAGCGGCGGCGCGATCGTCATCCAGCGTTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGCCGCTGTTCCTGCAACACGCGAGTATCCAt                                                                                               <  2:367600/250‑1 (MQ=255)
                                                                                                                                                      cggcggcgCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTc                                                                                           <  1:26910/251‑1 (MQ=255)
                                                                                                                                                            cgcgATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGCTTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGGATCCATCTTCGCCAGc                                                                                     >  1:4673/1‑251 (MQ=255)
                                                                                                                                                            cgcgATCATCAACCACCTGTTGCGCAGAAATCGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTCTGACGTTTTGCTCCGGCGTACAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGCGTTTCTCTTTGTCGTTGCTCCAGCAAAGCCTGCACGGCAGCAACGGCGCTTTGTGCCTCTTTGCTATGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCCAGTATcaaccttctccacc                                                                                     >  2:237594/1‑238 (MQ=255)
                                                                                                                                                                      aTCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCtt                                                                           >  2:381300/1‑251 (MQ=255)
                                                                                                                                                                            cTGTTGCGCAGAAATAGCCCCTCGATGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCAGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTcccgc                                                                      >  2:190958/1‑246 (MQ=255)
                                                                                                                                                                              gTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAAc                                                                    >  2:150299/1‑250 (MQ=255)
                                                                                                                                                                              gTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGaa                                                                     <  1:212622/249‑1 (MQ=255)
                                                                                                                                                                                     aGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGc                                                            >  2:443390/1‑251 (MQ=255)
                                                                                                                                                                                        aaTAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCctt                                                         >  1:9306/1‑251 (MQ=255)
                                                                                                                                                                                           aGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCcttctt                                                      >  1:337426/1‑251 (MQ=255)
                                                                                                                                                                                                  cGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAg                                               <  1:146218/251‑1 (MQ=255)
                                                                                                                                                                                                  cGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAg                                               <  1:162533/251‑1 (MQ=255)
                                                                                                                                                                                                   gTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCAcc                                                 >  1:185056/1‑248 (MQ=255)
                                                                                                                                                                                                     tGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCAcc                                                 <  1:407106/246‑1 (MQ=255)
                                                                                                                                                                                                     tGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCAcc                                                 >  2:407106/1‑246 (MQ=255)
                                                                                                                                                                                                      gggCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGcgcg                                                                                                                                                                            <  2:100284/122‑1 (MQ=255)
                                                                                                                                                                                                      gggCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGcgcg                                                                                                                                                                            >  1:100284/1‑122 (MQ=255)
                                                                                                                                                                                                          cAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATggagt                                       >  1:300290/1‑248 (MQ=255)
                                                                                                                                                                                                             tGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCgc                                    >  1:402888/1‑250 (MQ=255)
                                                                                                                                                                                                              gAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGAt                                   <  1:392938/251‑1 (MQ=255)
                                                                                                                                                                                                                aCGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATAc                                 <  1:116819/251‑1 (MQ=255)
                                                                                                                                                                                                                     aaCGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTTTCGATACGCCCg                            <  1:150299/251‑1 (MQ=255)
                                                                                                                                                                                                                      aCGCGTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCgg                           >  1:31234/1‑251 (MQ=255)
                                                                                                                                                                                                                       cgcgTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGc                          <  1:360487/251‑1 (MQ=255)
                                                                                                                                                                                                                       cgcgTATGACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGAGAGAATCCATCTTCGCCAGCACTTCACCTTCGCGAACAACCTTGCCTTCTTTCACCAGAATGGAGTCGATACGCCCGGc                          >  1:234389/1‑251 (MQ=255)
                                                                                                                                                                                                                              gACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAAtttt                    <  1:95567/250‑1 (MQ=255)
                                                                                                                                                                                                                              gACGTTTTGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTg                   <  1:359033/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tttGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCGAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGc              >  1:13139/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tttGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGc              >  2:150335/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ttGCTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGc              >  1:273615/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cTACGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGACCAAACTTGCCTTCTTTCCCCAGAAGGGTGTCGATACTCCCGGCACTTTTGCCGGCAat           >  2:53215/1‑250 (MQ=255)
                                                                                                                                                                                                                                         aCGGAGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAata          <  1:450699/249‑1 (MQ=255)
                                                                                                                                                                                                                                             aGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACACACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACt    >  2:269712/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGCATGGTGTCGCTACGCCCGGCAATTTTGCTg                >  2:4875/1‑239 (MQ=255)
                                                                                                                                                                                                                                             aGTCCAGTTTTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTg                <  1:4875/239‑1 (MQ=255)
                                                                                                                                                                                                                                                ccAGTTTTGCCTGGCGTTTATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGGCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTc  <  2:194169/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CGATATCTGCGGCAATGCGCCGTTCAGTGGCTTGTGCCGCTTCGACGCGGGTTTGCGCCTGAATGATATTGGTGCGTGCCGCTTCTATAGCCGCTTTAGAAGCCGATACCTGAGCTTTCGCCGATTCCAGCGCAGCTCGGGCGCTCTCAGCGGCGGCGCGATCGTCATCCAGCTGTTGCGCAGAAATAGCCCCTCGTTGGGCCAGTGAACGGGAACGCGTATGACGTTTTGCTACGGAGTCCAGTTCTGCCTGGCGTTGATTAACCAGCGACTGTGCGGCACGAGTTTCGCTTTGTCGTTGCTCCAGCAAAGCCTGCGCGGCAGCAACGGCGCTTTGTGCCTCTTTGATTTGCGCGATGGCTTCCAGTCGCTGTTCCTGCAACACGCGAGTATCCATCTTCGCCAGCACTTCACCTTCGCGAACAAACTTGCCTTCTTTCACCAGAATGGTGTCGATACGCCCGGCAATTTTGCTGGCAATATCCACTTC  >  NC_000913/3629975‑3630464

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: