Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,586,587 A→C 18.8% R538R (CGT→CGG ydeP ← acid resistance putative oxidoreductase YdeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077791,586,5870AC18.8% 128.1 / 2.8 70R538R (CGT→CGGydePacid resistance putative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base A (16/38);  new base C (9/7);  total (25/45)
Fisher's exact test for biased strand distribution p-value = 7.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGA  >  NC_007779/1586353‑1586822
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
cATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATGTCAGCATTACACCGGCGGGTTTTAACACGCCAcgcg                                                                                                                                                                                                                                            <  1:1731234/239‑1 (MQ=255)
cATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGCCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAAcacgccacgc                                                                                                                                                                                                                                             <  2:863519/238‑1 (MQ=255)
cATCCAGCGCCTGTCAGCAGCTGCATTTATCAGGTGAAATCCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGATCTCTGGCAGCCCAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCCCCAGATACTCCCAGGCTGCCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACCCTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCAAcgcg                                                                                                                                                                                                                                            <  1:1343378/239‑1 (MQ=255)
                    cTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTCACACGCCCCGCGAGGCATGAATCATAGACAtt                                                                                                                                                                                                                       >  1:1871279/1‑239 (MQ=255)
                      gCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGa                                                                                                                                                                                                                     <  1:1335988/239‑1 (MQ=255)
                          ttATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATcc                                                                                                                                                                                                                 <  1:1880335/239‑1 (MQ=255)
                            aTCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATGGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTc                                                                                                                                                                                                               <  2:2115685/239‑1 (MQ=255)
                            aTCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTc                                                                                                                                                                                                               <  1:1604524/239‑1 (MQ=255)
                            aTCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACACCACGCGAGGCATGAATCATCGACATTGAATCCTc                                                                                                                                                                                                               <  2:600311/239‑1 (MQ=255)
                             tCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCa                                                                                                                                                                                                              <  1:726992/239‑1 (MQ=255)
                              cAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCaa                                                                                                                                                                                                             <  1:796365/239‑1 (MQ=255)
                                 gTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAAAACGCCACGCACGGCATGAAACATCGACATTTACACCTCAACgg                                                                                                                                                                                                          >  2:242347/1‑239 (MQ=255)
                                     aacccaccGGGATGCCGGTTGCGCTGGTTATGGTCGGCTAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACGCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTAc                                                                                                                                                                                                      <  2:1806523/236‑1 (MQ=255)
                                                 tGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTTACACCCCACCCGAGGCATGAAACATCCACCTTTAATCCTCCACGGGTACCGCCCgcgccc                                                                                                                                                                                           >  2:1069082/1‑236 (MQ=255)
                                                         gcgcTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTTACACCCCACCCCAGGCATGAATCATACAACATTAATCCTCCACGGGTACCCCCCGCGCACCCgtttttt                                                                                                                                                                                  >  2:739446/1‑239 (MQ=255)
                                                           gcTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATCCTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACCCTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGtttcttcg                                                                                                                                                                                 <  1:2247410/238‑3 (MQ=255)
                                                           gcTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTg                                                                                                                                                                                 <  1:868734/238‑1 (MQ=255)
                                                                   taGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAAt                                                                                                                                                                        <  1:261750/239‑1 (MQ=255)
                                                                      tCGGCAAACTCTGGCAGCACAGCTTCAATGTCATTGCGGATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCTTCGACATTGAATCCTCAACGGTTACCGCCTgcg                                                                                                                                                                                              <  1:2084707/214‑1 (MQ=255)
                                                                         gCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACCCCCCCCCCGGCATGAATCATCGACCATGAAATCTCAACGGTTTACCCCCGCGCACCCGTTTTTTGCATGTTAAACCCGc                                                                                                                                                                   >  2:1242253/1‑238 (MQ=255)
                                                                                 gggCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTg                                                                                                                                                                                 <  1:2041267/215‑1 (MQ=255)
                                                                                        ccaGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACCCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCCCACTCTGATTTCAGCATTACACCGGCGGGTTTTAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                                                                                                                                   <  1:1464172/238‑1 (MQ=255)
                                                                                        acaGCTTCATTGTAATTGCGAATGCTATCATAATCTTCCCCCAGGTACTCCCGGGCTACCCGGCTTTGGGGTAGTGCTGCCTGGGCGTTTCCGGCGGCCACCGCACACTCTGATTTCAGCATTCCACCGGCGGTTTTTAACCCGCCACGCGAGGCATGAATCATCGACATTGCATCCTCAACGTTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                                                                                                                                   <  2:1707964/239‑1 (MQ=255)
                                                                                        acaGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAc                                                                                                                                                      <  1:283383/236‑1 (MQ=255)
                                                                                        acaGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTAATTTCAGCATTACCCCGGCGGGTTTGAACACGCCCCGCGAGGCATGAATCTTCGACATTGAATCCTCACCGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAc                                                                                                                                                      <  1:1069082/236‑1 (MQ=255)
                                                                                        acaGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACGCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGTTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGAc                                                                                                                                                      <  1:1242253/236‑1 (MQ=255)
                                                                                         caGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTTACACGCCACGCGAGGCCTTAATCATCCACCTTGAATCCTCAACGGTTACCGCCCGCGCCCCCGTTTTTTTCTTGTTAAACTCCCCACGTCCCCCGAcccg                                                                                                                                                   >  2:474158/1‑236 (MQ=255)
                                                                                            cTTCAATGTCATTGCGAATGCGATCATAATCTTCCCCCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGCCCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTAAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                                                                                                                                   <  1:293262/235‑1 (MQ=255)
                                                                                             ttCAATGTCATGGGGAATGCGATCTTAATCTGCCCCCAGATCCGCCCAGGCTACCCCGCTCTGGGGTAGTGCGGCCTTCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGaa                                                                                                                                              <  2:1972895/239‑1 (MQ=255)
                                                                                                 aTGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTTACACCCCACGCGAGGCGTGAATCATCGACACTGAATCCTCAAAGGTTACCCCCCGCGCACCGTTTTTTTTCATGTCAATATCGCCACCGCCCAGGACCGGCATAAcata                                                                                                                                          >  2:571801/1‑239 (MQ=255)
                                                                                                 aTGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACGCCCCGCCAGGCATGAATCATCGACCTTGGATCCTCAACGGTTACCGCCCGGGCACCGTTTTTTTTCATTTCTAACTCCCTCCCCCCCACGACCGGCCGAAaata                                                                                                                                          >  2:2081413/1‑239 (MQ=255)
                                                                                                      aTTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGCATCCTCAACGGTTACCGCCTCCGCCCCGTTTTTTTTCATTTCAAACTCCCCACGCCCCCGGGCCCGCAGAATATAgtcat                                                                                                                                     >  2:1255147/1‑235 (MQ=255)
                                                                                                      aTTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACGCCACCCGCGGGATGAAACATCGAAAATGAATACTCAAAGGGTAACGCCCGCGCACCCTTTTTTTGGCTGTTAAATCCCCCACGTCCCACGGCCGGCAggaaattgccct                                                                                                                                     >  2:293262/1‑228 (MQ=255)
                                                                                                          cGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGa                                                                                                                                                                                                                                    >  1:2013180/1‑141 (MQ=255)
                                                                                                          cGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGa                                                                                                                                                                                                                                    <  2:2013180/141‑1 (MQ=255)
                                                                                                              tGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGTCCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCg                                                                                                                            <  1:2225939/239‑1 (MQ=255)
                                                                                                              tGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCCACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCCACGTCCCAGGACCGGCAGAATATTGCTATGCCGTGCg                                                                                                                            >  1:1219907/1‑239 (MQ=255)
                                                                                                              tGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACCCCCCCCGAGGCCTGAATCCTCCACACTGTATCCCCAACTGGTACCGCCCCCGCCCCGCTTTTTTTCCTGTCAATCTCGCCACATCCCCGGACCGGCCGAATATAGCTCTGGCggggg                                                                                                                            >  2:925509/1‑235 (MQ=255)
                                                                                                                cGATCATAATCTTCCACCAGATACTCCCAGGCTCCCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGt                                                                                                                          <  1:571801/239‑1 (MQ=255)
                                                                                                                        aTCTTCCCCCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGCTTCCCGCGCCCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACCCGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGAt                                                                                                                  <  1:739446/239‑1 (MQ=255)
                                                                                                                        aTCTTCCCCCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGAt                                                                                                                  <  2:671226/239‑1 (MQ=255)
                                                                                                                             ccaccaGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACCCCCCGCCAGGCATGAATCATCGACATTTAATCCTCAACGGCTACCCCCCGCGCACCGTTTTTTTGCATGTCAATTTCGCCACGTCCCCGGACCGGCCCAAAATAGATATGCCGTGCGGTCAACAGATGagag                                                                                                             >  2:869982/1‑239 (MQ=255)
                                                                                                                                ccaGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACTCGCCCCGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCgg                                                                                                                                                   <  1:1553020/199‑1 (MQ=255)
                                                                                                                                     tACTCCCGGGCTCCCACGCTTTGGGGTGGTGCTGCCTGGGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGtt                                                                                                        <  2:1219907/236‑1 (MQ=255)
                                                                                                                                       cTCCCAGGCTACCACGATCTGGGGTAGTGGTGCTCGCGCGATTCCCGCGACCACTGCCCGCTCTGATTGCAGCTTTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGGCATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCt                                                                                                   <  2:1461514/239‑1 (MQ=255)
                                                                                                                                         cccAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTa                                                                                                 <  1:869982/239‑1 (MQ=255)
                                                                                                                                          ccAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTTACACGCCACGCGAGGGCTGAATCATCGACACTGAATCCTCCACGGGTACCCCCTGCGCACCGTTTTTTTTCATGTCAATCTCCCCACGTCCCAGGACCCGCAGAATATAGGTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAg                                                                                                >  2:1940803/1‑239 (MQ=255)
                                                                                                                                              gCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACCCCCCGCGAGGCATGAAACATCGACATTTAAACCTCAAACGTTACCGCCCGCGCCCCGGTTTTTTTCATGTCAAACTCGCTACCTCCCCGGACCCGCAGAAAATATCTATGCCCTGCGGTCACAAGATTAGAGGGGTTAAGCCTAGTgg                                                                                             >  2:2104532/1‑238 (MQ=255)
                                                                                                                                               cTACCACGCCCTGGGGTAGGCCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTCCACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGc                                                                                            <  1:2140713/238‑1 (MQ=255)
                                                                                                                                                 aCCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACGCCACCCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCCGCGCACCGTTTTTTTTCATGTCAATCTCCCCACGTCCCCGGACCGGCCGAATATAGCTATGCCGTGCGGGCAACAAAAGAAAGCGGGTAAGGTTAAGGGCtc                                                                                          >  2:1002561/1‑237 (MQ=255)
                                                                                                                                                      gCTCTGGGTTAGTGCTGCCTGCCCGATTCCCGCGACCACGGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCCCGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTa                                                                                    <  2:746103/239‑1 (MQ=255)
                                                                                                                                                      gCTCTGGGGTAGTGCTCCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGTTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTa                                                                                    <  1:242347/239‑1 (MQ=255)
                                                                                                                                                                 gTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTTACACGCCACCCGAGGCATTAATCATCGGACTTGAATCCTCCACGGTTAACCCCCGCGCCCCGGTTTTTTTCATGTCAATATCCCTACCGCCCAGGCCCGGCCGAAAAAAACCATGTCCGCGGGTCTACACATGAGAGAGGGTGCGATCAGTGGGCCCGCGTACCCCCACATCt                                                                        >  2:1228522/1‑239 (MQ=255)
                                                                                                                                                                     tgcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAAtt                                                                    <  1:1255147/239‑1 (MQ=255)
                                                                                                                                                                     tgcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAAt                                                                     <  1:1806643/238‑1 (MQ=255)
                                                                                                                                                                      gcctgcGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACCCCCCGCGCGGGATGAATCATCGACATTGAATCCTCCACGGTTACCGCCTGCGCACCGCTTTTTTTCATGTCAATCTCGCTACGTCCCCGGACCGGCCGAAAATAGCTACTCCCTGCGGTCCAAACATGAGAGCGGGTAAGCTTTATGGCTACCTGTACCCCCCAATCTCAtt                                                                    >  2:393285/1‑238 (MQ=255)
                                                                                                                                                                            gcgATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCCCCCGCGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTTCATGTCAATCTCGCTACGTCCCCGGACCGGCAGAATAAAACTATGCCGTGCGGTCCACAGATGAGAGGGGTTAAGCTTAGTGGCCACGTGTACCCCCCAATCCAAATTCCGTTa                                                              >  2:1168175/1‑239 (MQ=255)
                                                                                                                                                                            ccgATTCCCGCGCCCCCTGCACCCTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTaa                                                             <  2:961175/238‑1 (MQ=255)
                                                                                                                                                                                        cccACTGCACACTCTGATTTCAGCATGCCACCGGCGGGTTTTACCCCGCCCCGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGcgct                                                 <  1:2081413/238‑1 (MQ=255)
                                                                                                                                                                                         cccctgccaacTCTGATTTCAGCATTACCCCGGCGGTTTTTAACCCGCCCCGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTCCCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGcgctc                                                <  1:2219574/230‑1 (MQ=255)
                                                                                                                                                                                          cACTGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACGCCCCGCGCGGCCTGAAACATCGACATTGAATCCTCAACCGTTACCGCCCCCGCACCGTTTTTTTGCATGTGAAACTCGCTACGTCCCAGGACCGGCAGAATATATCTATGCCGCGCGGTGAACAGATGGGAGCGGGTAAGCGTAGTGGCTACGTGTTACGCCCAAACTAAATGCGTTTACGGTACAgcgcca                                                >  2:71901/1‑236 (MQ=255)
                                                                                                                                                                                             tGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGAt                                                                                                                  <  1:37367/170‑1 (MQ=255)
                                                                                                                                                                                             tGCACACTCTGATTTCAGCATTACACCGGCGGGGTTTAACACCCCACCCGAGGCATGAATCATCGACCTTGAATCCTCCACGGGTACCCCCCGCGCCCCGTTTTTTTGCATGTTAATCCCGCTACGTTCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGGGGGTAAGGTTAGGGGGAACGTGTAACCCCAAAACTAAATGTGTTTACGGGTCAGcgcccgcc                                             >  2:1051959/1‑237 (MQ=255)
                                                                                                                                                                                                accctctGATTTCAGCATTACCCCGGCGGGTTTTAACACGCCCCGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAATTCTAATTGCGTTAACGGTACAGCGCTCGCTTcc                                          <  1:2063572/235‑1 (MQ=255)
                                                                                                                                                                                                   ctctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGAt                                      <  1:1002561/239‑1 (MQ=255)
                                                                                                                                                                                                   ctctGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGACATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGAt                                      <  1:474158/239‑1 (MQ=255)
                                                                                                                                                                                                   ctctGATTTCAGCATTACACCGGCGGGGTTTAACACCCCACCCCAGGCATTAATCATCGACATTGAATCCTCAACGGTTACCGCCCGCGCCCCGTTTTTTTTCATGTCAATCTCCCTACGTCCCAGGACCGGCCGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGGTAAGCTTTGTGGCTACGTGTACCCCCCAAACTAAATGCGTTAAAGGTACAACGCTCGCCTccccc                                       >  2:524505/1‑236 (MQ=255)
                                                                                                                                                                                                       gATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATTAATCATAGACATTGAATACTCAAACGTTACCGCCCGCGCACCGCTTTTTTGCATGTTCATCTCGCCTCGTCCCAGGACCGGCCGGATATAGCCATGCCCTGCGGGCAACAGAGGAGAgggggt                                                                                                        >  2:37367/1‑168 (MQ=255)
                                                                                                                                                                                                                         gcggGGTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCCACGGGTACCGCCTGCGCACCGCTTTTTTTCATGTCAAATCCCCCACATCCCCGGACCCGCCGAAAATAACTATTCCGTGCGGTCCACCGATGAGAGGCGGTAAGCTTAGTGGCTACGTGTCACGCCAAATCCAAATGTGCTTACGGGACATCGCTCGCTCCCCCACCTGGCCTTGCCAACTCaaaact                >  2:2017986/1‑237 (MQ=255)
                                                                                                                                                                                                                          cggGTTTTAACACGCCACGTGAGGAATAAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAAtt                <  1:1015978/238‑1 (MQ=255)
                                                                                                                                                                                                                               tttAACCCGCCCCGCGTGGAATGAATCATCGACATTGAATCCTCAACGGTTGCCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGcccc           <  2:194286/238‑1 (MQ=255)
                                                                                                                                                                                                                                       gccacgcGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGa  >  1:172079/1‑239 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                              
CATCCAGCGCCTTTCAGCAGCTGCATTTATCAGGTGAAAACCACCGGGATGACGGATGCGCTGGTTATAGTCGGCGAACTCTGGCAGCACAGCTTCAATGTCATTGCGAATGCGATCATAATCTTCCACCAGATACTCCCAGGCTACCACGCTCTGGGGTAGTGCTGCCTGCGCGATTCCCGCGACCACTGCACACTCTGATTTCAGCATTACACCGGCGGGTTTTAACACGCCACGCGAGGCATGAATCATCGACATTGAATCCTCAACGGTTACCGCCTGCGCACCGTTTTTTTGCATGTCAATCTCGCTACGTCCCAGGACCGGCAGAATATAGCTATGCCGTGCGGTCAACAGATGAGAGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGA  >  NC_007779/1586353‑1586822

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