Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,171,8820CG44.2% 58.4 / 22.9 43intergenic (+126/‑175)rrfB/murB5S ribosomal RNA of rrnB operon/UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (19/5);  new base G (0/19);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 5.37e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTA  >  NC_000913/4171732‑4172023
                                                                                                                                                      |                                                                                                                                             
gCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                                                                                                               <  1:73735/151‑1 (MQ=31)
  gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCgg                                                                                                                                             >  2:832217/1‑151 (MQ=31)
        aGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGt                                                                                                                                        <  1:157915/150‑1 (MQ=31)
         gggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTgg                                                                                                                                      >  2:676431/1‑151 (MQ=31)
              cTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGc                                                                                                                                 >  2:855235/1‑151 (MQ=31)
                 ccAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                              <  2:416976/151‑1 (MQ=31)
                 ccAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                              <  1:535253/151‑1 (MQ=31)
                      cATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcc                                                                                                                         >  1:770974/1‑151 (MQ=31)
                       aTCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgccc                                                                                                                        <  2:1017337/151‑1 (MQ=31)
                          aaaTAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgcc                                                                                                                     <  2:426302/151‑1 (MQ=31)
                             tAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAt                                                                                                                   >  1:451047/1‑150 (MQ=31)
                             aaaaaAGAAAGCCTCAGCCAAAGGCCTGGGCTTTTCGTTTTATCTGTTGTTTGTGGGTGACCGCTCTCCTGAGTAGGAAAACCCCGCGGGGAGCGAATTTAAACGTTGGGAAGCAACGGCGGGGGGGGTGGGGGGAAGGACGCCCGCCATa                                                                                                                  <  1:989811/150‑1 (MQ=2)
                                  cGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt                                                                                                              >  1:917540/1‑150 (MQ=31)
                                   gAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTg                                                                                                             >  1:665003/1‑150 (MQ=31)
                                    aagggCTCAGTCAAAAGACTGGGCTTTTCTTTTTATCTTTTTTTTGTGGGTGACCCCTCTCCTGAGTAGGACAAACCCGCGGGGGGCGGATTTGAACGTTGCAAAGAAACGGCCGGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt                                                                                                              <  2:24976/145‑1 (MQ=11)
                                    aaaGGCACAGGGGAAAGACGGGGCCTTTCTTTTTTTTGGTTGTTTGTCGGTGAACGCTCTCCTGAGAAGAAAAAACCCGCGGGGAGCGGTTTTAAACGTTGAAAAGCAACGCCCGGGGGGGGGGGGGGCGGGACGCCCGCCATAAACt                                                                                                              <  1:995177/148‑1 (MQ=2)
                                    aaaGCCTCAGTGAAAAGTCTGGGCCTTTCGTTTTACCTTTTGTTTGTGGGTGAACCCTCTCGGGAGAAGAACAAACCCGCGGGGAGGGTATTTGAACGTTGCAAAGCAACGGCGGGGAGGGGGGGGGGAAGGACGCCCGCCATAAACTg                                                                                                             <  1:106901/149‑1 (MQ=2)
                                      aGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGcc                                                                                                           >  1:291913/1‑149 (MQ=31)
                                                 aaaGACGGGGCTTTTTTTTTTTTTTGTTTTTTGTCGGGAAACCCTCTCTGGAGTAGGAAAACCCCGCGGGGAGGGGTTTTAAACGTTGCGAAGCAACGCCCGGGGGGGGGGCGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAAt                                                                                                <  2:765196/149‑1 (MQ=2)
                                                                                       gAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGttt                                                        >  1:847241/1‑151 (MQ=32)
                                                                                        accGCTCTCCTGGGTAGGCAAACCCCGCGGGGGGCGGTTTTGAACTTTGAAAAGAAAGGGCCGGGGGGGGGGGGGGGAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGttt                                                        <  2:503861/148‑1 (MQ=2)
                                                                                          cGCTCTCCTGAGTAGGAAAAACCCGGCGGGAGCGGTTTTAAAGTTTGCAAAGAAAGGGCCGGGGGGGTGGCGGTAAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGttt                                                        <  2:1011306/148‑1 (MQ=11)
                                                                                             tctcgggAGTGGGACAACCCCGCGGGGGGGGATTTTAAACTTTGCGAAGCACGGGCGGGGGGGGGGGCGGGCGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTc                                                       <  1:3867/140‑1 (MQ=2)
                                                                                               tgaggAGTAGGAAAAACCCGCGGGGGGGGGTTTGGAATTTTGCGAAGAAGCGGGCGGGGGGGGGGCGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTAc                                                    <  2:745633/143‑1 (MQ=11)
                                                                                               tcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAc                                                >  2:997036/1‑151 (MQ=38)
                                                                                               tcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGACAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAc                                                >  1:605615/1‑151 (MQ=38)
                                                                                                ccTGAGTAGAACAAACCCGCGGGGAGCGTATTTAAACGTTGGGAAGCAACGGCCGGGGGGGGGGGGGGAAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAct                                               <  2:681956/151‑1 (MQ=14)
                                                                                                     gagggACAAACCCGCCGGGCGGGTTTTTAAACTTTGCAAACAAAGGGCCGGGGGGGTGGGGGGCGGGCCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAct                                               <  2:638970/143‑1 (MQ=14)
                                                                                                        gaaCAAATCCGCGGGGAGGGGATTAGAGCGTTGCAAAGCAACGGCCGGGGGGGGGGCGGGCAGGACCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTg                                        <  2:347718/148‑1 (MQ=14)
                                                                                                          aaaaaCCCGCGGGGAGGGAATTTGAACGTTGCGAAGAACGGGCCGGGGGGGGGGGGGGCGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAAct                                               <  2:129424/139‑1 (MQ=14)
                                                                                                             aaTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTgtcgtc                                   >  2:429644/1‑150 (MQ=255)
                                                                                                               tCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCat                                 >  1:546331/1‑150 (MQ=255)
                                                                                                                ccgcggGGAGCGGATTTGAACGTTGCGAAGAAACGGCCGGGGGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCat                                 <  2:84540/149‑1 (MQ=255)
                                                                                                                         gCGGATTTGAACGTTGCGAAGCAACGGCGGGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGcc                      <  2:681217/151‑1 (MQ=255)
                                                                                                                             aTTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcc                  >  1:848061/1‑151 (MQ=255)
                                                                                                                             aTTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGAAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcc                  >  1:787342/1‑151 (MQ=255)
                                                                                                                              tttGACGGTTGGGAAGCAGGGGCCGGGGGGGGGGGGGGAGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATcc                  <  1:826857/150‑1 (MQ=18)
                                                                                                                              tttGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCCGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTATTCCTGTCGTCATATCTACAAGCCATccc                 >  2:330582/1‑151 (MQ=255)
                                                                                                                                tGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATccccc               >  2:269986/1‑151 (MQ=255)
                                                                                                                                    cGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcac            >  1:176369/1‑150 (MQ=255)
                                                                                                                                      ttgaaaagaaccGGCCGGGGGGGGGGGGGGGGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCcaca           <  1:210789/138‑1 (MQ=21)
                                                                                                                                           aagaaccGGCCGGGGGGGGGGCGGGCGGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATa       <  2:863441/142‑1 (MQ=21)
                                                                                                                                           aaGCAACGGCCGGGGGGGGGGCGGGAGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCca             <  1:308420/142‑1 (MQ=21)
                                                                                                                                              cAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTa  >  1:240945/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
GCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTA  >  NC_000913/4171732‑4172023

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