Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,857,452 G→T 5.9% intergenic (‑129/+168) 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,4520GT 5.9% 207.8 / 2.0 73intergenic (‑129/+168)recR/ybaBrecombination protein RecR/nucleoid‑associated broad specificity regulator
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base T (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GATGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTT  >  NC_002947/4857303‑4857601
                                                                                                                                                      |                                                                                                                                                     
gATGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCg                                                                                                                                                       <  1:6523163/150‑1 (MQ=255)
 aTGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGa                                                                                                                                                      <  2:3307103/150‑1 (MQ=255)
  tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTCTACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGaa                                                                                                                                                     <  1:3798980/150‑1 (MQ=255)
  tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGaa                                                                                                                                                     <  2:1023636/150‑1 (MQ=255)
  tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGaa                                                                                                                                                     <  2:2114566/150‑1 (MQ=255)
  tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGaa                                                                                                                                                     <  1:3696431/150‑1 (MQ=255)
  tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGaa                                                                                                                                                     <  1:3771050/150‑1 (MQ=255)
  tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGaa                                                                                                                                                     <  1:5224222/150‑1 (MQ=255)
  tGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGaa                                                                                                                                                     <  1:5224163/150‑1 (MQ=255)
    agagGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAagc                                                                                                                                                   <  2:5553202/150‑1 (MQ=255)
          cTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGcccc                                                                                                                                             <  2:3570098/150‑1 (MQ=255)
           tGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCg                                                                                                                                            >  1:6705432/1‑150 (MQ=255)
             aaGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAagc                                                                                                                                                   >  2:5081467/1‑141 (MQ=255)
             aaGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAagc                                                                                                                                                   <  1:5081467/141‑1 (MQ=255)
             aaGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGccctttc                                                                                                                                          <  1:6052976/150‑5 (MQ=255)
              aGCTCTTGGGAAAGGCCTGACAATTCAATAAGTGGCGGTTTATACCCAGGGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCa                                                                                                                                         <  2:6705432/150‑1 (MQ=255)
              aGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCTTCGCAGccctttcg                                                                                                                                         <  2:5988285/150‑6 (MQ=255)
                 tCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCActc                                                                                                                                      >  2:5392496/1‑150 (MQ=255)
                             ccTGACAATACATTACGTGGGGTTTTATCCCCCGGGGGGGGGGTGGCGTAAAAACACCCCGCCCACCGCCGGCACGGTCTGCTTTTTTGGGGGCGGGGCCTTGTGTGGCG‑AATGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAg                                                                                                                          <  2:2254527/150‑1 (MQ=255)
                             ccTCACAATAAATTAAGCGGCGTTTTATACCCCGGGGGGGGGGTGGCGTAAAACCCCCGCGCCCAAAGCCTGAATGGTCTGTTTTTTTGTGGGAGCGGCGTTGTGTCGCG‑AACGGGCGGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAg                                                                                                                          <  2:4750795/150‑1 (MQ=255)
                              cttaaatttCAATCCGTGGGGGTTTATACCCTGGGGGGGGGGTGGCGTAAAATCACCTCGCCAAACGCCTGCATGGTCTGATTTTTTGTGGTAGGGGCCTTGTGTCGCT‑AACGGGCTGCGACGCCCCCCCGGCACTCTGTGCTGCGAAGc                                                                                                                         <  2:534124/143‑1 (MQ=255)
                              cTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCg                                                                                                                                                       <  1:1384485/120‑1 (MQ=255)
                              cTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCg                                                                                                                                                       >  2:1384485/1‑120 (MQ=255)
                                 aCAATCCAATTCGAGGCGGTTTCTCCCCCGGGGGGGGGGTGGCGTCAAAACACCGCGCCCAACGCCTCCACGGTCGGCTTTTTTGTGGGAGCGGCCTTGGGCCGCC‑ACCGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGa                                                                                                                      <  2:3031172/150‑1 (MQ=255)
                                  aaaTTCAATACGCGGCTGTTTAACCCCCGGGGGGGGGGTGGCGTCAAACCACCGCGCCCAACGCCTGCATGGTCTGCTTTTTGGTGGGAGCGGCCTGGTGGGGCG‑AACGGGCGGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGaa                                                                                                                     <  2:1694575/149‑1 (MQ=255)
                                   aaTTCATTAGGCGGGGTTTTATCCCCGGGGGGGGGGGTGGCTAAAAACCACCCCCCCAACCGCCTGCTTGTTGTGTTTTTTTGTGGGCGGGGCGTTGTGGGGCT‑ACTGGGGTGCGAAGCCCCCCCGGCACTCTGTGCTGCGAAGCTGaaa                                                                                                                    <  2:309543/150‑1 (MQ=255)
                                       cAATACGTGGGGTTTTATACCCCGGGGGGGGGTTGGCTAAAACTCCCCTCCCCCAACCCCTGCATGTTGTGCTTTTTTGTGGGTGCGGCCTTGGCTCGCT‑AACGGGCTGCTAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATccc                                                                                                                <  2:6299392/150‑1 (MQ=255)
                                          aaactggCGTTTTAAACCCCGGGGGGGGGGTGGCCACAAAACACCCCCCAAAAACACTCAAAGCTCTTATTTTTTGTGGTAGCGCACTTTTTTCGCC‑AATGGGCTGCAAACAACCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCggg                                                                                                             <  2:1368288/145‑1 (MQ=255)
                                            agagGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AACGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGc                                                                                                           <  2:6811992/150‑1 (MQ=255)
                                            agagGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAACTCCCGGGGc                                                                                                           >  1:2790088/1‑150 (MQ=255)
                                                      tatcCCCGGGGGGGGGGGTGGCGTCAACGCACCGCGCCCAACGCCGGCACGGTCTGCTTTTTTGTGGGCGCGGCGTTGGGGGGCG‑ACCGGGGTGCGACGCCGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCgcgc                                                                                                 <  2:5416718/150‑1 (MQ=255)
                                                       atcccccGCGGGGGCGGTGGCGCCAAAACCCCCCGCCAAACCCCTGCAAGTTCTGCTTTTTTGGGGGTGGGGCCTTGTGCGGCG‑AACGGGCGGCCTAGCCGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGcc                                                                                                <  2:2941609/147‑1 (MQ=255)
                                                       agaCCCCGGGGGGGGGGGGGCGTCAAACCACCCCCCCCAACGCCGCCCTGGTCGGGTTTTTGGGGGGTGCGGCCTTGTGTCGCG‑AAGGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGcc                                                                                                <  2:4044398/148‑1 (MQ=19)
                                                              ggggggggggggGCGTCAAACCACCCCCCCCAAAGCCCGATCGGTGTGCTTTTTGGGGGGGGCGGCCGTGTGACGCG‑ACCGGGCTGCGACGCACCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTcg                                                                                         <  2:199909/148‑1 (MQ=255)
                                                              gCGGGGGCGGTGGCTTCAAAGCCCCGCGCCAAACGCCTGAAAGGAGTGTTTTTTGGTGGGCGCGGCCTTGT‑TCCCGTAACGGGCTGCGACGACGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTcg                                                                                         <  2:6488833/150‑1 (MQ=255)
                                                                ggggggggTGGCGTAAACCCACCGCGCCCAACGCCGGCCGGGTGTGC‑TTTTTGGGGGAGCGGCCTTGTGTCGCT‑AACGGGCGGCCACGCCGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGa                                                                                      <  2:6278732/150‑1 (MQ=255)
                                                                 ggggggGTGGCGCCAAACCACCTCGCCCAACGCCTGCACGGTCTGCTTTTTTTTGGGAGCGGCCTGGTGTCGCG‑AACGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcga                                                                                      >  2:1161752‑M2/1‑139 (MQ=255)
                                                                 ggggCGGTGGGGACAAATCACCTCGCCCAACGCCTCCATTGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcga                                                                                      >  2:4218339‑M2/1‑139 (MQ=255)
                                                                 ggggCGGTGGCGTCAAATCACCTCGCCCAACGCCTGCATGGTCTGCTTTTTTGTGGGAGGGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcga                                                                                      >  2:4218292‑M2/1‑139 (MQ=255)
                                                                 ggggCGGGGGCTTCAAAGCACCGCGCCCAACGCCTGCAAGGGCGGCTTTTTTGGGGGCGCGGCGTTGGGGCGCG‑ACGGGGCTGCGACGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGa                                                                                      <  2:5957906/150‑1 (MQ=255)
                                                                      gggggCGTCAAAACCCCCCCCCCAACGCCTGAATGGTGTTTTTTTTTGTGGGTGCGGCCTTGTGTCGCC‑AAGGGGCTGCCAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACaa                                                                                 <  2:4011593/147‑1 (MQ=21)
                                                                                                                              gCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAc                         >  2:4541672/1‑150 (MQ=255)
                                                                                                                               ccTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAAcc                        >  2:5710648/1‑150 (MQ=255)
                                                                                                                                 ttGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGc                      >  1:3865354/1‑150 (MQ=255)
                                                                                                                                    tgtCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCccctttcgcgacacaaggccgctcccacaaaaaagcagaccttgcaggctttgggcgaggtgatttgacgccaccgcc                   <  2:4800796‑M2/150‑79 (MQ=255)
                                                                                                                                    tgtCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCAc                   >  2:2471043/1‑150 (MQ=255)
                                                                                                                                     gtCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACt                  >  1:6213179/1‑150 (MQ=255)
                                                                                                                                      tCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCAcgg                 >  2:5558060/1‑148 (MQ=255)
                                                                                                                                      tCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtg                 >  1:702088/1‑150 (MQ=255)
                                                                                                                                      tCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtg                 >  1:783684/1‑150 (MQ=255)
                                                                                                                                        gcg‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgtt               >  2:4384402/1‑150 (MQ=255)
                                                                                                                                         cg‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:1089578/1‑150 (MQ=255)
                                                                                                                                         cg‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:4826858/1‑150 (MQ=255)
                                                                                                                                         cg‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:5940800/1‑150 (MQ=255)
                                                                                                                                         cg‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:5940958/1‑150 (MQ=255)
                                                                                                                                         cg‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttg              >  1:6467298/1‑150 (MQ=255)
                                                                                                                                          g‑AAAGGGCTGTGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttgt             >  1:950831/1‑150 (MQ=255)
                                                                                                                                          g‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttgt             >  2:6690673/1‑150 (MQ=255)
                                                                                                                                          g‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACtgttgt             >  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)
                                                                                                                                                      |                                                                                                                                                     
GATGAGAGGGCTGAAGCTCATGGGAAAGGCCTGACAATTCAATAAGAGGCGGTTTATACCCAGCGGGGGCGGTGGCGTCAAATCACCTCGCCCAAAGCCTGCAAGGTCTGCTTTTTTGTGGGAGCGGCCTTGTGTCGCG‑AAAGGGCTGCGAAGCAGCCCCGGCACTCTGTGCTGCGAAGCTGAAATCCCGGGGCCGCTTCGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAAAAAAGCAGACCGCAACTGGGTTGAAATCTCCGCAAAAAAAACCGCCACTGTTGTCAGCGGCGGTT  >  NC_002947/4857303‑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: ‑