Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,152,000 A→T 7.5% intergenic (+33/‑364) PP_3654 → / → PP_3655 leucine‑responsive regulatory protein/cytosine/purine/uracil/thiamine/allantoin permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,152,0000AT 7.5% 403.0 / 14.5 134intergenic (+33/‑364)PP_3654/PP_3655leucine‑responsive regulatory protein/cytosine/purine/uracil/thiamine/allantoin permease family protein
Reads supporting (aligned to +/- strand):  ref base A (59/64);  new base T (4/6);  total (64/70)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

CCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCC  >  NC_002947/4151854‑4152148
                                                                                                                                                  |                                                                                                                                                    
ccGGCGTGTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCACTGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAC                                                                                                                                                   <  2:2982602‑M1/150‑1 (MQ=255)
  ggCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCC                                                                                                                                                 <  2:1257014‑M1/150‑1 (MQ=255)
    cGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGC                                                                                                                                               <  1:3252659‑M1/150‑1 (MQ=255)
    cGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTAGTCCCACTGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTACCCCCGC                                                                                                                                               <  2:1141436‑M1/150‑1 (MQ=255)
    cGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTAGTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGC                                                                                                                                               <  2:3962739‑M1/150‑1 (MQ=255)
      tttCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCgcgc                                                                                                                                             <  2:877432/150‑1 (MQ=255)
       ttCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaa                                                                                                                                            <  2:2840368/150‑3 (MQ=255)
       ttCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaa                                                                                                                                            <  1:4888299/150‑3 (MQ=255)
       ttCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaa                                                                                                                                            <  1:4888269/150‑3 (MQ=255)
       ttCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaa                                                                                                                                            <  1:3684473/150‑3 (MQ=255)
        tCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGcca                                                                                                                                           >  1:2076007/1‑149 (MQ=255)
        tCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCt                                                                                                                                           <  1:1055820/150‑1 (MQ=255)
         cAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaagg                                                                                                                                          <  1:131261/150‑5 (MQ=255)
         cAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTTACCcgcgaagg                                                                                                                                          <  1:185262/150‑5 (MQ=255)
         cAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTg                                                                                                                                          <  2:3539996/150‑1 (MQ=255)
              gCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAg                                                                                                                                     >  2:4996320/1‑150 (MQ=255)
                tGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGt                                                                                                                                   <  1:3228681/150‑1 (MQ=255)
                tGTATGACGAGCACCTGACCAGTAGTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGt                                                                                                                                   <  2:2864199/150‑1 (MQ=255)
                  tATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGa                                                                                                                                 <  1:3397248/150‑1 (MQ=255)
                   aTGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGAc                                                                                                                                <  2:5097238/150‑1 (MQ=255)
                    tGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACg                                                                                                                               >  2:387458/1‑150 (MQ=255)
                      aCGAGCTCCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTa                                                                                                                             <  1:4649463/150‑1 (MQ=255)
                       cGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAg                                                                                                                            <  1:3418501/150‑1 (MQ=255)
                        gAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCGGCAAAGGTGACGTAGc                                                                                                                           >  2:4757012/1‑150 (MQ=255)
                        gAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGAAAACCCCGCCGGCCAAGGTGACGTTGc                                                                                                                           >  1:4834956/1‑150 (MQ=255)
                             ccTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGCCCCCTTCGCGGGTAAAACCGCGCCTGCAAACGTGACGTAGCCCGTc                                                                                                                      >  1:4636180/1‑150 (MQ=255)
                              cTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttg                                                                                                                     <  1:1450687‑M1/150‑27 (MQ=255)
                              cTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCg                                                                                                                     >  2:668998/1‑150 (MQ=255)
                               tGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGc                                                                                                                           >  2:3076434/1‑143 (MQ=255)
                                gACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCa                                                                                                                   <  2:4427893/150‑1 (MQ=255)
                                gACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCa                                                                                                                   <  1:1919177/150‑1 (MQ=255)
                                   cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGa                                                                                                                <  1:3248774/150‑1 (MQ=255)
                                   cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGa                                                                                                                >  2:2421055/1‑150 (MQ=255)
                                   cAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGAAAACCCCCCCGGCAAAGGTGTCCTAGCCCCTCGCAGGa                                                                                                                >  1:1507943/1‑150 (MQ=255)
                                   cAGAATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGa                                                                                                                >  2:2421001/1‑150 (MQ=255)
                                      tATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCCCCAGCAAAGGTGTCGGAGCCCCGCGCCGGAgcg                                                                                                             <  1:1087493‑M2/150‑4 (MQ=255)
                                       agtccCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                            <  1:3035699‑M1/148‑36 (MQ=255)
                                       aTTCCCAATGTGACACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCCTGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                            <  2:1575956‑M1/150‑36 (MQ=255)
                                       aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCTCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                            <  2:2140661‑M1/150‑36 (MQ=255)
                                       aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                            <  2:2977523‑M1/150‑36 (MQ=255)
                                       aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                            <  2:4783325‑M1/150‑36 (MQ=255)
                                       aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgg                                                                                                            >  1:1178932‑M2/1‑146 (MQ=255)
                                       aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgg                                                                                                            >  2:2026418‑M2/1‑146 (MQ=255)
                                       aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgg                                                                                                            >  2:2542815‑M2/1‑146 (MQ=255)
                                       aTTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGGGCCGGCCTCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                            <  2:2238732‑M1/150‑36 (MQ=255)
                                        ttCCCACTGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGGTGCTCCCCATGTAGCTGTTTCCTCGTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcggg                                                                                                           >  2:3459729‑M2/1‑145 (MQ=255)
                                        ttCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTTGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcggg                                                                                                           >  2:3826030‑M2/1‑145 (MQ=255)
                                        ttCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcggg                                                                                                           <  2:3661661‑M1/150‑37 (MQ=255)
                                        ttCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcggg                                                                                                           <  2:3912324‑M1/150‑37 (MQ=255)
                                        ttCCCAACGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcggg                                                                                                           <  2:2985962‑M1/150‑37 (MQ=255)
                                            cAATGTGAAACGCTTGAGCTCGACCGTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCAGCAAACGTGACGTTGCCCGGCGCAGGAgcgggggtg                                                                                                       <  2:2654791‑M2/150‑10 (MQ=255)
                                             aaTGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgggtttac                                                                                                      >  2:1172677‑M2/1‑140 (MQ=255)
                                             aaTGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgggtttac                                                                                                      >  2:1333810‑M2/1‑140 (MQ=255)
                                               tgtgAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCCCCTGCAAAGGTGACGTAGCCCCGCGCAGGAgcgggggtaccc                                                                                                    <  2:3989104‑M2/150‑13 (MQ=255)
                                               tgtgAAACGCTTGAGCTCAACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgggtttaccc                                                                                                    <  1:3586097‑M1/150‑44 (MQ=255)
                                               agtgaAACGCTTGAGCTCGACCTTGGTGATGAAGGAACTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgggtttaccc                                                                                                    >  2:651064‑M2/2‑138 (MQ=255)
                                                   aaaCGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACACGCGCCGGCAAAGGTGACGTAGCCCCTCGCAGGAgcgggtttacccgcgc                                                                                                <  2:867777‑M2/150‑17 (MQ=255)
                                                    aaCGCTTGAGCTCGACCTTGGTGTTGAAGGAAGTGATCCAGGATCGGTTGCTCCCCCTGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgggtttacccgcgaa                                                                                               >  1:1183601‑M2/1‑133 (MQ=255)
                                                       gCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTAACACGCGCCTGGAAATGTGTGGTTTGCCGTGGGAGGTTGGGGTTTTTCCGCGAagag                                                                                            >  1:3607465/1‑150 (MQ=255)
                                                        cTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCCCCTGCAAACGTGACGTAGCCCCTCGCAGGAGCGGGTTTACCCGCGAAGAgg                                                                                           >  2:3675864/1‑150 (MQ=255)
                                                        cTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCAGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAgg                                                                                           >  2:3678351/1‑150 (MQ=255)
                                                        cTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCCCGCCGGCAAAGGTGGCGGAGCCCGTCGCAGGAGCGGGTTTAACCCCGAAGAgg                                                                                           >  1:2458894/1‑150 (MQ=255)
                                                              cTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGGGAAACCCCCCCCGCCAAGGTGTCGTAGCCCCCCCCAGGCGCGGGTTTTACCCCCCAGAggcccgcc                                                                                     >  1:1419025/1‑150 (MQ=255)
                                                               tCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                            >  2:4218211‑M1/1‑91 (MQ=255)
                                                               tCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgg                                                                                                            <  1:4218211‑M1/126‑36 (MQ=255)
                                                                  aCCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCAGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAgg                                                                                 >  1:2935750/1‑150 (MQ=255)
                                                                  aCCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAAACCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAgg                                                                                 >  2:2493000/1‑150 (MQ=255)
                                                                     ttGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCaa                                                                              <  2:3752789/150‑1 (MQ=255)
                                                                           atgaCGGAAGTGACCCCGGATCGGTGGCGCCCCATGTCGCTGTTTCCTCTTGTGCCGGCCCCTCCGCGGGTAAACCCGCGCCGCCAATGGGGACGTCGCCGGTCGCAGGCGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTc                                                                        <  2:174594/150‑1 (MQ=37)
                                                                           atgaAGTAAGTGAACAAGTTTGGGTTGCTCCCCTTGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAACCCCGCGCCTGAAATGGTGCCTTCGCCCGTCGCTGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTc                                                                        <  1:826974/150‑1 (MQ=37)
                                                                            ttaaggaagTGAGCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTGGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCa                                                                       <  2:3278360/148‑1 (MQ=255)
                                                                            tgaAGTATGTGATCCAGTATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCa                                                                       <  2:3485723/150‑1 (MQ=255)
                                                                                ggaagTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCACGCAACCGTTCATCTg                                                                   >  1:4646488/1‑150 (MQ=255)
                                                                                       aTCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCGGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCACGCAACCGTTCAAATGCCCTGaa                                                            >  2:1806452/1‑150 (MQ=39)
                                                                                        tCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGGTTACCCGCGAATAGGCCGGCCCAGGCAACCTTTCATCTGCCCAGaaa                                                           >  2:2098044/1‑150 (MQ=255)
                                                                                            ggttCGGTTGCTCCCCATGTAGCTTTTTCCTCTTGTGCCGGCCCCGTCGCGGTTAACCCCGCGCCTGCACCGGTGACGTAGCCCGTCGCTGGTGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAAcgc                                                       <  2:2249440/147‑1 (MQ=255)
                                                                                                   ttGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTTAACCCGCGCCGGCCAAGGTGACGTAGCCCCGCGCCGGAGCGGGGTTACCCCCGAAGAGGCCGGCCCCGGCCACCGTTCATCTGCCCGGGCAACGCGcacaac                                                >  1:1225724/1‑148 (MQ=255)
                                                                                                      cTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCa                                             <  2:4329179/150‑1 (MQ=255)
                                                                                                        ccccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCAAATGCCCTGAAAAAGCGCACATCGCAcc                                           >  2:363890/1‑150 (MQ=255)
                                                                                                         cccATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCTGCAAAGGTGGCGGAGCCCGTCGCCGGAGCGGGTTTACCCGCGAACAGGCCCGCCCACGCCACCCTTCAAAAGCCCTGAAAACGCGCCCATCGCACCa                                          >  2:3484872/1‑150 (MQ=255)
                                                                                                              gTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAATGGTGACGTAGCCCGTCGCTGGAGCGGTTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCa                                     <  2:1348443/150‑1 (MQ=255)
                                                                                                              gTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCCCATCGCACCAACCCa                                     >  2:1591085/1‑150 (MQ=255)
                                                                                                                   tGTTTCCTTTTGTGCCGGCCCGTGGGGGTTTAAACCGGCGCGTTCAATGTTGAGGTAGCCCGTCGCTGGTGGTGTTTAACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCa                                <  1:495085/150‑1 (MQ=255)
                                                                                                                    gTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCCCGCCGGCCAAGGTGACCGAGCCCCGCGCAGGAGCGGGGTTACCCGCGAAGCGGCCGGCCCCGGCAACCAATCAAATGCCCAGAAAACCCGCACCACCCACCAAcccccccccc                               >  2:2677399/1‑143 (MQ=255)
                                                                                                                     tttCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCAcc                              >  2:3299662/1‑150 (MQ=255)
                                                                                                                         ctctTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGa                          <  2:2077081/150‑1 (MQ=255)
                                                                                                                         ctctTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCATTCATCAGCCCTGAAAACGCGCACAACGCACCAACCCAGCCCACCGTGa                          >  2:4687169/1‑150 (MQ=255)
                                                                                                                          tctTGTGCCGGCCCCTTCGCGGGTTAAACCCCGCCAGCAAAGGTGGCGGAGCCCCGCGCAGGAGCGGGGTTAACCCCCCACAGGCCCGCCCACGCCACCAATAAACTGCCCCGAAAAAGCGCACAAAGCACCCAACCACcccccccggaa                         >  2:41752/1‑145 (MQ=255)
                                                                                                                          tctTGTGCCGGCCCCTTCGCGGGTAAAACCGCCCCAGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCCCGAAGAGGCCCGCCCAGGCAACCATTCATAATCCCGGAAAACGCGCACCAAGCACCAACCCCCCCCAcccagaa                         >  2:3543080/1‑145 (MQ=255)
                                                                                                                           ctTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATc                        <  2:4776513/150‑1 (MQ=255)
                                                                                                                            ttGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGGCGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCCACCTTTCATCTGCCCTGACAACGCGCACATCGCACCAACCCCGCCCACCGTGCTcc                       >  2:4224556/1‑150 (MQ=38)
                                                                                                                                       ccTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCAccc            >  2:1673870/1‑150 (MQ=255)
                                                                                                                                        cTTCGCGGGTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgt           >  2:2468035‑M1/1‑18 (MQ=255)
                                                                                                                                              gggTACCCCCGCtcctgcgcggggctacgtcacctttgcaggggggggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcaca     <  2:3989104‑M1/150‑139 (MQ=255)
                                                                                                                                              gggTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcagttccttcatcaccaaggtcgagctcaagcgtttcact     >  2:651064‑M1/1‑12 (MQ=255)
                                                                                                                                               ggTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCgccgcc    >  2:3537991/1‑150 (MQ=255)
                                                                                                                                                gTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacatt   >  2:1333810‑M1/1‑10 (MQ=255)
                                                                                                                                                gTAAACCCGCtcctgcgacgggctacgtcacctttgcaggcgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggtcgagctcaagcgtttcacatt   >  2:1172677‑M1/1‑10 (MQ=255)
                                                                                                                                                gTAAACCCGCGCCTGCAAAGGAGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCTCCAACCCAGCCCACCGTGATCCGGTTGCCACCCCcgccgcac   >  2:4883141/1‑148 (MQ=255)
                                                                                                                                                 tAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAgcgggtttacccgcgaaggggccggcacaagaggaaacagctacatggggagcaaccgatcctggatcacttccttcatcaccaaggttgagctcaagcgtttcacattg  >  2:3586097‑M2/1‑40 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCGGCGTTTCAGAAGCTGTATGACGAGCACCTGACCAGTATTCCCAATGTGAAACGCTTGAGCTCGACCTTGGTGATGAAGGAAGTGATCCAGGATCGGTTGCTCCCCATGTAGCTGTTTCCTCTTGTGCCGGCCCCTTCGCGGGTAAACCCGCGCCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCC  >  NC_002947/4151854‑4152148

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