Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,156 G→A 18.9% intergenic (+65/‑82) purH → / → purD bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,490,1560GA18.9% 390.9 / 84.3 197intergenic (+65/‑82)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (70/79);  new base A (14/21);  total (84/113)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-02

Only 100 of 215 total aligned reads displayed.

GATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACC  >  NC_002947/5490007‑5490304
                                                                                                                                                       |                                                                                                                                                    
  tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                                    <  2:5387670/150‑1 (MQ=255)
      tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGACAACGAAGGTGGATTCAc                                                                                                                                                >  2:3441826/1‑150 (MQ=255)
       gATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAcc                                                                                                                                               <  2:6303503/150‑1 (MQ=255)
         tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgc                                                                                                                                             <  1:7828503/150‑1 (MQ=255)
         tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgc                                                                                                                                             <  2:6243870/150‑1 (MQ=255)
          gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                            <  2:1617129/150‑1 (MQ=255)
          gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                            <  2:326679/150‑1 (MQ=255)
              aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgtat                                                                                                                                        <  2:2555797/150‑8 (MQ=37)
               aGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATCCGAGGGGGGGGTACCCGCCGCCtt                                                                                                                                       >  2:7834743/1‑150 (MQ=255)
                 tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTc                                                                                                                                      <  1:1598103/150‑1 (MQ=14)
                 tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcg                                                                                                                                     <  1:4756537/150‑1 (MQ=255)
                  cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                    <  1:294762/150‑1 (MQ=37)
                  cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                    <  1:1413379/150‑1 (MQ=255)
                  cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                    <  1:5320083/150‑1 (MQ=255)
                   aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                   <  1:1085852/150‑1 (MQ=255)
                   aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                   <  2:6542779/150‑1 (MQ=255)
                   aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAAGTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgcg                                                                                                                                   <  2:4059509/150‑1 (MQ=25)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                  <  1:4566020/150‑1 (MQ=37)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                  <  1:3740925/150‑1 (MQ=37)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                  <  2:2329503/150‑1 (MQ=37)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                  <  2:3030040/150‑1 (MQ=37)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                  <  1:4959384/150‑1 (MQ=255)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                  <  1:3501845/150‑1 (MQ=255)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                  <  1:8520336/150‑1 (MQ=255)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                  <  2:2443434/150‑1 (MQ=255)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                  <  2:2898800/150‑1 (MQ=255)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                  <  2:6589221/150‑1 (MQ=255)
                    tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                  <  1:7683758/150‑1 (MQ=255)
                      gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCCAAAAAGGTGGTGGGTTGACCGCCCCATTCGCGGGc                                                                                                                                >  1:7300578/1‑150 (MQ=255)
                          ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCGCAAACGCTGGTGGGTTCGCCGCCGCCTTCGCGGGCCAGc                                                                                                                            >  2:1268905/1‑150 (MQ=255)
                          ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                            <  1:804776/150‑1 (MQ=255)
                          ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                            <  2:3155524/150‑1 (MQ=255)
                          ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                            <  1:804520/150‑1 (MQ=255)
                            gcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGccc                                                                                                                          <  2:5095102/150‑1 (MQ=25)
                                acgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTc                                                                                                                      <  2:2148509/150‑1 (MQ=255)
                                    ggCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGCATACGATGGGGGGTTGGACCCCCCAGTCGCGGGGCAGCCCGCGccccc                                                                                                                  >  2:5569400/1‑148 (MQ=255)
                                     gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGAAAAAGATGGTGGCTTGAACGCCCCCATGGCGGGGAAGCCCCCGCCcaca                                                                                                                 >  2:3417639/1‑150 (MQ=255)
                                     gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATTGTGGATTCACCGCCGCATTCGCGGGCCAGCCCCCGCCcacc                                                                                                                 >  2:8017805/1‑149 (MQ=255)
                                     gCCGGCAGGGCGTTGGTGTTCACGGGCTTGCGCAACGCCCGCCATTAATTAGGCGTAGCCCGCGGAAAGTCGTGTTCCTGTGGGCGGGGGCTTGCCCGCTAATACTTTGGTGGTTTCACCGCCGCATTCGGGGGCAAGCCCGCTCCcaca                                                                                                                 <  2:8537888/150‑1 (MQ=255)
                                           accgcgTTGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGCCAAGCCGAGATCCTGTGGGAGCGGGCTGGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTc                                                                                                           <  2:7129464/148‑1 (MQ=255)
                                           aTCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCCCTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTc                                                                                                           <  2:5831701/150‑1 (MQ=255)
                                                gATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCGCGCATACGATGGTGGTTTTACCCCCGCCTTAGCGGGGCAGCCCCCGCCCCCCGTGTTGGgcgc                                                                                                      >  2:6948225/1‑150 (MQ=255)
                                                  tGGTCTTCCCGGGCAGGCGCCACTCCCCCCACTAATCCCGGGTAGCCCGCGGAACGTCTCTACCGTGTGGGAGCGGGCTGGCCCGCTACTACGTTGGTGGATACCCCGCCGCTTGCGGGGGAACGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                                    <  2:1838181/150‑1 (MQ=255)
                                                  gggccTTCACGGCCCTGCCCCACGCCCGCCACTCCCTCCGCGGACCCCGCGGAAAGGCGACATCGGGGGGGCGGGGGCTTCCCCGCTAATCCGAGGGGGGATCCCCCGCCCAATCCGCGGGCACCCCCGCTCCCACAGTGTTCGGCGCaa                                                                                                    <  2:352845/146‑1 (MQ=255)
                                                       ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGCATAAGATTGTGGATTCACCGCCGCCATCGCCGGCCAGCCCCCACCCCCCGGGTTCGGCGCAAccccc                                                                                               >  2:4932462/1‑148 (MQ=255)
                                                        tCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGTCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCGAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACg                                                                                              >  1:7078128/1‑150 (MQ=255)
                                                            tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGAt                                                                                          <  2:1435255/150‑1 (MQ=255)
                                                                 tGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGGGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACtt                                                                                     <  2:5710862/150‑1 (MQ=255)
                                                                        cTTCCGCCACTAATTACGCGTATCCCGAGCCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGTTGGTGGATTCACCGCCGAATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                              <  2:4636013/150‑1 (MQ=38)
                                                                         ttCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCGACAAACCTTGGGGATTCCCCGCCGCAGTCGCGGGGCAGCCCGCGCCCACCCTGTTGGGCGCCAGCCCCCGGGTCCCATGAATCGAg                                                                             >  2:7031442/1‑150 (MQ=255)
                                                                         ttCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCCCGAATACGGTGGTGGATTAACCGCCGCAATAGCGGGCAAGCCCGCGCCCACAGTGGTCGGCGCAAGCCAACGGATAACATGAAACgcg                                                                             >  2:3592407/1‑148 (MQ=255)
                                                                              ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTTCCCGCCGAAACCATTGGGGGTTAAACGCCCCCATAGCGGGCAAGCCCGCTCCCACAGGGTTAGGCGCAACCCACCGGGTCAAGTGTAACGAGGtttt                                                                        >  2:2787024/1‑150 (MQ=255)
                                                                                   aaTTACGCGGATCCGGAGGCACGGCGAGATCCGGGGGGTGCGGGCTTGCCCGCGAATACGTTGGGGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAt                                                                   <  2:8357454/150‑1 (MQ=255)
                                                                                           ggATCCCGAGGCAAGTCGAGATCCTGTGGGTGGGGGCTTGCCCGCTAATACGATGGTGTATCCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGttt                                                           <  2:2870570/150‑1 (MQ=255)
                                                                                            gATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                          <  2:6517969/150‑1 (MQ=255)
                                                                                              tCCCGAGGCATTTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGTGGGTGTATTCACCGCCGCATTGGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                        <  2:1666062/150‑1 (MQ=38)
                                                                                                              gaTCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg                                        <  2:7877582/150‑1 (MQ=255)
                                                                                                                   gtgtgGGCGGGGGGTGGCCCGCCAATCCTTTGGTGGTTCCCCCGCCGAATTCGGGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                                    <  1:6835746/149‑1 (MQ=255)
                                                                                                                     gtgGGCGCGGGCTTGCCCGCTAATACGTTGGTGGATACCCCGCCGAATTCGCGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                  <  1:2782259/150‑1 (MQ=255)
                                                                                                                     gtgGGAGCGGGCTTGCCCGCTATTACGATGGTGGATTCCCCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                  <  2:7486615/150‑1 (MQ=255)
                                                                                                                     gtgGGAGCGGGCTTGCCCGCGAATACGTGGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                  <  2:8223009/150‑1 (MQ=255)
                                                                                                                         tcgCGGGCAAGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc                              >  2:2244089/3‑150 (MQ=255)
                                                                                                                         gggCGGGCTGGCCCGCGAATCCGTTGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCAc                              <  1:5784774/148‑1 (MQ=255)
                                                                                                                          cgCGGGCTTGCCCGCGAATACGAGGGTGGATCCACCGCCGAATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACg                             <  2:3933045/149‑1 (MQ=255)
                                                                                                                                  ggCCCGCTAATACGTGGGTGGATTCCCCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcc                     <  2:858172/149‑1 (MQ=255)
                                                                                                                                   gCCCGCTAATCCGCGGGTTCATCCCCCCCCGCTTTGGGGGGCACCCCCGCTCCCACAGTGTTCGGCACAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                    <  1:8293649/150‑1 (MQ=255)
                                                                                                                                   gCCCGCGAATCCGCTGGTGGATCCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                    <  2:1203231/150‑1 (MQ=255)
                                                                                                                                    cccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                   >  1:8318618/1‑150 (MQ=37)
                                                                                                                                     ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                  >  2:3360556/1‑150 (MQ=37)
                                                                                                                                     ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGAGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                  >  2:3323686/1‑150 (MQ=25)
                                                                                                                                     ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGGAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                  >  1:6342507/1‑150 (MQ=255)
                                                                                                                                     ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                  >  2:5606189/1‑150 (MQ=255)
                                                                                                                                     ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                  >  1:1225699/1‑150 (MQ=255)
                                                                                                                                     ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                  >  1:3065847/1‑150 (MQ=255)
                                                                                                                                     ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                  >  1:7356208/1‑150 (MQ=255)
                                                                                                                                     ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                  >  1:562841/1‑150 (MQ=255)
                                                                                                                                     ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                  >  1:6160541/1‑150 (MQ=255)
                                                                                                                                      cgcgAATGCGGCGGTGAATCCACCATCATATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                 >  2:2885116/1‑150 (MQ=37)
                                                                                                                                      cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCGGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                 >  2:2702946/1‑150 (MQ=255)
                                                                                                                                      cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                 >  2:2704288/1‑150 (MQ=255)
                                                                                                                                       gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                >  1:844193/1‑150 (MQ=37)
                                                                                                                                       gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                >  1:4659428/1‑150 (MQ=37)
                                                                                                                                       gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                >  1:8409872/1‑150 (MQ=255)
                                                                                                                                       gcgAACACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                >  1:6513192/1‑150 (MQ=255)
                                                                                                                                        cgAATACTCTGGTGGATTCCCCGCCGCATGCGGGGCCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa               <  2:1688991/150‑1 (MQ=255)
                                                                                                                                        cgAATACGTTGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa               <  2:3217386/150‑1 (MQ=255)
                                                                                                                                         gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg              >  1:2279449/1‑150 (MQ=255)
                                                                                                                                         gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg              >  1:5126534/1‑150 (MQ=255)
                                                                                                                                         gAATACGATGGTGGATTCACCGCCACATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg              >  2:7408706/1‑150 (MQ=255)
                                                                                                                                              aCGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc          >  2:4749087/1‑150 (MQ=255)
                                                                                                                                               ggtgggtggtTTCCCCGCCGCCGGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGACAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc         <  2:6714287/146‑1 (MQ=255)
                                                                                                                                                gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc        >  1:2636130/1‑150 (MQ=255)
                                                                                                                                                gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc        >  2:4418509/1‑150 (MQ=255)
                                                                                                                                                gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc        <  2:8661653/150‑1 (MQ=255)
                                                                                                                                                gATGGTGGATACACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc        <  1:6584005/150‑1 (MQ=255)
                                                                                                                                                   ggtgaATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAgg     >  1:7895332/1‑150 (MQ=39)
                                                                                                                                                    gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa    >  1:3918940/1‑150 (MQ=255)
                                                                                                                                                    gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa    >  1:8144835/1‑150 (MQ=255)
                                                                                                                                                     tggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAc   >  1:2969180/1‑150 (MQ=255)
                                                                                                                                                      ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAcc  >  2:8420306/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACC  >  NC_002947/5490007‑5490304

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