Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,105,058 A→T 5.6% intergenic (+20/‑170) phhB → / → PP_5694 pterin‑4‑alpha‑carbinolamine dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,105,0580AT 5.6% 491.3 / 5.5 167intergenic (+20/‑170)phhB/PP_5694pterin‑4‑alpha‑carbinolamine dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (80/75);  new base T (3/9);  total (83/84)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

GCTGACTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  NC_002947/5104912‑5105207
                                                                                                                                                  |                                                                                                                                                     
gctgACTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCagag                                                                                                                                                    <  1:4861205/150‑1 (MQ=255)
 ctgactgaATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGTGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAgg                                                                                                                                                   >  1:2107796/1‑150 (MQ=255)
    actgaATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGccc                                                                                                                                                <  1:1358422/150‑1 (MQ=255)
      tgaATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTACACACCAGAAGGAAGGGGCCAGAGGCCCCt                                                                                                                                              >  1:4111772/1‑150 (MQ=255)
      tgaATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCt                                                                                                                                              >  1:1015234/1‑150 (MQ=255)
         aTGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCtttt                                                                                                                                           >  1:54943/1‑150 (MQ=255)
             ggCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCa                                                                                                                                       <  1:5212264/150‑1 (MQ=255)
               cAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATg                                                                                                                                     <  1:4969731/150‑1 (MQ=255)
                 aGGTCACCGTGACCTGGTGGAGCCACTCGCTCAAAGGCCTGCACCGCAACGCCTTCATCATGTGTGCCCGCACTGACAAGGTGGCGGAACCGGCTGACGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATggg                                                                                                                                   <  2:5428619/150‑1 (MQ=255)
                 aGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATggg                                                                                                                                   <  2:5341334/150‑1 (MQ=255)
                 aGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATggg                                                                                                                                   <  1:4197662/150‑1 (MQ=255)
                 aGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATggg                                                                                                                                   <  2:4688978/150‑1 (MQ=255)
                    tCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCa                                                                                                                                >  1:1143978/1‑150 (MQ=255)
                     cACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCat                                                                                                                               <  2:4522819/150‑1 (MQ=255)
                         gTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGAGAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATAtgt                                                                                                                           <  2:2286129/150‑1 (MQ=255)
                         gTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACTAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATAtgt                                                                                                                           <  2:4407557/150‑1 (MQ=255)
                             ccTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTgcgc                                                                                                                       <  1:1766786/150‑1 (MQ=255)
                              cTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCt                                                                                                                      >  2:2740429/1‑150 (MQ=255)
                              cTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCt                                                                                                                      >  2:2740419/1‑150 (MQ=255)
                                  tggAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTc                                                                                                                  <  1:1453135/150‑1 (MQ=255)
                                    gAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGc                                                                                                                <  1:4449947/150‑1 (MQ=255)
                                     aGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGcc                                                                                                               <  2:4059694/150‑1 (MQ=255)
                                      gCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGccc                                                                                                              <  1:806453/150‑1 (MQ=255)
                                       ccACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCg                                                                                                             >  2:5230393/1‑150 (MQ=255)
                                         aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCa                                                                                                           >  2:4777075/1‑150 (MQ=255)
                                         aCTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCa                                                                                                           >  1:197217/1‑150 (MQ=255)
                                            cGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGccc                                                                                                              <  1:2302001/144‑1 (MQ=255)
                                            cGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCaaaa                                                                                                        <  1:610393/150‑1 (MQ=255)
                                             gATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAAc                                                                                                       <  2:2181459/150‑1 (MQ=255)
                                             gATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAAc                                                                                                       <  2:1388895/150‑1 (MQ=255)
                                             gATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAAc                                                                                                       <  1:3636294/150‑1 (MQ=255)
                                                 aaaGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGAc                                                                                                   <  2:5367278/150‑1 (MQ=255)
                                                 aaaGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGAc                                                                                                   >  1:528776/1‑150 (MQ=255)
                                                  aaGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACa                                                                                                  >  2:3551740/1‑150 (MQ=255)
                                                   aGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACaa                                                                                                 <  2:3975573/150‑1 (MQ=255)
                                                     gCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAAc                                                                                               <  1:4777075/150‑1 (MQ=255)
                                                     gCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAAc                                                                                               >  1:5404405/1‑150 (MQ=255)
                                                       cTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGa                                                                                             >  2:1574089/1‑150 (MQ=255)
                                                       cTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGa                                                                                             >  1:3910683/1‑150 (MQ=255)
                                                            ccGCAACGACTTCATCTTGTGTGAACGCACTGACATGGTGGCGTAAACGGCTGATGGCGGTAATTAAACACCATATGTATGGGGCCTGAGGCCCCTTTTTTCTTGGGGCATATGTGGGCTTTTCGCCGGAAAAACCGACAAACGACCggg                                                                                        <  2:628440/150‑1 (MQ=255)
                                                            ccGAAACTACTTCTTATTGTTTGCACGCACTGACACGGGGGGGTAATCGGCTTATGGCGGGAATTAAACACCTGAAGTATGGGGCCTGTGGCCCATTTTTTATTGGTGATTTTGTGCGGTTTTCGCCCGCAAAACCGACAAACGACCggg                                                                                        <  2:628444/150‑1 (MQ=255)
                                                                     cTTCATCATGTGTGCACGCACTTACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACACACGACCGGGAAAAAAAcc                                                                               <  2:4185334/150‑1 (MQ=255)
                                                                      ttCATCTTGTGTGCACGCACTGACACGGTGGCGGAAACGGCTGAAGGCCGTAAGTAAACACCAGAAGGAAGGGGCCAGTGGCCCATTTTTTCATGGGGCATTTGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCa                                                                              <  1:3004586/150‑1 (MQ=255)
                                                                      ttCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCa                                                                              >  2:3916062/1‑150 (MQ=255)
                                                                      ttCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCa                                                                              >  1:5152764/1‑150 (MQ=255)
                                                                      ttCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCa                                                                              >  1:5152636/1‑150 (MQ=255)
                                                                      ttCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCa                                                                              >  2:3801682/1‑150 (MQ=255)
                                                                        catcctGTGTGCACGCAGTGACACGGTGGCGTAAACGGCTGACGGCCGGAATTAAACACCAGACGGATGGGGCCAGAGGCCCCTTTTTTCAGGGGGCATTTGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCaaa                                                                            <  1:2427254/150‑1 (MQ=255)
                                                                        catcatGTGTGCACGCACTTACATGGTGGGGGAAACGGCTGACGGCGGGAATTAAACCCCAGATGGAAGGGGCCAGAGGCCCCTTTTTTCTTGGGGAATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCaaa                                                                            <  2:480557/150‑1 (MQ=255)
                                                                        aatcatGTGTGACCGCACTGACACGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGAATATGTGCGCTGTTCGCCCGCAAAACCTACAAACGACCGGGAAAAAAACCaaa                                                                            <  2:2340188/149‑1 (MQ=255)
                                                                         atcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAc                                                                           >  2:2163975/1‑150 (MQ=255)
                                                                         atcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAAc                                                                           >  2:3079160/1‑150 (MQ=255)
                                                                          tcttgtgtgCACGCACTGACAAGGTGGCGTAAAGGGCTGACGGCCGTAAGTAAACACCAGAAGGACGGGGCCAGAGGCCCATTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACt                                                                          <  1:872241/147‑1 (MQ=255)
                                                                          tcatGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACt                                                                          >  1:3465298/1‑150 (MQ=255)
                                                                            atGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGt                                                                        >  1:11859/1‑150 (MQ=255)
                                                                                gtgCACGCACTGACACGGTGGCGGAAACGGCTTAAGGCCGTAAGTAAACACCAGATGGAAGGGGCCAGAGGCCCATTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATc                                                                    <  1:2125695/150‑1 (MQ=255)
                                                                                gtgCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATc                                                                    >  2:2578563/1‑150 (MQ=255)
                                                                                    acgcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGa                                                                >  2:194296/1‑150 (MQ=255)
                                                                                    acgcacTCAAAAGGTGCCTGAAACGCCTTAAGCCCGCAAGTAAACACCACAAGCAAGGGCCCACAGCCCCATTTTTTCATGGGGCATATGCGCGCTGTTCGCCCGCAAAACCCACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGa                                                                <  1:4798544/150‑1 (MQ=255)
                                                                                      gcacTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACt                                                              >  1:4236614/1‑150 (MQ=255)
                                                                                       cacTGACCAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGACGGAAGGGGCCAGAGGCCCATTTTTTCTTGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACtt                                                             <  2:3697316/150‑1 (MQ=255)
                                                                                         cTGACATGGTGGCGGAAACGGCTGAAGGCCGTATTTAAACACCATAAGTATGGGGCCTGCGGCCCCTTTTTTATTGGGGCATTTGTGCGCTGTTCGCCCGAAAACCCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                           <  1:5187846/150‑1 (MQ=255)
                                                                                         cTGACAGGGTGGCGTAAAGGGCTGATGGCCGGAATTAAACACCTTATGGATGGGGCCAGAGGCCCTTTTTTTCATGGGGAATATGGGCGCTTTTCGCCCCCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTtgt                                                           <  1:3658659/150‑1 (MQ=255)
                                                                                            acacGGTGGCGGAAACGGCTGACGGCCGGAAGTAAACACCAGACGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACACACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTa                                                        <  2:1961637/150‑1 (MQ=255)
                                                                                            aCAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTa                                                        <  2:1961622/150‑1 (MQ=255)
                                                                                              aaGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATc                                                      >  2:79532/1‑150 (MQ=255)
                                                                                                 gTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTg                                                   >  1:1168515/1‑150 (MQ=255)
                                                                                                      gtaaaCGGCTGAAGGCCGGAAGTAAACACAAGAAGGAAGGGGCCAGAGGCCCATTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCCACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGa                                              <  1:3012377/148‑1 (MQ=255)
                                                                                                      ggAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGa                                              >  2:4981775/1‑150 (MQ=255)
                                                                                                      ggAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGa                                              >  2:3233140/1‑150 (MQ=255)
                                                                                                      ggAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGa                                              >  1:4377500/1‑150 (MQ=255)
                                                                                                          aCGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtt                                          >  2:4213963/1‑150 (MQ=255)
                                                                                                          aCGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCtt                                          >  2:4213828/1‑150 (MQ=255)
                                                                                                             ggtGAAGGCCGTAAGTAAACACCAGAATGAAGGGTCCATCGGCCCCTTTTTTCTTGGGGCTTATGTGCGCTGTTCGCCCTCAAAACCGACAAACGACCGGGAAAAAAACCAATCTGTCATCCAGACTTGTGTATCCTGCCCTAGCttttt                                       <  2:1097193/148‑1 (MQ=255)
                                                                                                                taaGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGAAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGaa                                    <  2:302680/149‑1 (MQ=255)
                                                                                                                 aaGGCCGGATGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAg                                   <  1:4017428/150‑1 (MQ=255)
                                                                                                                   ggCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCt                                 >  1:4655786/1‑150 (MQ=255)
                                                                                                                    gCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTa                                <  2:2107796/150‑1 (MQ=255)
                                                                                                                     ccGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTGTTTGAAGCTAc                               >  2:918919/1‑150 (MQ=255)
                                                                                                                      cGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACa                              <  2:2574673/150‑1 (MQ=255)
                                                                                                                       ggAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACaa                             >  1:3578158/1‑150 (MQ=255)
                                                                                                                        gAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACaaa                            <  2:1015234/150‑1 (MQ=255)
                                                                                                                               acacCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCACTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCt                     >  1:4042094/1‑150 (MQ=255)
                                                                                                                                    aGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGaa                                    <  2:5142738/130‑1 (MQ=255)
                                                                                                                                    aGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGgcgc                >  2:3376345/1‑150 (MQ=255)
                                                                                                                                     gaaggaagGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGa                                                                <  2:4804327/101‑1 (MQ=255)
                                                                                                                                     gaaggaagGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCa               >  1:1608368/1‑150 (MQ=255)
                                                                                                                                       aggaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAACCGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGa             <  2:528776/150‑1 (MQ=255)
                                                                                                                                         gaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCaaaa                                                                                                        <  1:2439817/57‑1 (MQ=255)
                                                                                                                                         gaaggGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACt           >  1:5014820/1‑150 (MQ=255)
                                                                                                                                          aaggGGCCTGAGGCCACTTTTTTCTTGGGGCATATGTGCGCGGTGCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtt          <  1:2795323/150‑1 (MQ=255)
                                                                                                                                          aaggGGCCTCTGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtt          >  2:2869806/1‑150 (MQ=255)
                                                                                                                                            ggggCCAGCGGCCCCTTTTTTCTTGGGGCATATGGGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtttt        <  1:2003788/150‑1 (MQ=255)
                                                                                                                                            ggggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtttt        <  1:4470470/150‑1 (MQ=255)
                                                                                                                                            ggggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTATATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACtttt        <  2:2965501/150‑1 (MQ=255)
                                                                                                                                             gggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtc       >  2:3394970/1‑150 (MQ=255)
                                                                                                                                             gggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtc       <  2:252167/150‑1 (MQ=255)
                                                                                                                                              ggCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtca      <  1:4018005/150‑1 (MQ=255)
                                                                                                                                               gCCAGAGGCCCCTTTTTTCTTGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAACCGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcac     <  2:54943/150‑1 (MQ=255)
                                                                                                                                                  agagTCCCATTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCTGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTtcactca  <  2:4002303/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCTGACTGAATGGGGCAAGGTCACCGTGACCTGGTGGAGCCACTCGATCAAAGGCCTGCACCGCAACGACTTCATCATGTGTGCACGCACTGACAAGGTGGCGGAAACGGCTGAAGGCCGGAAGTAAACACCAGAAGGAAGGGGCCAGAGGCCCCTTTTTTCATGGGGCATATGTGCGCTGTTCGCCCGCAAAACCGACAAACGACCGGGAAAAAAACCAAACTGTCATCCAGACTTGTGTATCCTGCCCGAGCTTTTTGAAGCTACAAACGCGCCTGGCGCAGACTTTTCACTCA  >  NC_002947/5104912‑5105207

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: