Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,857,454 A→C 20.9% intergenic (‑131/+166) 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,857,4540AC20.9% 141.2 / 8.5 64intergenic (‑131/+166)recR/ybaBrecombination protein RecR/nucleoid‑associated broad specificity regulator
Reads supporting (aligned to +/- strand):  ref base A (33/16);  new base C (3/12);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 2.22e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTT  >  NC_002947/4857305‑4857601
                                                                                                                                                      |                                                                                                                                                   
tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTCTACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGaa                                                                                                                                                     <  1:3798980/150‑1 (MQ=255)
tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGaa                                                                                                                                                     <  2:2114566/150‑1 (MQ=255)
tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGaa                                                                                                                                                     <  1:3696431/150‑1 (MQ=255)
tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGaa                                                                                                                                                     <  1:3771050/150‑1 (MQ=255)
tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGaa                                                                                                                                                     <  2:1023636/150‑1 (MQ=255)
tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGaa                                                                                                                                                     <  1:5224222/150‑1 (MQ=255)
tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGaa                                                                                                                                                     <  1:5224163/150‑1 (MQ=255)
  agagGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAagc                                                                                                                                                   <  2:5553202/150‑1 (MQ=255)
        cTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGcccc                                                                                                                                             <  2:3570098/150‑1 (MQ=255)
         tGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCg                                                                                                                                            >  1:6705432/1‑150 (MQ=255)
           aaGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAagc                                                                                                                                                   <  1:5081467/141‑1 (MQ=255)
           aaGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAagc                                                                                                                                                   >  2:5081467/1‑141 (MQ=255)
           aaGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGccctttc                                                                                                                                          <  1:6052976/150‑5 (MQ=255)
            aGCTCTTGGGAAAGGCCTGACAATTCAATAAGTGGCGGTTTATACCCAGGGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCa                                                                                                                                         <  2:6705432/150‑1 (MQ=255)
            aGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCTTCGCAGccctttcg                                                                                                                                         <  2:5988285/150‑6 (MQ=255)
               tCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCActc                                                                                                                                      >  2:5392496/1‑150 (MQ=255)
                           ccTGACAATACATTACGTGGGGTTTTATCCCCCGGGGGGGGGGTGGCGTAAAAACACCCCGCCCACCGCCGGCACGGTCTGCTTTTTTGGGGGCGGGGCCTTGTGTGGCGAATGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAg                                                                                                                          <  2:2254527/150‑1 (MQ=255)
                           ccTCACAATAAATTAAGCGGCGTTTTATACCCCGGGGGGGGGGTGGCGTAAAACCCCCGCGCCCAAAGCCTGAATGGTCTGTTTTTTTGTGGGAGCGGCGTTGTGTCGCGAACGGGCGGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAg                                                                                                                          <  2:4750795/150‑1 (MQ=255)
                            cttaaatttCAATCCGTGGGGGTTTATACCCTGGGGGGGGGGTGGCGTAAAATCACCTCGCCAAACGCCTGCATGGTCTGATTTTTTGTGGTAGGGGCCTTGTGTCGCTAACGGGCTGCGACGCCCCCCCGGCACTCTGTGCTGCGAAGc                                                                                                                         <  2:534124/143‑1 (MQ=255)
                               aCAATCCAATTCGAGGCGGTTTCTCCCCCGGGGGGGGGGTGGCGTCAAAACACCGCGCCCAACGCCTCCACGGTCGGCTTTTTTGTGGGAGCGGCCTTGGGCCGCCACCGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGa                                                                                                                      <  2:3031172/150‑1 (MQ=255)
                                aaaTTCAATACGCGGCTGTTTAACCCCCGGGGGGGGGGTGGCGTCAAACCACCGCGCCCAACGCCTGCATGGTCTGCTTTTTGGTGGGAGCGGCCTGGTGGGGCGAACGGGCGGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGaa                                                                                                                     <  2:1694575/149‑1 (MQ=255)
                                 aaTTCATTAGGCGGGGTTTTATCCCCGGGGGGGGGGGTGGCTAAAAACCACCCCCCCAACCGCCTGCTTGTTGTGTTTTTTTGTGGGCGGGGCGTTGTGGGGCTACTGGGGTGCGAAGCCCCCCCGGCACTCTGTGCTGCGAAGCTGaaa                                                                                                                    <  2:309543/150‑1 (MQ=255)
                                     cAATACGTGGGGTTTTATACCCCGGGGGGGGGTTGGCTAAAACTCCCCTCCCCCAACCCCTGCATGTTGTGCTTTTTTGTGGGTGCGGCCTTGGCTCGCTAACGGGCTGCTAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATccc                                                                                                                <  2:6299392/150‑1 (MQ=255)
                                        aaactggCGTTTTAAACCCCGGGGGGGGGGTGGCCACAAAACACCCCCCAAAAACACTCAAAGCTCTTATTTTTTGTGGTAGCGCACTTTTTTCGCCAATGGGCTGCAAACAACCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCggg                                                                                                             <  2:1368288/145‑1 (MQ=255)
                                          agagGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAACGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGc                                                                                                           <  2:6811992/150‑1 (MQ=255)
                                          agagGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAACTCCCGGGGc                                                                                                           >  1:2790088/1‑150 (MQ=255)
                                                    tatcCCCGGGGGGGGGGGTGGCGTCAACGCACCGCGCCCAACGCCGGCACGGTCTGCTTTTTTGTGGGCGCGGCGTTGGGGGGCGACCGGGGTGCGACGCCGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCgcgc                                                                                                 <  2:5416718/150‑1 (MQ=255)
                                                     atcccccGCGGGGGCGGTGGCGCCAAAACCCCCCGCCAAACCCCTGCAAGTTCTGCTTTTTTGGGGGTGGGGCCTTGTGCGGCGAACGGGCGGCCTAGCCGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGcc                                                                                                <  2:2941609/147‑1 (MQ=255)
                                                     agaCCCCGGGGGGGGGGGGGCGTCAAACCACCCCCCCCAACGCCGCCCTGGTCGGGTTTTTGGGGGGTGCGGCCTTGTGTCGCGAAGGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGcc                                                                                                <  2:4044398/148‑1 (MQ=19)
                                                            ggggggggggggGCGTCAAACCACCCCCCCCAAAGCCCGATCGGTGTGCTTTTTGGGGGGGGCGGCCGTGTGACGCGACCGGGCTGCGACGCACCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTcg                                                                                         <  2:199909/148‑1 (MQ=255)
                                                            gCGGGGGCGGTGGCTTCAAAGCCCCGCGCCAAACGCCTGAAAGGAGTGTTTTTTGGTGGGCGCGGCCTTGTTCCCGTAACGGGCTGCGACGACGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTcg                                                                                         <  2:6488833/150‑1 (MQ=255)
                                                              ggggggggTGGCGTAAACCCACCGCGCCCAACGCCGGCCGGGTGTGCTTTTTGGGGGAGCGGCCTTGTGTCGCTAACGGGCGGCCACGCCGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGa                                                                                      <  2:6278732/150‑1 (MQ=255)
                                                               ggggggGTGGCGCCAAACCACCTCGCCCAACGCCTGCACGGTCTGCTTTTTTTTGGGAGCGGCCTGGTGTCGCGAACGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcga                                                                                      >  2:1161752‑M2/1‑139 (MQ=255)
                                                               ggggCGGTGGGGACAAATCACCTCGCCCAACGCCTCCATTGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcga                                                                                      >  2:4218339‑M2/1‑139 (MQ=255)
                                                               ggggCGGTGGCGTCAAATCACCTCGCCCAACGCCTGCATGGTCTGCTTTTTTGTGGGAGGGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcga                                                                                      >  2:4218292‑M2/1‑139 (MQ=255)
                                                               ggggCGGGGGCTTCAAAGCACCGCGCCCAACGCCTGCAAGGGCGGCTTTTTTGGGGGCGCGGCGTTGGGGCGCGACGGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGa                                                                                      <  2:5957906/150‑1 (MQ=255)
                                                                    gggggCGTCAAAACCCCCCCCCCAACGCCTGAATGGTGTTTTTTTTTGTGGGTGCGGCCTTGTGTCGCCAAGGGGCTGCCAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACaa                                                                                 <  2:4011593/147‑1 (MQ=21)
                                                                                                                            gCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAc                         >  2:4541672/1‑150 (MQ=255)
                                                                                                                             ccTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAcc                        >  2:5710648/1‑150 (MQ=255)
                                                                                                                               ttGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGc                      >  1:3865354/1‑150 (MQ=255)
                                                                                                                                  tgtCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcc                   <  2:4800796‑M2/150‑79 (MQ=255)
                                                                                                                                  tgtCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCAc                   >  2:2471043/1‑150 (MQ=255)
                                                                                                                                   gtCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACt                  >  1:6213179/1‑150 (MQ=255)
                                                                                                                                    tCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCAcgg                 >  2:5558060/1‑148 (MQ=255)
                                                                                                                                    tCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtg                 >  1:702088/1‑150 (MQ=255)
                                                                                                                                    tCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtg                 >  1:783684/1‑150 (MQ=255)
                                                                                                                                      gcgAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgtt               >  2:4384402/1‑150 (MQ=255)
                                                                                                                                       cgAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:1089578/1‑150 (MQ=255)
                                                                                                                                       cgAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:4826858/1‑150 (MQ=255)
                                                                                                                                       cgAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:5940800/1‑150 (MQ=255)
                                                                                                                                       cgAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:5940958/1‑150 (MQ=255)
                                                                                                                                       cgAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:6467298/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGTGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttgt             >  1:950831/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttgt             >  2:6690673/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttgt             >  1:4487314/1‑150 (MQ=255)
                                                                                                                                          ccggggCTGCTTCGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTGCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTc            >  2:4538635/4‑150 (MQ=38)
                                                                                                                                           aaGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGGTGTcc           >  2:4008475/1‑149 (MQ=255)
                                                                                                                                            ggggCTGCTTCGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAg          >  2:5642135/2‑150 (MQ=255)
                                                                                                                                             gggCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggt         <  1:2940962‑M2/150‑89 (MQ=255)
                                                                                                                                              ggatGCGAAGCAGCACCGGCACTCTGTGCAGCGAAGCTGAACTCCCGGGGCCGCTTCGCGCACCTTTCGCGACCCAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAACGCCACGGTTGGCAccg        >  2:6725635/4‑150 (MQ=37)
                                                                                                                                              ggCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggta        <  1:6746546‑M2/150‑90 (MQ=255)
                                                                                                                                              ggCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggta        <  2:4599114‑M2/150‑90 (MQ=255)
                                                                                                                                                cTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggc      >  1:5294271/1‑150 (MQ=255)
                                                                                                                                                 tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggtataa     <  2:1572842‑M2/150‑93 (MQ=255)
                                                                                                                                                 tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcccccgctgggtataa     <  2:2458098‑M2/150‑93 (MQ=255)
                                                                                                                                                 tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg     >  2:1664759/1‑150 (MQ=255)
                                                                                                                                                 tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg     >  2:1910325/1‑150 (MQ=255)
                                                                                                                                                 tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg     >  2:1910476/1‑150 (MQ=255)
                                                                                                                                                 tGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggcg     >  2:5885367/1‑150 (MQ=255)
                                                                                                                                                  gCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggc      >  1:114724/1‑148 (MQ=255)
                                                                                                                                                  gCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAgcggc      <  2:114724/148‑1 (MQ=255)
                                                                                                                                                   cGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGt   >  2:3515154/1‑150 (MQ=255)
                                                                                                                                                    gAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGcggcgggt  >  2:6992053/1‑148 (MQ=255)
                                                                                                                                                    gAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACGGGTGTCAGCGGCGGtt  >  2:4709080/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTT  >  NC_002947/4857305‑4857601

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: