Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,184,035 T→G 12.6% intergenic (‑123/+31) PP_4560 ← / ← csbD ribonuclease BN/stress response protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,184,0350TG12.6% 403.6 / 49.7 161intergenic (‑123/+31)PP_4560/csbDribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base T (92/43);  new base G (11/9);  total (103/58)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

Only 100 of 164 total aligned reads displayed.

AAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCT  >  NC_002947/5183887‑5184184
                                                                                                                                                    |                                                                                                                                                     
aaGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTc                                                                                                                                                      >  1:3393668/1‑150 (MQ=255)
 aGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCa                                                                                                                                                     >  2:27930/1‑150 (MQ=255)
   ccGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAcc                                                                                                                                                   >  2:2448611/1‑150 (MQ=255)
   ccGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAcc                                                                                                                                                   <  2:751100/150‑1 (MQ=255)
   ccGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAcc                                                                                                                                                   >  2:4032372/1‑150 (MQ=255)
    cgcgCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                  <  1:1892898/150‑1 (MQ=255)
    cgcgCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                                  >  1:2586886/1‑150 (MQ=255)
      cgcAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCga                                                                                                                                                <  1:3565192/150‑1 (MQ=255)
       gcAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                               <  1:1031786/150‑1 (MQ=255)
         aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaa                                                                                                                                             <  2:1655025/150‑1 (MQ=255)
         aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaa                                                                                                                                             <  1:2739473/150‑1 (MQ=255)
         aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaa                                                                                                                                             <  1:2748032/150‑1 (MQ=255)
         aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                             <  2:4024641/150‑1 (MQ=255)
         aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                             <  2:1168808/150‑1 (MQ=255)
           gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaaga                                                                                                                                           <  2:3742693/150‑1 (MQ=255)
           gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaaga                                                                                                                                           <  1:3785208/150‑1 (MQ=255)
           gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCACGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaaga                                                                                                                                           <  1:3807995/150‑1 (MQ=255)
              ggggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGc                                                                                                                                        <  2:3206781/150‑1 (MQ=255)
              ggggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGc                                                                                                                                        <  1:1093106/150‑1 (MQ=255)
               gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGGcc                                                                                                                                       <  1:642923/150‑1 (MQ=255)
               gggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGcc                                                                                                                                       <  2:1007087/150‑1 (MQ=255)
                ggAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCg                                                                                                                                      >  1:1392883/1‑150 (MQ=255)
                ggAAAATCAGAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGGCCg                                                                                                                                      <  1:1264242/150‑1 (MQ=39)
                      tCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCg                                                                                                                                >  2:3908770/1‑150 (MQ=255)
                         tAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCActtatt                                                                                                                             <  1:3635439/150‑6 (MQ=38)
                           gCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtca                                                                                                                           >  2:2316425/1‑150 (MQ=255)
                            cGTCTCCTTTCGCCGATGTATCTTTACACCGCCATGTCTAATAGACCGGGCCCGAGCTGTGCTGCTCCCGGCCGTGTTTTATTAATTTGCTTACTAATACGTGCGGGCCTCTTCGGGGGTTCACCCGATATGAGGCCGGCACCGCCtcaa                                                                                                                          <  2:2273105/150‑1 (MQ=255)
                             gTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaag                                                                                                                         >  1:568341/1‑150 (MQ=255)
                              tctcCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaagt                                                                                                                        >  1:3736363/1‑150 (MQ=255)
                                  cTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaagtcaag                                                                                                                    >  1:1999646/1‑150 (MQ=255)
                                   tttCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAgg                                                                                                                   >  1:2257339/1‑150 (MQ=255)
                                   tttCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAgg                                                                                                                   >  1:1568391/1‑150 (MQ=255)
                                        ccgcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttct                                                                                                              >  1:355151/1‑150 (MQ=255)
                                         cgcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttctt                                                                                                             >  1:2783719/1‑150 (MQ=255)
                                         cgcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttctt                                                                                                             >  1:2783829/1‑150 (MQ=255)
                                           cTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCa                                                                                                           >  2:1810985/1‑150 (MQ=255)
                                               aTCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGc                                                                                                       >  2:4198045/1‑150 (MQ=255)
                                                tCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCa                                                                                                      >  2:612005/1‑150 (MQ=255)
                                                       actggaaaGTATATTAGCCCGTGCCCTTGCTGTGCTGCTCCGGGCCGTTTTTTAATATTTTGCTTACTAATATTTGCCGGCCTCTTCGGGGTTTCACCCGAGAAGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTca                                                                                               <  2:523871/144‑1 (MQ=255)
                                                       acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTca                                                                                               >  2:3332258/1‑150 (MQ=255)
                                                       acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTca                                                                                               >  2:348864/1‑150 (MQ=255)
                                                           gCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACAtt                                                                                           >  1:587467/1‑150 (MQ=255)
                                                                 atAATAGACCGTGCCCTAGGTGTGCTGCTCCGGGCGGTGTTTTAATATTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGTTTCACCCTAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCtt                                                                                     <  2:1599286/149‑1 (MQ=255)
                                                                  tAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTg                                                                                    >  2:891008/1‑150 (MQ=255)
                                                                    aTAGACGGTGCCGGCGGTGTGGTGCGCCGGGCGGTGGTTGAAGATTTTGCTTACTATTACGTGCCGGCCTGTTCGGGGGTTCACCCGAGAAGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGAc                                                                                  <  1:3714882/150‑1 (MQ=255)
                                                                        aCCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACAGTTCCCTTGACctgc                                                                              >  1:1298838/1‑150 (MQ=255)
                                                                              cccGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcc                                                                        >  1:1044626/1‑150 (MQ=255)
                                                                              cccGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcc                                                                        >  1:2320152/1‑150 (MQ=255)
                                                                                cGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcctc                                                                      <  2:4000274/150‑1 (MQ=255)
                                                                                cGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcctc                                                                      >  2:3729726/1‑150 (MQ=255)
                                                                                     tgtgCTGCTCCCGGCCGTGGTTGAAGAATTTGCGTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCGGCTTCACCGCATCCTGCCCCTGTCCCTGGACCTGCTGCGCCTCGCCtt                                                                 >  1:830419/1‑150 (MQ=255)
                                                                                      gtgCTGCTCCCGGCCGTGTTTTAAGATTGTGCTTACGAATAAGTGCCGGCCTCTTCGCGGTTACACCCGAGAAGTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCttt                                                                <  2:2642773/150‑1 (MQ=255)
                                                                                          tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGc                                                            >  1:2289421/1‑150 (MQ=255)
                                                                                                cGGCCGTGGTTGAATAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTAACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTgc                                                      <  1:3710887/150‑1 (MQ=255)
                                                                                                 ggCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTgcg                                                     >  2:3678119/1‑150 (MQ=255)
                                                                                                   ccGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgcc                                                   >  2:4284621/1‑150 (MQ=255)
                                                                                                   ccGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTACACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgcc                                                   >  2:3792859/1‑150 (MQ=255)
                                                                                                     gTGGTTTACTAATTTGCTTACTAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgcctt                                                 <  1:2251249/150‑1 (MQ=255)
                                                                                                           gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTc                                           >  2:3466935/1‑150 (MQ=255)
                                                                                                           gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTc                                           >  2:2998222/1‑150 (MQ=255)
                                                                                                           gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTc                                           >  1:1829791/1‑150 (MQ=255)
                                                                                                            aagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCg                                          >  2:273926/1‑150 (MQ=255)
                                                                                                                 tttGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCg                                     >  1:888785/1‑150 (MQ=255)
                                                                                                                 tttGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGGCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCg                                     >  1:672707/1‑150 (MQ=255)
                                                                                                                      ttACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCt                                >  1:319743/1‑150 (MQ=255)
                                                                                                                          taaTAAGTGCCGGCCTCTTCGCGGGTACACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                            <  1:1845639/149‑1 (MQ=255)
                                                                                                                          gAATACGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                            <  2:2384835/150‑1 (MQ=255)
                                                                                                                          gAATAAGTGCCGGCCTCGTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                            <  2:2986921/150‑1 (MQ=255)
                                                                                                                          gAATAAGGGCCGGCCTCTTCGCGGGTTCACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                            <  1:1846173/150‑1 (MQ=255)
                                                                                                                          gAAAAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGACGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                            <  2:3222243/150‑1 (MQ=255)
                                                                                                                               aGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGt                       >  1:2596743/1‑150 (MQ=255)
                                                                                                                                gTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTc                      <  2:4319412/150‑1 (MQ=255)
                                                                                                                                gTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTc                      <  2:2573990/150‑1 (MQ=255)
                                                                                                                                 tGCCGGCCTCTTCGGGGTTACCCCCGATATGCGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                     <  1:2448611/150‑1 (MQ=255)
                                                                                                                                 tGCCGGCCTCTTCGGGGGTTCACCCGAGAAGCGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                     <  1:2957661/150‑1 (MQ=255)
                                                                                                                                 tGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCg                     >  2:2900097/1‑150 (MQ=255)
                                                                                                                                  gCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCgg                    <  2:1392883/150‑1 (MQ=255)
                                                                                                                                   ccGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGt                   >  1:196834/1‑150 (MQ=255)
                                                                                                                                     ggCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCa                 >  2:834616/1‑150 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCAc                >  1:4090945/1‑150 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCAc                >  1:3912264/1‑150 (MQ=255)
                                                                                                                                        ctctTCTCGGGTGAACCCGCGAAGAGGCCAGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCt              >  1:2145601/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCtt             >  2:1546722/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCtt             >  2:4052122/1‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACAGCATCCTTCACATTTCCCTTGACATCCTGCGCCTCGCCTATGAGCTCCTGAGCCTTGCCTACGGACCGCAGCTTTTCGATCTCGGTCACCtt             >  1:3480097/1‑150 (MQ=255)
                                                                                                                                           ttCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc           >  1:2030551/1‑150 (MQ=255)
                                                                                                                                           ttCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc           >  1:2194714/1‑150 (MQ=255)
                                                                                                                                           ttCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc           <  2:2194714/150‑1 (MQ=255)
                                                                                                                                           ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc           >  2:3505183/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc           >  1:4382183/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc           >  1:1152226/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc           >  1:1629393/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc           >  1:275238/1‑150 (MQ=255)
                                                                                                                                           aactcgGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc           >  1:1974835/5‑150 (MQ=255)
                                                                                                                                              gcgGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCaa        >  1:911971/1‑150 (MQ=255)
                                                                                                                                                gggTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAaca      >  1:3224803/1‑150 (MQ=255)
                                                                                                                                                  gTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACAcc    >  1:4158466/1‑150 (MQ=255)
                                                                                                                                                   ttCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                            <  1:1439419/125‑1 (MQ=255)
                                                                                                                                                   ttCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                            >  2:1439419/1‑125 (MQ=255)
                                                                                                                                                    tCACCCGAGAAGCGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCt  <  1:3335393/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGCCGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCT  >  NC_002947/5183887‑5184184

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: