Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,960,101 C→A 13.4% R109S (CGC→AGC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,960,1010CA13.4% 77.6 / 3.1 35R109S (CGC→AGC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (18/11);  new base A (3/2);  total (21/14)
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

TTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTT  >  CP009273/3959965‑3960246
                                                                                                                                        |                                                                                                                                                 
tttCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGATTCCCCTAGCCAAATCCGCAGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCACCGAAAGAAGGTGAACGCTATTTTgcgc                                                                                                                                       >  1:121115/1‑149 (MQ=255)
    ctccGTTCCGCAGACAGCTCCTCCCTCGCCGGTCCTGAAGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGACAGACGGGGAACGCTATGTTGCgctgct                                                                                                                                   >  1:81506/1‑149 (MQ=255)
    ctccGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTACAAACATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCGGGTGAGAATCGCCCGCCTAAAGAAGAGGAACGCTATATTGCgctgct                                                                                                                                   >  1:104160/1‑149 (MQ=255)
         ttCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCGATTTTGCGCTGCTGaaaa                                                                                                                              >  2:282099/1‑148 (MQ=255)
           ccGCAGCCGCCTCCTACCCCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCCCCGCACCGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                                                                                            <  2:21059/149‑1 (MQ=255)
                  ccgCTCCTACCTCGCCGGTCCTGAGGACATCTACCTTTCCTCTACCCAATGCCGCCGTTTCAACCCCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                                                                                     <  2:184631/147‑1 (MQ=255)
                  catcTCCTGCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAACTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGGTGACGAAg                                                                                                                     <  2:201564/146‑1 (MQ=255)
                  cAGCTCCTACCTCGCCGGTCCTTATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTCGTAAGATTCGCCCGCCGAAAGAAGGTGAAAGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                                                                                     <  2:66579/149‑1 (MQ=255)
                  cAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCTCGCCGAAAGAAGGTGAAAGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                                                                                                     <  1:228866/149‑1 (MQ=255)
                              cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGGAGGTGAACGCTATTTTGCGCTGCGGACAGTTAACGAAGTTAACTTCGaga                                                                                                         >  2:209527/1‑147 (MQ=255)
                                ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAATGTTTTTTTGCGCTGCTGAAAGTGAACGAAGTTAACTTCGACaaa                                                                                                       <  1:19651/149‑1 (MQ=255)
                                   gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCAGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACCAACCt                                                                                                    >  1:266266/1‑149 (MQ=255)
                                      cTGATGACATCTACGTTTCCCCTGGCCAACTCAGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaa                                                                                                 <  2:168722/149‑1 (MQ=255)
                                                      ttCCCCTAGCCAAATCAGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACaaaa                                                                                 <  1:189681/149‑1 (MQ=255)
                                                             aGCCAAATTCGCCGTGTCAAACTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAAAGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt                                                                          >  2:252230/1‑149 (MQ=255)
                                                               ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGATAACGAAGTTAACTCCGACAAACCAGAAACCGCCAGCAACAAAAACACCTTTg                                                                        >  2:291793/1‑149 (MQ=255)
                                                                 aaaTCCGCCGTTTCAACCTCCGCCCTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGGAAACGCCCGCAACAAAGTCCTCTTTgag                                                                      >  1:277682/1‑149 (MQ=255)
                                                                    tCCGCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGATGTTAATTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc                                                                   <  1:19381/149‑1 (MQ=255)
                                                                    tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTTTCGTTGTTTATTTTGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAAc                                                                   <  1:19544/149‑1 (MQ=255)
                                                                     ccgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGATATTTTGCTCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAAGAAAATCCTCTTTGAGAACt                                                                  >  1:41231/1‑149 (MQ=255)
                                                                       gccgTTTCAACCTCCGCACTGGTGATACCATCTCTCGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTTACAGTTAACGAAGTTAACTTCGACAAACCTGAAGACGCCCGAAACCAAATCCACTTTGAGAACTTa                                                                >  2:202865/1‑149 (MQ=255)
                                                                           tttCAACCTCCGCACTGGTGATACCATCACTGGTAAGATTGGCCCGCCGAAGGAAGGAAGACGCTATTTTGCGTTGATGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCACCAAAATCCTCTTTGAGAACTTAAccc                                                            >  2:24684/1‑149 (MQ=255)
                                                                             tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCg                                                          <  2:278019/149‑1 (MQ=255)
                                                                               aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAAAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGATAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCCCt                                                        >  1:67225/1‑149 (MQ=255)
                                                                                         cTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAAGTGAAAGAAATTTTACGATGCTGAAAGTTAAAGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATACTCTTTGAGAACTTAACCCCGCGGCACGCAAAc                                              >  1:2600/1‑149 (MQ=255)
                                                                                           ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAActc                                            >  1:231233/1‑149 (MQ=255)
                                                                                           ggTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGAACGCAAActc                                            >  1:231231/1‑149 (MQ=255)
                                                                                                             aaGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGCGAACTTACCCCCGCTGCACGCAAACTCTCGTCTGCGtgcgcaacg                          >  2:273145/1‑141 (MQ=255)
                                                                                                              aGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGAACGCAAACTCTCGTCTGCGTATGGAACGt                         >  2:75754/1‑149 (MQ=255)
                                                                                                                     cccgccGAAAGAAGGTGAAATCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCCAACTCTCGTCTGCGCATGGCACGTggccacg                  >  2:297797/1‑144 (MQ=255)
                                                                                                                      ccgccgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACgg                 >  1:97973/1‑149 (MQ=255)
                                                                                                                           gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTa            <  2:121115/149‑1 (MQ=255)
                                                                                                                           gAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTGAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGCATGGAACGGGGTAACGGtttgc            >  2:73853/1‑146 (MQ=255)
                                                                                                                                 agGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGTTGTACGCTATCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAg      <  1:261030/149‑1 (MQ=255)
                                                                                                                                     gAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCTTATTTTCGTCTGCGTATGGAACGTGGGAACGGTTCTACTGAAGAttt  <  1:5011/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
TTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTT  >  CP009273/3959965‑3960246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: