Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,948,205 T→C 15.0% N90S (AAC→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,2050TC15.0% 89.8 / 8.8 41N90S (AAC→AGC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (17/17);  new base C (2/4);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

CAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCG  >  NC_002947/5948057‑5948346
                                                                                                                                                    |                                                                                                                                             
cAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTCCCGACGATGGTGTCGCCGGTACGCAGGt                                                                                                                                               <  1:272217/149‑1 (MQ=255)
   cggcgTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGCTGa                                                                                                                                            <  1:45659/149‑1 (MQ=255)
          aGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAAcggc                                                                                                                                      <  3:86093/148‑1 (MQ=255)
            gTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGa                                                                                                                                   >  5:118682/1‑149 (MQ=255)
               ttCGAACAGGTTCTTGTTCTTCTCGTTTTGCGGGCGGTCGATGTTGATGGTATCGACCTTCAGCAGGGCGATGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGGTGAAACGTCGGATc                                                                                                                                 <  8:283634/148‑1 (MQ=255)
                   aaCAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATctggct                                                                                                                            >  6:220278/1‑149 (MQ=255)
                   aaCAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGATGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATctggct                                                                                                                            >  2:122714/1‑149 (MQ=255)
                      aGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGtt                                                                                                                                              >  7:114697/1‑128 (MQ=255)
                      aGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGtt                                                                                                                                              <  8:114697/128‑1 (MQ=255)
                       ggATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGc                                                                                                                        <  5:243542/149‑1 (MQ=255)
                       ggAGCTTGTTCTTCGCGTTTTCCGCGCTGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGTTACGCAGGTTGAAACGGCGGATCTGGCTGGGc                                                                                                                        <  2:134594/149‑1 (MQ=255)
                              gTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGAc                                                                                                                   <  8:202752/147‑1 (MQ=255)
                               ttcttcGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATa                                                                                                                <  5:220278/149‑1 (MQ=255)
                               ttcttcGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGATGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATa                                                                                                                <  1:122714/149‑1 (MQ=255)
                                  ttcGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGa                                                                                                              >  5:322002/1‑148 (MQ=255)
                                     gcgTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGTTGGGCGAGACATAGATATc                                                                                                          <  6:118682/149‑1 (MQ=255)
                                       gTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTTGGCGGATCTTGCCGAAGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGt                                                                                                        >  5:196477/1‑149 (MQ=255)
                                                    tCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAg                                                                                           <  7:144435/149‑1 (MQ=255)
                                                              tGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGCTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCa                                                                                 <  1:80948/149‑1 (MQ=255)
                                                                   tCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCAGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCCGc                                                                            >  8:224538/1‑149 (MQ=255)
                                                                        cTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGCTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAAc                                                                       <  3:324909/149‑1 (MQ=255)
                                                                               agGTCGAAGTAACGCTCCCCTTCCTTCGGTTGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGTGGGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGga                                                                 <  7:224538/148‑1 (MQ=255)
                                                                                gggCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGgaaa                                                               >  8:124234/1‑149 (MQ=255)
                                                                                gggCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGgaa                                                                <  6:17726/148‑1 (MQ=255)
                                                                                gggCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGgaa                                                                >  5:17726/1‑148 (MQ=255)
                                                                                 ggCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCAt                                                                                <  5:178246/131‑1 (MQ=255)
                                                                                 ggCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCAt                                                                                >  6:178246/1‑131 (MQ=255)
                                                                                          aCGCTCCACTTCCCTCGGTGGGCGGATCTTGACGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCAACAGCCGAACGCAGGAAACCGAAACCAt                                                     >  8:182222/1‑149 (MQ=255)
                                                                                             cTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGCTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCt                                                  >  3:213747/1‑149 (MQ=255)
                                                                                                cccTTGCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGAGGCATCAGCCGCACGCAGGAGGCCGAAAGCGTCCTg                                                 >  4:156953/1‑147 (MQ=255)
                                                                                                cccTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGGAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGga                                               <  6:196477/149‑1 (MQ=255)
                                                                                                cccTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGga                                               <  7:36228/149‑1 (MQ=255)
                                                                                                  cttccttcGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGgag                                              >  3:13126/1‑148 (MQ=255)
                                                                                                       gtcgggcggCGGATCTTGCCGACGACGGTGTCGCCGGTACGCAGGCTGAACCGGCGGAGCTGGCTCGGCGAGACATAGACATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGACACCGCAACCATCCTGGAGAAtctc                                        <  4:213747/142‑1 (MQ=255)
                                                                                                            tGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCa                                   >  4:325865/1‑149 (MQ=255)
                                                                                                                 ggATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCg                              >  3:274087/1‑149 (MQ=255)
                                                                                                                          cgacgaTGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCAccc                        >  2:10888/1‑146 (MQ=255)
                                                                                                                                  gtgtCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCtt              <  6:322002/148‑1 (MQ=255)
                                                                                                                                       gCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTcgccgc        >  7:63794/1‑149 (MQ=255)
                                                                                                                                       gCCGGTACGCAGGCTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTcgccgc        >  1:159869/1‑149 (MQ=255)
                                                                                                                                            tACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTc   <  4:274087/149‑1 (MQ=255)
                                                                                                                                              cGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCg  <  8:63794/148‑1 (MQ=255)
                                                                                                                                                 aGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggcc                                                                                                 >  4:357740/1‑50 (MQ=255)
                                                                                                                                                 aGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggcc                                                                                                 <  3:357740/50‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
CAGCGGCGTCAGGTTTTCGAACAGGATCTTGTTCTTCGCGTTTTCCGGGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCG  >  NC_002947/5948057‑5948346

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: ‑