Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,854,478 G→A 6.4% intergenic (‑182/+115) recR ← / ← ybaB recombination protein RecR/nucleoid‑associated broad specificity regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,854,4780GA 6.4% 186.3 / 2.7 62intergenic (‑182/+115)recR/ybaBrecombination protein RecR/nucleoid‑associated broad specificity regulator
Reads supporting (aligned to +/- strand):  ref base G (44/14);  new base A (2/2);  total (46/16)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCA  >  NC_002947/4854329‑4854627
                                                                                                                                                     |                                                                                                                                                     
ttATACCCCGGGGGGGCGGTGGCGTCAAATCCCCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGGGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTcgcg                                                                                                                                                       <  2:5322304/150‑1 (MQ=255)
            ggggggggggCTTCAAAACCCCCCCCCAAACGCCGGAAGGTTCTGCTTTTTTGTGGGCGCGGGTTTGTGCCGCTAAGGGGGTGGCACCCCGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGa                                                                                                                                           <  2:1480614/150‑1 (MQ=255)
                   gggcGTCAAATCACCGCGCCCAACGCCTGCATGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAACGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaagg                                                                                                                                    >  2:4510817‑M2/2‑132 (MQ=255)
                                                                         gCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTATGTGCTGCGAAGCTGAAAAACCGGGGCCGCTTCGCGCcccattagcgacacaaggccgctcccacaaaaaatcagaccatgcacgctttgggcgaggtgatttgacgcc                                                                              <  1:4462073‑M2/150‑73 (MQ=255)
                                                                           cTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCg                                                                            >  2:3966051/1‑150 (MQ=255)
                                                                           cTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACGGGGTTGACATCTCCGCAAAAAAAACCg                                                                            >  2:3202911/1‑150 (MQ=255)
                                                                                   gcgAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctgtgggcgaggtgatttgacgccaccgcccccg                                                                    <  2:4108413‑M2/150‑83 (MQ=255)
                                                                                     gAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgct                                                                  <  1:73325‑M2/150‑85 (MQ=255)
                                                                                       cggggCTGCTTCGCAGCCCTGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCa                                                                >  2:1058205/3‑150 (MQ=255)
                                                                                         gggCTGCTTCGCAGCCCTGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgc                                                              >  2:351706/1‑150 (MQ=255)
                                                                                          ggCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggta                                                             <  1:5328530‑M2/150‑90 (MQ=255)
                                                                                           gCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggtat                                                            <  1:1657697‑M2/150‑91 (MQ=255)
                                                                                           gCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACCAGGCCGCGCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcgg                                                            >  1:4553138/1‑150 (MQ=255)
                                                                                           gCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcgg                                                            >  1:2703429/1‑150 (MQ=255)
                                                                                             tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAATAAACCGCCACTGTTGTCAgcggcg                                                          >  2:3966014/1‑150 (MQ=255)
                                                                                             tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg                                                          >  2:2166328/1‑150 (MQ=255)
                                                                                             tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg                                                          >  2:2106606/1‑150 (MQ=255)
                                                                                             tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg                                                          >  2:2106580/1‑150 (MQ=255)
                                                                                             tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg                                                          >  2:1803046/1‑150 (MQ=255)
                                                                                             tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg                                                          >  1:943454/1‑150 (MQ=255)
                                                                                             tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGATGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg                                                          >  2:1784941/1‑150 (MQ=255)
                                                                                             tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCAGccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggtataa                                                          <  1:3810025‑M2/150‑93 (MQ=255)
                                                                                               cGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgtt                                                                    <  2:3333783/138‑1 (MQ=255)
                                                                                               cGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgtt                                                                    >  1:3333783/1‑138 (MQ=255)
                                                                                                 aaGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGttt                                                      >  2:1939408/1‑150 (MQ=255)
                                                                                                 aaGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGttt                                                      >  2:168423/1‑150 (MQ=255)
                                                                                                  agcagcCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggtataaaccgc                                                     <  1:1744572‑M2/150‑98 (MQ=255)
                                                                                                   gcagcCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACTAGGCCGCTCCCACAAGAAAGCAGGCCGCCAGGGGGTTGACATCTCCGCAAAGGAAAGCGCCCCTGTTGTCGGCGGGGGtggtg                                                    >  2:2909787/1‑146 (MQ=25)
                                                                                                      gcCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTggg                                                                                                    <  2:481353/99‑1 (MQ=255)
                                                                                                      gcCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTggg                                                                                                    >  1:481353/1‑99 (MQ=255)
                                                                                                       ccccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagca                                                                                                                 <  2:3602979‑M2/85‑38 (MQ=255)
                                                                                                       ccccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagca                                                                                                                 >  1:3602979‑M2/1‑48 (MQ=255)
                                                                                                       ccccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTa                                                >  2:794941/1‑150 (MQ=255)
                                                                                                       ccccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTa                                                >  1:1969715/1‑150 (MQ=255)
                                                                                                       ccccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTa                                                >  1:5113506/1‑150 (MQ=255)
                                                                                                         ccGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTGGAGTATc                                              >  2:4633190/1‑150 (MQ=255)
                                                                                                          cGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCg                                             >  2:4761609/1‑150 (MQ=255)
                                                                                                           ggCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGc                                            >  2:2854110/1‑150 (MQ=255)
                                                                                                           ggCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGc                                            >  2:368344/1‑150 (MQ=255)
                                                                                                            gCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCa                                           >  1:4479045/1‑150 (MQ=255)
                                                                                                            gCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACAGTTGTCAGCGGCGGTTTTGAGTATCGCa                                           >  1:2068415/1‑150 (MQ=255)
                                                                                                            gCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCAGccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcga                                                                                            <  2:3075308‑M2/101‑59 (MQ=255)
                                                                                                            gCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCAGccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcga                                                                                            >  1:3075308‑M2/1‑43 (MQ=255)
                                                                                                            gCACTCTGTGCTGCGAAGCTGAAAGCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggtataaaccgcctcttattga                                           >  2:4462073‑M2/1‑43 (MQ=255)
                                                                                                                tctGTGCTGCGAAGCAGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggtataaaccgcctcttagtgaattg                                       <  2:3494842‑M2/150‑112 (MQ=255)
                                                                                                                  tgtgCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGccc                                     >  1:4356120/1‑150 (MQ=255)
                                                                                                                       tgcGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGc                                                                           <  1:2370353/107‑1 (MQ=255)
                                                                                                                       tgcGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGc                                                                           >  2:2370353/1‑107 (MQ=255)
                                                                                                                        gcGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGaa                               >  1:4438745/1‑150 (MQ=255)
                                                                                                                              cTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCt                                    >  2:3290704/1‑139 (MQ=255)
                                                                                                                              cTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCt                                    <  1:3290704/139‑1 (MQ=255)
                                                                                                                              cTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCAt                         >  2:4711004/1‑150 (MQ=255)
                                                                                                                               tGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATc                        >  2:4636210/1‑150 (MQ=255)
                                                                                                                                  gcgaaaggggCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTg                     >  2:532627/7‑150 (MQ=255)
                                                                                                                                  gcgaaaggggCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGGTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTg                     >  2:532631/7‑150 (MQ=255)
                                                                                                                                  aaTCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTg                     >  1:4348007/1‑150 (MQ=255)
                                                                                                                                     cccGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGa                                                   >  2:4475094/1‑117 (MQ=255)
                                                                                                                                     cccGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGa                                                   <  1:4475094/117‑1 (MQ=255)
                                                                                                                                       cGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGcc                >  2:1523927/1‑150 (MQ=255)
                                                                                                                                       cGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGcc                >  1:4006351/1‑150 (MQ=255)
                                                                                                                                              gCTTCGCGCCCCTTTCGCGACACAAGGCCGCTACCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGAGGCGGTTTTGAGTATCGCAATGCCCATAGAACGGCATCTTGAAGCCCggcggc         >  1:3485701/1‑150 (MQ=255)
                                                                                                                                               cTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCa        >  2:4656066/1‑150 (MQ=255)
                                                                                                                                                ttCGCAGCCCTTTCGCGACACAAGGCCCCACCCACAAAAAAGCAGACCGCAACTGGGTAGAAATCTCCGCAAAAAAAACCGCCACTGTAGTCAGCGCCGGTTATGAGTATCCCAATGCCCATAGAACGGCAACTTGAAGCCCGGCGGCAg       >  2:1173776/1‑150 (MQ=37)
                                                                                                                                                    cgcCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGc   >  1:874711/1‑150 (MQ=255)
                                                                                                                                                     gcCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCa  >  2:4437108/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTTTTGAGTATCGCAATGCCCTTAGAACGGCATCTTGAAGCCCGGCGGCAGCTGCA  >  NC_002947/4854329‑4854627

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: