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

GAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCG  >  NC_000913/4040614‑4040912
                                                                                                                                                     |                                                                                                                                                      
gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCgg                                                                                                                                                      >  2:15338/1‑151 (MQ=31)
   agTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAgg                                                                                                                                                   <  2:657130/151‑1 (MQ=31)
    gTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAggg                                                                                                                                                  >  1:907741/1‑151 (MQ=31)
         gAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCg                                                                                                                                             <  2:322934/151‑1 (MQ=31)
           aCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCggg                                                                                                                                           >  1:144857/1‑151 (MQ=31)
                cAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAgg                                                                                                                                       >  1:626107/1‑150 (MQ=31)
                    cATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGAcgcc                                                                                                                                  >  1:224571/1‑151 (MQ=31)
                         aaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCa                                                                                                                             <  1:384604/151‑1 (MQ=31)
                           tAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCAt                                                                                                                            >  1:447051/1‑150 (MQ=31)
                                  aaaGGCCCAGTCGGAAGACTGGCCTTTTCTTTTTATTTGTTTTTTGTCGGTGACCCCTCTCCTGAGAAGAAAAAACCCGCGGGGAGCGGATTTAAACGTTGCGAA‑GCAACGGCCGGGAGGGTGGCGGGCAGGACGCCCGCCATAAACt                                                                                                                       <  2:331206/148‑1 (MQ=11)
                                      gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGc                                                                                                                >  2:582168/1‑151 (MQ=31)
                                       cTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGc                                                                                                                >  2:719406/1‑150 (MQ=31)
                                            tCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaa                                                                                                           >  2:904332/1‑150 (MQ=31)
                                                                                        cGCTCTCCGGAGTGGGACAACCCCGCGGGGGGGGGTTTTGAACGTTGCGAA‑GAACCGGCCGGGAGGGGGGGGGGAAGGCCCCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCAtt                                                                  <  2:850422/147‑1 (MQ=11)
                                                                                        cGCTCCCCGGAGTAGGAAAAACCCGCGGGGGGCGGTTTTGAACTTTGCGAA‑GCACCGGCCGGGGGGGGGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgcg                                                              <  1:274108/151‑1 (MQ=11)
                                                                                         gCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgcgc                                                             >  2:66272/1‑151 (MQ=35)
                                                                                           tctcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGgcgc                                                             >  2:841689/1‑149 (MQ=35)
                                                                                            ctcGTGAGTAGAAAAAACCCGCGGGGTGCGGTTTTGAGCG‑TGCGAAGGCAGCGCCCGGGGGGGTGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAg                                                           <  2:54656/150‑1 (MQ=11)
                                                                                             tcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGa                                                          >  2:108347/1‑150 (MQ=38)
                                                                                             tcCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGa                                                          >  1:855859/1‑150 (MQ=38)
                                                                                              ccTGAGTAGGAAAAACCGGGCGGGGGGGGATTTGAACTTTGCGAA‑GCACCGGCCCGGAGGGTGGCGGGCGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaa                                                         <  1:378816/150‑1 (MQ=14)
                                                                                              ccTGAGAAGAAAAACCCCGCGGGGGGCGGATTGGAACGTGGCGAA‑GCAACGCCCGGGGGGGTGGGGGGCGGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaa                                                        <  1:692110/151‑1 (MQ=14)
                                                                                                tGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGaaaa                                                       >  1:818528/1‑150 (MQ=40)
                                                                                                          aaccccgcggGGGGGGGTTTTGAACTTTGAAAA‑GCACGGGCCGGGGGGGGGGGGGGCGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATg                                                  <  2:308312/141‑1 (MQ=14)
                                                                                                          aaaCCCGCGGGGAGCGGATTTGAACGTTGCGAA‑GCAGCGGCCGGGGGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGa                                             <  2:750325/150‑1 (MQ=17)
                                                                                                           aacccgcggGGAGGGGATTTAACCGTTGGGAA‑GAACCGGCCGGGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGAt                                            <  1:207655/147‑1 (MQ=17)
                                                                                                                gcgggggGCGGATTTGAACGTTGCGAA‑GCAACGGCCGGGGGGGTGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGAc                                       <  1:787858/147‑1 (MQ=17)
                                                                                                                  gggggggggttttGAACGTGGCAAA‑GAACCGGCGGGGGGGGGGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGAc                                       <  2:396837/138‑1 (MQ=17)
                                                                                                                   gggggggtttttGAACGTTGCAAA‑GCAGGGGCCGGGGGGGTGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACg                                      <  1:415008/139‑1 (MQ=17)
                                                                                                                   gggAGCGGATTGGACCTTTGCGAA‑GCACCGGCCGGGGGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGAc                                       <  2:274769/147‑1 (MQ=17)
                                                                                                                    gggggggATTTGAACGTGGCGAA‑GCACCGGCCGGGGGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTgg                                  <  2:777008/146‑1 (MQ=18)
                                                                                                                          gATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGt                             >  2:571749/1‑150 (MQ=255)
                                                                                                                             ttaaaCTTTGCGAA‑GCACGGGCCGGGGGGGGGGGGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCtt                          <  1:374451/147‑1 (MQ=18)
                                                                                                                                  cGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATcc                      >  2:221472/1‑149 (MQ=255)
                                                                                                                                                 ggCCGGGGGGGGGGGGGGAGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTAct              <  1:50394/143‑1 (MQ=18)
                                                                                                                                                 ggCCGGGGGGGGGGGGGGAGGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCt        <  1:903946/149‑1 (MQ=21)
                                                                                                                                                 ggCCGGGAGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCtt           <  1:418008/146‑1 (MQ=255)
                                                                                                                                                     cGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCg  >  1:306544/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAA‑GCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATGCCGGATGCGACGCTGGCGCGTCTTATCCAGCCTACTCTTGCTTATCCG  >  NC_000913/4040614‑4040912

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