Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,040,7620CG42.2% 78.5 / 21.4 45intergenic (+126/+144)rrfA/mobB5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base C (24/2);  new base G (0/19);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 8.61e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTAT  >  NC_000913/4040612‑4040909
                                                                                                                                                       |                                                                                                                                                   
gCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGccc                                                                                                                                                     <  1:363078/151‑1 (MQ=31)
     agTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg                                                                                                                                                >  2:506734/1‑151 (MQ=31)
            aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCgg                                                                                                                                         >  2:604009/1‑151 (MQ=31)
                   aGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAc                                                                                                                                  <  2:1170215/151‑1 (MQ=31)
                    ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcg                                                                                                                                 <  2:501998/151‑1 (MQ=31)
                     gCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgc                                                                                                                                >  2:172183/1‑151 (MQ=31)
                      cATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcc                                                                                                                               >  2:151330/1‑151 (MQ=31)
                         cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgc                                                                                                                            >  1:347228/1‑151 (MQ=31)
                                           cAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc                                                                                                          >  2:595941/1‑151 (MQ=31)
                                            aGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc                                                                                                          >  2:866560/1‑150 (MQ=31)
                                            aGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc                                                                                                          >  2:1000216/1‑150 (MQ=31)
                                               cGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAt                                                                                                      >  2:741348/1‑151 (MQ=31)
                                                ggAAGCCTGGCCT‑TTTCGTTTTATCTGTTTTTGGTGGGTAACCGCTCTCCGGAGTAGGAAAAACCGGCGGGGGGGGGTTTTAAAGTTTGCGAAGCAAGGGCCGGGGGGGGGGGGGGCGGGCCCCCCGCCATAAACTGCCAGGCAt                                                                                                           <  1:457445/145‑1 (MQ=2)
                                                 gaaaggcgggCCTTTTCTTTTTATCTTTTTTTTGTGGTTGACCCCCCTCGGGAGAAGAAAAAACCCGCCGGGCGCGTATTTGAAGGTGGCAAAGCAACGGCCGGGGGGGGGGGGGGCAGAACGCCCGCCATAAACTGCCAGGCATCAAAt                                                                                                      <  1:1020057/143‑1 (MQ=2)
                                                 aaagagggggCC‑TTTCTTTTTTTTTGTTTTTTGTGGGTGAACGCTCTCCTGAGTGGGACAAACCCGCGGGGAGCGGATTAAAACTTGGCAAAGCAACGGCCGGGAGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTa                                                                                                    <  1:608168/144‑1 (MQ=11)
                                                                                     ggtaacccccctcCTGAGAAGAAAAAACCCGCGGGGGGGGAATTTGAGGGTTGGGAAGCAACGGCCGGGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCa                                                                 <  2:471598/141‑1 (MQ=11)
                                                                                      gTGACCCCTCCCCTGAGAAGAAAAACTCCGCGGGGGGGGGATTGGAACGTGGCGAAGAACCGCCGGGGGGGGGGGGGGGGAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCAt                                                                <  1:1068823/151‑1 (MQ=2)
                                                                                         accGCTCTCTGGAGTAGAACAAACCCGCGGGGGGGGGATTTGAACGTTGAGAAGCACCGGCCGGGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTg                                                              <  2:936321/148‑1 (MQ=11)
                                                                                         accGCTCTCCGGAGTAGGAAAAACCCGCGGGGGGGGGTTTTGAATGTTGCAAAGCACCGGCCGGGGGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCAtt                                                               <  2:1099365/147‑1 (MQ=11)
                                                                                           cGCTCTCCGGGGAAGAAAAACCCCGCGGGGAGCGGATTGAAACGTGGCAAAGCACCGGCCGGGGGGGGGGGGGGCGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTgg                                                             <  1:106855/149‑1 (MQ=11)
                                                                                            gCTCTCGTGAGAAGGAAAAACCCGCGGGGGGGGTATTTGATCGTGGGGAAGAACCGGCCGGGGGGGGGGCGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTg                                                              <  1:563059/147‑1 (MQ=11)
                                                                                                   tGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaa                                                    >  2:327257/1‑150 (MQ=40)
                                                                                                     aGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaaaa                                                  >  2:712860/1‑150 (MQ=255)
                                                                                                     aGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaaaa                                                  >  2:587147/1‑150 (MQ=255)
                                                                                                       tAGGAAAAACCCGCCGGGAGCGGATATAAGGGTCCGGAAGCAGGGGCGGGGGGGGTGGCGGGAGGGACCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAAt                                                <  1:198741/150‑1 (MQ=14)
                                                                                                        aGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGc                                              >  1:178215/1‑151 (MQ=255)
                                                                                                           aaaaccccgcggGGAGCGATTTTGAAGGTTGAGAAGCAGCGGCCGGGGGGGGGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGa                                                       <  2:553946/133‑1 (MQ=14)
                                                                                                           aaaaaCCCGCGGGGGGGGGATTTGACCTTTGCGAAGCACCGGCCGGGGGGGGGGGGGGAGGGACCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCgg                                           <  2:903089/149‑1 (MQ=14)
                                                                                                           aCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCgg                                           >  1:747347/1‑151 (MQ=255)
                                                                                                             aaaTCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGCGGCGGGCAGGACGCCCGCCATAAACTACCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGAt                                         >  1:290443/1‑151 (MQ=255)
                                                                                                                 ccgcggGGGGGGGTTTTGAACTTTGCGAAGCACGGGCCGGGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGc                                       <  1:797724/149‑1 (MQ=17)
                                                                                                                 ccgccgGGGGGGGATTTGAAGTTTGCAAAGCACGGGCCGGGGGGGGGGGGGGAAGGCCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGc                                       <  2:1250769/149‑1 (MQ=17)
                                                                                                                 ccgccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCg                                      >  1:115868/1‑150 (MQ=255)
                                                                                                                    ccgGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGc                                  >  2:952404/1‑151 (MQ=255)
                                                                                                                      gggAGCGTTTTTGAGGGTGGCGAAGAAGCGGCCGGGGGGGGGGGGGGCGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCt                                 <  1:437034/150‑1 (MQ=17)
                                                                                                                       ggAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTgg                               >  1:469332/1‑151 (MQ=255)
                                                                                                                        ggggggATTGAAAGGTTGGGAAGCACCGGCCGGGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGgc                              <  1:1068008/147‑1 (MQ=18)
                                                                                                                              aTTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCt                        >  2:875979/1‑151 (MQ=255)
                                                                                                                                 tGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTAt                     >  2:582783/1‑151 (MQ=255)
                                                                                                                                 tGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTAt                     >  2:489391/1‑151 (MQ=255)
                                                                                                                                       ttGCGAAGCACCGCCCGGGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAg                 <  2:625229/149‑1 (MQ=21)
                                                                                                                                       ttGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGcc               >  1:63426/1‑151 (MQ=255)
                                                                                                                                       ttGAGAAGCACGGGCCGGGGGGGGGGCGGGGAGAACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGc                <  1:2696/150‑1 (MQ=18)
                                                                                                                                             aGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTActct         >  2:999472/1‑151 (MQ=255)
                                                                                                                                               caccGGCCGGGGGGGTGGGGGGCAGGGGGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCtt        <  2:291048/147‑1 (MQ=21)
                                                                                                                                                    gCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTAt  >  1:846387/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
GCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCC‑TTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTAT  >  NC_000913/4040612‑4040909

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