Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,872,187 A→G 8.6% D657G (GAT→GGT) ‡ xdhB → xanthine dehydrogenase subunit XdhB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,872,1870AG 8.6% 943.1 / 53.0 313D657G (GAT→GGT) ‡xdhBxanthine dehydrogenase subunit XdhB
Reads supporting (aligned to +/- strand):  ref base A (148/138);  new base G (10/17);  total (158/155)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 315 total aligned reads displayed.

CACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCC  >  NC_002947/4872038‑4872335
                                                                                                                                                      |                                                                                                                                                    
    ccGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCg                                                                                                                                                  <  1:9010109/150‑1 (MQ=255)
      gAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCt                                                                                                                                                <  1:2296075/150‑1 (MQ=255)
           tcttctACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAAcc                                                                                                                                           <  2:3900780/150‑1 (MQ=255)
             ttctACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCt                                                                                                                                         >  2:872596/1‑150 (MQ=255)
             ttctACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCt                                                                                                                                         >  1:7405986/1‑150 (MQ=255)
             ttctACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCt                                                                                                                                         >  2:4763489/1‑150 (MQ=255)
               ctACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGc                                                                                                                                       <  1:126836/150‑1 (MQ=255)
                tACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCg                                                                                                                                      >  1:2378944/1‑150 (MQ=255)
                      cGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGAc                                                                                                                                >  1:6300398/1‑150 (MQ=255)
                             aGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGcc                                                                                                                         >  2:8636188/1‑150 (MQ=255)
                               gCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAg                                                                                                                       <  1:6128799/150‑1 (MQ=255)
                               gCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAg                                                                                                                       >  1:4202607/1‑150 (MQ=255)
                                 ccGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGt                                                                                                                     >  2:3193469/1‑150 (MQ=255)
                                      gccggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAgg                                                                                                                <  2:3604322/150‑1 (MQ=255)
                                      gccggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAgg                                                                                                                <  1:5006300/150‑1 (MQ=255)
                                      gccggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAgg                                                                                                                <  1:3821195/150‑1 (MQ=255)
                                         ggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGggtg                                                                                                             <  2:7204052/150‑1 (MQ=255)
                                          gccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGggtg                                                                                                             >  2:7543581/1‑149 (MQ=255)
                                          gccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGggtg                                                                                                             <  1:7543581/149‑1 (MQ=255)
                                            cgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGtt                                                                                                          <  2:3219360/150‑1 (MQ=255)
                                             gTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGttt                                                                                                         <  2:82976/150‑1 (MQ=255)
                                              ttCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTc                                                                                                        <  1:2905837/150‑1 (MQ=255)
                                              ttCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTc                                                                                                        <  2:6192552/150‑1 (MQ=255)
                                                 tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                                                     >  2:6358699/1‑150 (MQ=255)
                                                 tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                                                     >  1:4745619/1‑150 (MQ=255)
                                                 tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                                                     <  2:2341844/150‑1 (MQ=255)
                                                 tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                                                     <  2:8991969/150‑1 (MQ=255)
                                                  actactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATcc                                                                                                    >  1:332040/1‑150 (MQ=255)
                                                  actactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATcc                                                                                                    >  2:3033908/1‑150 (MQ=255)
                                                  actactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATcc                                                                                                    <  2:8007029/150‑1 (MQ=255)
                                                   ctactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCa                                                                                                   >  1:5841226/1‑150 (MQ=255)
                                                    tattatgCCTATGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAg                                                                                                  <  1:1079026/144‑1 (MQ=255)
                                                       ttCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGc                                                                                               >  2:1784679/1‑150 (MQ=255)
                                                       ttCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGc                                                                                               >  2:6006113/1‑150 (MQ=255)
                                                          gCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATg                                                                                            >  1:3194184/1‑150 (MQ=255)
                                                          gCCTATGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATg                                                                                            >  2:4380682/1‑150 (MQ=255)
                                                                ggCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctgg                                                                                      <  2:7965731/150‑1 (MQ=255)
                                                                ggCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGggctgg                                                                                      <  2:6770610/150‑1 (MQ=255)
                                                                ggCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGggctgg                                                                                      <  1:8578027/150‑1 (MQ=255)
                                                                 gcgcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctggc                                                                                     <  2:4739628/150‑1 (MQ=255)
                                                                   gcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGgctggctg                                                                                   >  2:7264614/1‑150 (MQ=255)
                                                                    ctgcctgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGa                                                                                  <  1:4948919/150‑1 (MQ=255)
                                                                       cctgcTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGacca                                                                               >  1:6636680/1‑150 (MQ=255)
                                                                        ctgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccac                                                                              >  2:2502554/1‑150 (MQ=255)
                                                                         tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                                                             >  2:3387591/1‑150 (MQ=255)
                                                                         tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                                                             >  1:2985206/1‑150 (MQ=255)
                                                                         tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                                                             <  2:4097920/150‑1 (MQ=255)
                                                                            gTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGaa                                                                          <  1:8919332/150‑1 (MQ=255)
                                                                            gTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGaa                                                                          >  1:2174833/1‑150 (MQ=255)
                                                                            gTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGaa                                                                          >  1:2062911/1‑150 (MQ=255)
                                                                             tGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAag                                                                         >  1:7295968/1‑150 (MQ=255)
                                                                                aaGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGc                                                                      <  2:6542291/150‑1 (MQ=255)
                                                                                aaGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGc                                                                      >  1:870060/1‑150 (MQ=255)
                                                                                    tatcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCttg                                                                   <  2:466748/145‑3 (MQ=255)
                                                                                     atcggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTgt                                                                  <  2:226803/146‑1 (MQ=255)
                                                                                     aTCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtg                                                                 >  1:4836630/1‑150 (MQ=255)
                                                                                      tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGa                                                                                  >  1:1183611/1‑132 (MQ=255)
                                                                                      tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                                                                <  1:5197165/150‑1 (MQ=255)
                                                                                      tcgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGgtgt                                                                <  1:100998/150‑1 (MQ=255)
                                                                                         gtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTgg                                                              <  1:5706503/150‑1 (MQ=255)
                                                                                         gtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTgg                                                              >  1:6407132/1‑150 (MQ=255)
                                                                                           cgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTggga                                                            <  1:8121797/150‑3 (MQ=255)
                                                                                           cgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGaa                                                            <  1:7966632/150‑1 (MQ=255)
                                                                                            gATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAAc                                                           <  1:4398215/150‑1 (MQ=255)
                                                                                              tACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGc                                                         <  1:6923063/150‑1 (MQ=255)
                                                                                              tACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGc                                                         <  1:5527107/150‑1 (MQ=255)
                                                                                              tACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGc                                                         <  2:1532399/150‑1 (MQ=255)
                                                                                                 gCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAagag                                                                       >  2:7140149/1‑133 (MQ=255)
                                                                                                 gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                                                      >  2:1273190/1‑150 (MQ=255)
                                                                                                 gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                                                      <  2:5601442/150‑1 (MQ=255)
                                                                                                 gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                                                      >  2:1678681/1‑150 (MQ=255)
                                                                                                 gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                                                      >  2:2993028/1‑150 (MQ=255)
                                                                                                 gCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaa                                                      >  1:7057959/1‑150 (MQ=255)
                                                                                                  cTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaaa                                                     <  2:2558482/150‑1 (MQ=255)
                                                                                                  cTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaaa                                                     <  1:2363578/150‑1 (MQ=255)
                                                                                                  cTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCaaa                                                     <  1:3530184/150‑1 (MQ=255)
                                                                                                   tGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTGGGATGACAAAg                                                    >  1:265132/1‑150 (MQ=255)
                                                                                                   tGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAg                                                    <  1:1626203/150‑1 (MQ=255)
                                                                                                     aCCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTGGGATGACAAAgga                                                  <  2:1051847/150‑2 (MQ=255)
                                                                                                           gAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtg                                            >  2:4203661/1‑150 (MQ=255)
                                                                                                                cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc                                       >  2:5018957/1‑150 (MQ=255)
                                                                                                                cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc                                       >  1:1855955/1‑150 (MQ=255)
                                                                                                                cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc                                       >  2:3693358/1‑150 (MQ=255)
                                                                                                                cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc                                       >  2:6250591/1‑150 (MQ=255)
                                                                                                                 aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAATGGCAAGCTGATGAcc                                      >  1:1455973/1‑150 (MQ=255)
                                                                                                                 aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                      >  2:846153/1‑150 (MQ=255)
                                                                                                                 aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                      >  2:3296867/1‑150 (MQ=255)
                                                                                                                 aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                      >  1:2537062/1‑150 (MQ=255)
                                                                                                                 aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                      >  2:4038036/1‑150 (MQ=255)
                                                                                                                 aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                      >  1:3114530/1‑150 (MQ=255)
                                                                                                                  aGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCa                                     >  1:7629960/1‑150 (MQ=255)
                                                                                                                     tGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcg                                  >  2:8552975/1‑150 (MQ=255)
                                                                                                                           gcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggcccgg                            >  2:3535228/1‑150 (MQ=255)
                                                                                                                           gcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggcccgg                            >  2:3436673/1‑150 (MQ=255)
                                                                                                                                     ccTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACaa                  <  2:537264/150‑1 (MQ=255)
                                                                                                                                      cTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAg                 >  2:1155517/1‑150 (MQ=255)
                                                                                                                                       tGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGa                >  1:7694251/1‑150 (MQ=255)
                                                                                                                                         cACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCAGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATc              <  2:5232686/150‑1 (MQ=255)
                                                                                                                                             acgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCgg          <  2:6975209/150‑1 (MQ=255)
                                                                                                                                                tGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCg       >  2:6812771/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CACGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGCC  >  NC_002947/4872038‑4872335

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