Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_pgi380,0130G.20.0% 73.9 / 17.4 30intergenic (‑132/+56)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (24/0);  new base . (6/0);  total (30/0)
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.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAG  >  NC_000913_3_hsa_pgi/379874‑380155
                                                                                                                                            |                                                                                                                                              
cTGGGCATTTCGCACCACATCATCTTCATTCAATTAAAATTAGGGGGGTATTCTATATGTCAATTCAAACCCCCCTATAGTATATAACCCTGGAGAATAAAATTTATCCGGTGAAAGTGGCCGGAAAACAAACAGGAAA‑gggggggggg                                                                                                                                       <  2:407526/149‑1 (MQ=255)
 tGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGTGGGGGGGc                                                                                                                                      >  2:1330268/1‑149 (MQ=255)
                           atgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGATAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGa                                                                                                            >  2:1403311/1‑149 (MQ=255)
                            tgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGATTGTGGTCGGAAAACAAAGAGGATA‑GGGGGGGGGGGTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                           >  1:1272115/1‑149 (MQ=255)
                            tgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATA‑CCCCCTATAGTATATAACAGTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGTTa                                                                                                          >  1:756960/1‑149 (MQ=255)
                                        taGGGGGGGATTCTATATGTCAATGCATACCCCCCTATAGGATATAACACGGGAGCATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGggaggccg                                                                                               >  1:493695/1‑147 (MQ=255)
                                           gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATACGGTGAATGAGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGACCCGGAAACCGGGCGGCaaaaag                                                                                           >  2:1466444/1‑147 (MQ=255)
                                             gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGGATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGTCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGa                                                                                          >  2:70576/1‑149 (MQ=255)
                                                     tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGCAGGAATCAGaa                                                                                   >  1:1210827/1‑149 (MQ=255)
                                                     tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAt                                                                                  >  2:1338093/1‑149 (MQ=255)
                                                     tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                                                 >  2:24307/1‑149 (MQ=255)
                                                            cAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAGGCCCGCCCCGGAAAGCGGGCGGCAGAAGCAATAAAAATTACCAgt                                                                           >  1:300360/1‑148 (MQ=255)
                                                             aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGTGGGGGTAATCGGCAGGTAAGGCCGCCCCGGAAAGCGGGCGGCAGAGGGAATCAGAATTTCCAGGt                                                                          >  2:712989/1‑149 (MQ=255)
                                                                  aTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGAc                                                                     >  1:135045/1‑149 (MQ=255)
                                                                        cccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCCGAAGGAATCAGAATTACCAGGTCAGACGGGCTg                                                               >  1:919715/1‑149 (MQ=255)
                                                                         ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATACCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                              >  1:61033/1‑149 (MQ=255)
                                                                                tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCCGATAGCGGGCGGCAGAAGGAATCAGAATTTCCACGTCAGACGGGCTGCAAGTTGc                                                       >  1:137533/1‑149 (MQ=255)
                                                                                    taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGACAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGAcc                                                  >  1:1438142/1‑149 (MQ=255)
                                                                                          tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGATGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCCGACCGTTaa                                             >  1:377802/1‑149 (MQ=255)
                                                                                          tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTAAGACGGGCTGCAAGTTGCAGACCGTTaa                                             >  2:668884/1‑149 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCAAATCGGCATGGAAGGCCGCCCCGGATACCGGGGGGCAGAAGGAATAAGAATTTCCAGGTACTACCGGCGGCAATTTGCCAACAGTTAAAat                                          >  2:1516378/1‑149 (MQ=255)
                                                                                                 tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GTGGGGGGGGGTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCTGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTGGCCGACCGTTaaatacaac                                      >  1:681414/1‑143 (MQ=255)
                                                                                                 tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCATACCGTTAAAatcatc                                      >  2:53257/1‑149 (MQ=255)
                                                                                                    aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGCCCGTTAAAATCATCGGt                                   >  1:184479/1‑149 (MQ=255)
                                                                                                     aTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCCCCCCGCAAAGCGGGCGGCAGAAGGAATCAGAATTTCCCGTGCAGCCGGGCTGCAAGTTGCAGCCCGTTAAAATCATCGGtt                                  >  1:329957/1‑149 (MQ=255)
                                                                                                               tGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa                       >  1:11780/1‑149 (MQ=255)
                                                                                                                      ggTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACt                >  1:611676/1‑149 (MQ=255)
                                                                                                                       gTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGCCGGGCTGCAAGTTGCAGACCGTTAAAATCTTCGGTTGGGGTGTCGTACAACACt                >  1:769362/1‑149 (MQ=255)
                                                                                                                          ggAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCATAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTa             >  1:592768/1‑149 (MQ=255)
                                                                                                                                  aaGAGGAAA‑GGGGGGGGGGCTAATAGGCAGGGAAGGCCGCCCCGGATTGCGGGCGGCATAAGGAATCAGAATTTCAAGGTCAGACGGGCTGCAAGTTGCAGACCGTGAAAATAATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt     >  2:532583/1‑149 (MQ=255)
                                                                                                                                     agGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCAAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACATTACCTGCCGTCAg  >  1:548451/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                              
CTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAG  >  NC_000913_3_hsa_pgi/379874‑380155

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