Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,948,251 C→T G75S (GGC→AGC)  rho ← transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,948,2510CT92.9% 81.1 / ‑4.4 29G75S (GGC→AGC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (14/12);  minor base G (1/1);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCG  >  NC_002947/5948107‑5948386
                                                                                                                                                |                                                                                                                                       
ggTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggctggc                                                                                                                                      <  5:143728/148‑1 (MQ=255)
  tCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGATGGCCa                                                                                                                                    >  5:16218/1‑148 (MQ=255)
    gAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAgg                                                                                                                                  <  7:142748/148‑1 (MQ=255)
     aaGTTGATGGTATCGCCCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCGGGCCAg                                                                                                                                   >  1:181966/1‑146 (MQ=255)
     aaGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGGCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAgg                                                                                                                                  <  1:51936/147‑1 (MQ=255)
           aTGGTATCGACCTTCAGCAGGGCGAAGTACCGCTCCCCTTCCTGCGGTGGGCGGATCTTGCCGACGATGGTCTCGCCGGGCCGCAGCTTGACACGGCGGATCGGGCTGGGCGAGACAGAGATCACGTCTCGCCCGGCCAGGTAgtccg                                                                                                                           >  5:21802/1‑144 (MQ=255)
           aTGGTAGCGACCTTCAGCAGGGCGAAGTACCGCTCCCCTTCCTTCGGGGGGCGGAGCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATCGATCTCGTCGGGGCGGGCCAGGTAGGAc                                                                                                                            <  6:77983/147‑1 (MQ=255)
            tGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGc                                                                                                                          >  7:160184/1‑148 (MQ=255)
                 tCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAg                                                                                                                     >  7:113786/1‑148 (MQ=255)
                    aCCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGcc                                                                                                                   >  1:49342/1‑147 (MQ=255)
                               ggCGAGGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGGTGAAGCGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCCGCCGAACGCAGgaaa                                                                                                       >  8:12071/1‑148 (MQ=255)
                                gCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGgaaac                                                                                                      <  2:49342/148‑1 (MQ=255)
                                     gTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCAGCAGCCGAACGCAGgaaaccgaaa                                                                                                 <  7:131409/148‑1 (MQ=255)
                                      tAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGGGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGCCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGgaaaccgaaac                                                                                                <  1:126615/148‑1 (MQ=255)
                                          gCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATc                                                                                            >  3:65790/1‑148 (MQ=255)
                                            tCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCt                                                                                          <  4:64017/148‑1 (MQ=255)
                                            tCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTGCACACGGCGGATCTGGGGGGGCGAGACATAGATTTCGTCGGGGGTGGCCAGGTTGGCCGCAGCAGCGGAACGCAGgaacccgaaccctgcct                                                                                          >  4:57755/1‑143 (MQ=255)
                                                         gTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCa                                                                              <  1:118098/147‑1 (MQ=255)
                                                              cGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGc                                                                        >  5:110531/1‑148 (MQ=255)
                                                              cGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGc                                                                        >  1:9682/1‑148 (MQ=255)
                                                                         gacgaTGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGTTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCAcc                                                                 >  5:61595/1‑144 (MQ=255)
                                                                             aTGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGAtct                                                         <  2:153016/148‑1 (MQ=255)
                                                                                 tgtCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTc                                                     >  7:423/1‑148 (MQ=255)
                                                                                              cAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCa                                        <  8:423/148‑1 (MQ=255)
                                                                                                      aaCGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGcttct                                 >  3:63240/1‑147 (MQ=255)
                                                                                                      aaCGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGcttct                                 >  3:153421/1‑147 (MQ=255)
                                                                                                        cggcggATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTca                              <  4:90219/148‑1 (MQ=255)
                                                                                                                             aCATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGAAGTAAACCGAAACAATCCTGGATAATCTCCAGCACGCAGTCACCCGAGATCTCTTCGCCGCTTTTCGAATGCTTCTTCAGCAGGGCGAAAATCACGTCCt         <  8:113786/148‑1 (MQ=255)
                                                                                                                                   atatCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTgc   <  3:64123/148‑1 (MQ=255)
                                                                                                                                    tatCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTgcg  >  6:147214/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                       
GGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCG  >  NC_002947/5948107‑5948386

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: