Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,872,187 A→G 10.0% 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,1870AG10.0% 424.6 / 27.7 149D657G (GAT→GGT) ‡xdhBxanthine dehydrogenase subunit XdhB
Reads supporting (aligned to +/- strand):  ref base A (65/69);  new base G (4/11);  total (69/80)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

Only 100 of 149 total aligned reads displayed.

CGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGC  >  NC_002947/4872040‑4872334
                                                                                                                                                    |                                                                                                                                                   
cGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGAtt                                                                                                                                                   <  1:3644597/150‑1 (MQ=255)
 gCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTc                                                                                                                                                  <  2:571656/150‑1 (MQ=255)
 gCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTc                                                                                                                                                  <  1:2337443/150‑1 (MQ=255)
 gCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGAGGGCGATTc                                                                                                                                                  >  2:1255059/1‑150 (MQ=255)
   cGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGc                                                                                                                                                <  2:886249/150‑1 (MQ=255)
    gAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCt                                                                                                                                               <  2:4208927/150‑1 (MQ=255)
    gAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCt                                                                                                                                               <  2:1575564/150‑1 (MQ=255)
       gATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGaa                                                                                                                                            <  1:3865404/150‑1 (MQ=255)
           ttctACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCt                                                                                                                                        >  1:5043328/1‑150 (MQ=255)
               aCGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGa                                                                                                                                    >  1:3254713/1‑150 (MQ=255)
                cGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGGGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGAGTCGCTGAACCCTGCGag                                                                                                                                   >  2:3840784/1‑149 (MQ=255)
                cGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGAt                                                                                                                                   >  2:347753/1‑150 (MQ=255)
                 gACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATc                                                                                                                                  <  1:4844052/150‑1 (MQ=255)
                   ccGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGa                                                                                                                                >  1:5070891/1‑150 (MQ=255)
                     gCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGAGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGTACCCTGCGATCGACa                                                                                                                              >  2:692514/1‑150 (MQ=255)
                         ccAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCgg                                                                                                                          >  2:2613553/1‑150 (MQ=255)
                           aGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGcc                                                                                                                        <  2:891825/150‑1 (MQ=255)
                             gCCCGCGGCCGGCCGTTCTACTTCTTCTCCTTTGGCGCTGCCTGGGTGGAAGTTATCGTCGATACTCTGACCGGAGTATACAATTTGTTGCGCGACTACATCCTGAACTTCGTGGGCGATTCTCTGAACCCTGCTATCGACATCTGCCat                                                                                                                      <  2:3750104/150‑2 (MQ=255)
                             gCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAg                                                                                                                      <  2:4261645/150‑1 (MQ=255)
                              cccGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAgg                                                                                                                     <  1:4955439/150‑1 (MQ=255)
                               ccGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGt                                                                                                                    >  1:4040065/1‑150 (MQ=255)
                                cgcgGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTg                                                                                                                   <  2:4450914/150‑1 (MQ=255)
                                  cggccggcCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGa                                                                                                                 <  1:4971481/150‑1 (MQ=255)
                                  cggccggcCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGa                                                                                                                 <  2:1226805/150‑1 (MQ=255)
                                   ggccggccGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAg                                                                                                                <  1:1439986/150‑1 (MQ=255)
                                       ggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGggtg                                                                                                            >  1:3164814/1‑150 (MQ=255)
                                       ggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGggtg                                                                                                            >  1:1659092/1‑150 (MQ=255)
                                       ggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGggtg                                                                                                            >  2:1635843/1‑150 (MQ=255)
                                       ggccgTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGggtg                                                                                                            >  2:5068885/1‑150 (MQ=255)
                                       cgccgTTCTACGACTTCGCCTATGGCGCGGCCTGCGTGGAAGTGATCGTCGCTACGCGGCCCGGCGAATCCACGATGGTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCAGACCCCTGCGATCGACATCGGCCAGGTGGAGggtg                                                                                                            <  2:1907025/149‑1 (MQ=255)
                                              ctactaCTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCAt                                                                                                     >  2:1150578/1‑150 (MQ=255)
                                               tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                                                    <  2:42583/150‑1 (MQ=255)
                                               tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                                                    >  1:622673/1‑150 (MQ=255)
                                               tactacTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATc                                                                                                    >  2:4565970/1‑150 (MQ=255)
                                                actactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATcc                                                                                                   <  1:1000275/150‑1 (MQ=255)
                                                actactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATcc                                                                                                   >  2:110274/1‑150 (MQ=255)
                                                  tactTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAg                                                                                                 <  1:2970597/150‑1 (MQ=255)
                                                   attatgCCTATGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAgg                                                                                                <  1:4212150/145‑1 (MQ=255)
                                                   actTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAgg                                                                                                <  1:4034285/150‑1 (MQ=255)
                                                     ttCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGc                                                                                              >  1:1778793/1‑150 (MQ=255)
                                                     ttCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGc                                                                                              >  1:4398365/1‑150 (MQ=255)
                                                     ttCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGc                                                                                              >  2:413984/1‑150 (MQ=255)
                                                     ttCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGc                                                                                              >  1:1593008/1‑150 (MQ=255)
                                                         ccTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATgg                                                                                          <  2:1371299/150‑1 (MQ=255)
                                                           tttGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggc                                                                                        >  1:4316947/1‑150 (MQ=255)
                                                           tatGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGggc                                                                                        >  1:4884772/3‑150 (MQ=255)
                                                            atGGCGCGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGggct                                                                                       <  1:1449062/149‑1 (MQ=255)
                                                              ggCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctgg                                                                                     <  2:276532/150‑1 (MQ=255)
                                                               gcgcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctggc                                                                                    <  2:5043686/150‑1 (MQ=255)
                                                               gcgcGGCCTGCTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGggctggc                                                                                    <  1:2346997/150‑1 (MQ=255)
                                                                cgcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctggct                                                                                   <  1:4222061/150‑1 (MQ=255)
                                                                cgcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctggct                                                                                   <  1:2660438/150‑1 (MQ=255)
                                                                cgcTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGggctggct                                                                                   <  2:642366/150‑1 (MQ=255)
                                                                  cggcctgcTCGGAAGTATCGGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGa                                                                                 <  2:801568/148‑1 (MQ=255)
                                                                   ggcctgcTCGGAAGTATCTGTCGATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGac                                                                                <  2:3337755/149‑1 (MQ=255)
                                                                       tgcGTGGACGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                                                            <  2:3410497/150‑1 (MQ=255)
                                                                       tgcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGaccacc                                                                            >  2:198107/1‑150 (MQ=255)
                                                                        gcGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCg                                                                           <  2:3742509/150‑1 (MQ=255)
                                                                             gAAGTGATCGTCGATACGCTGACCGGCGAATACCAGATGTTGCGCGCCGACCTCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCCGGTGGAGGGTGGTTTCAGCCAGGGAATTGGCTGGCTGACCACCGAagag                                                                      <  2:4597913/150‑1 (MQ=255)
                                                                              aaGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGc                                                                     >  1:4431188/1‑150 (MQ=255)
                                                                                     cgtcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGtgtg                                                              <  1:2740517/150‑1 (MQ=255)
                                                                                      ggtcgATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGtctg                                                              <  2:3850024/149‑1 (MQ=255)
                                                                                        tcgATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGa                                                            >  2:3306091/1‑150 (MQ=255)
                                                                                          gATACGCTGACCGGCGAATATGTGGTGGAGCGCACCGATATTCTTCACGATACGGGCCGTTCGCTGAACCGTGTGATTGATATTGGTCAGGTGGAAGGCGGCTTCGTGCAGGGCATGGGCTGGCTGACCACGGAAGAGCTTGTCTgggat                                                          <  2:3478308/150‑4 (MQ=255)
                                                                                              cGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCa                                                      >  1:2478507/1‑150 (MQ=255)
                                                                                                     cGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCaa                                               >  1:1015006/1‑150 (MQ=255)
                                                                                                     cGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCaa                                               >  1:3059716/1‑150 (MQ=255)
                                                                                                      ggCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAg                                              <  1:837533/150‑1 (MQ=255)
                                                                                                       gCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGc                                             <  2:1674342/150‑1 (MQ=255)
                                                                                                       gCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGc                                             <  2:4472429/150‑1 (MQ=255)
                                                                                                       gCGAATACAAGAGGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGc                                             <  2:3969419/150‑1 (MQ=255)
                                                                                                           aTACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgat                                         <  1:1282876/150‑1 (MQ=255)
                                                                                                           aTACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgat                                         <  1:46484/150‑1 (MQ=255)
                                                                                                             aCAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCtgatga                                       >  2:3506861/1‑150 (MQ=255)
                                                                                                              cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc                                      <  2:835642/150‑1 (MQ=255)
                                                                                                              cAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAc                                      >  2:2543960/1‑150 (MQ=255)
                                                                                                               aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                     >  2:2185423/1‑150 (MQ=255)
                                                                                                               aaGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGAcc                                     >  1:2161798/1‑150 (MQ=255)
                                                                                                                   tGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcg                                 <  2:2354474/150‑1 (MQ=255)
                                                                                                                   tGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcg                                 <  2:2354539/150‑1 (MQ=255)
                                                                                                                    gTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGAGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcgg                                >  1:2466941/1‑150 (MQ=255)
                                                                                                                       gcgcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggccc                             <  2:3610609/150‑1 (MQ=255)
                                                                                                                         gcgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggcccgg                           >  1:2147508/1‑150 (MQ=255)
                                                                                                                          cgcCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAAcggcccggc                          <  1:519347/150‑1 (MQ=255)
                                                                                                                              gACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGc                      >  1:3254363/1‑150 (MQ=255)
                                                                                                                              gACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGc                      >  2:4031790/1‑150 (MQ=255)
                                                                                                                               aCATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCt                     >  1:106520/1‑150 (MQ=255)
                                                                                                                                cATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTa                    <  1:456076/150‑1 (MQ=255)
                                                                                                                                cATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTa                    >  1:4557788/1‑150 (MQ=255)
                                                                                                                                  tCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACa                  <  1:2977004/150‑1 (MQ=255)
                                                                                                                                     tGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAg                                                                                                 >  1:437832/1‑68 (MQ=255)
                                                                                                                                     tGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAg                                                                                                 <  2:437832/68‑1 (MQ=255)
                                                                                                                                     tGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGa               <  2:833533/150‑1 (MQ=255)
                                                                                                                                     tGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGa               <  1:1740272/150‑1 (MQ=255)
                                                                                                                                      gCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGAt              <  1:4885521/150‑1 (MQ=255)
                                                                                                                                        acgacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTTCAAGATcc            <  1:1255059/150‑1 (MQ=255)
                                                                                                                                          gacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCg          >  2:4815620/1‑150 (MQ=255)
                                                                                                                                          gacgTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCg          <  2:4862754/150‑1 (MQ=255)
                                                                                                                                                ggCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTg    >  2:2020261/1‑150 (MQ=255)
                                                                                                                                                  cGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGc  >  2:4408205/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCCGAAGATCTTCTACGACCGCAGCCAGGCCCGCGGCCGGCCGTTCTACTACTTCGCCTTTGGCGCTGCCTGCGTGGAAGTGATCGTCGATACGCTGACCGGCGAATACAAGATGTTGCGCGCCGACATCCTGCACGACGTGGGCGATTCGCTGAACCCTGCGATCGACATCGGCCAGGTGGAGGGTGGTTTCATCCAGGGCATGGGCTGGCTGACCACCGAAGAGCTGGTGTGGAACGCCAAAGGCAAGCTGATGACCAACGGCCCGGCCAGCTACAAGATCCCGGCCGTGGC  >  NC_002947/4872040‑4872334

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: