Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,960,789 C→G 18.3% T338R (ACA→AGA)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,960,7890CG18.3% 93.4 / 20.6 44T338R (ACA→AGA) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (14/22);  new base G (5/3);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

TGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATT  >  CP009273/3960642‑3960937
                                                                                                                                                   |                                                                                                                                                    
tgGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACa                                                                                                                                                     <  2:148872/149‑1 (MQ=255)
gggATGCCACCGCCCTGTATCTTCCGAAACGCTGCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCGTATCGATACCGGTTTTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACa                                                                                                                                                     <  2:208100/148‑1 (MQ=255)
 ggATGCCAACGCCCTGCATCGTCCGCAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATTGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGTAGTTATCTACGAAGAGTTTAAAGGTACAg                                                                                                                                                    <  2:532863/149‑1 (MQ=255)
   aTGCCACCGCCCGGCAGCGTCCGAAACGCTTCTTTGGTGCGGCGCTTAACGTGGACGAGGGCGGCAGCCTGACTATTATCGTGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGc                                                                                                                                                  <  2:534088/149‑1 (MQ=255)
     gCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTACCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTAGAGGCaa                                                                                                                                                <  2:337534/149‑1 (MQ=255)
        aaCGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGGCGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCGACAt                                                                                                                                             >  2:663831/1‑149 (MQ=255)
         aCGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGCGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTCCGAATAGTTTAAAGGTACAGGCAACATg                                                                                                                                            <  2:523098/149‑1 (MQ=255)
                gCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTAGAGGCAACATGGAACTGc                                                                                                                                     <  2:75459/149‑1 (MQ=255)
                 cATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCa                                                                                                                                    >  2:105561/1‑149 (MQ=255)
                   tCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCAcc                                                                                                                                  <  1:263997/149‑1 (MQ=255)
                    cGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTCGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAACGTACAGGCAACATGGAACTGCACct                                                                                                                                 <  2:485752/149‑1 (MQ=255)
                              gcgtcttTGGTGCTGCGCGTAAAGTGGAGGAGGGCGGCAGCATGACCATTATTGCGACGGCGATTACCGATACCGGGTGTAAAATGGACGAAGTTATATACGAAGAGTTTAAAGGTACAGGCAGAATGGAACTGCGCATCTCTCGaaac                                                                                                                       >  1:301854/4‑148 (MQ=255)
                               cttcttTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGa                                                                                                                      <  2:617508/149‑1 (MQ=255)
                                ttcttTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTAGAGGCAACATGGAACTGCACCTCTCTCGTAAGAt                                                                                                                     <  1:699324/149‑1 (MQ=255)
                                ttcttTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGAt                                                                                                                     <  2:350108/149‑1 (MQ=255)
                                 tcgttGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATGATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATc                                                                                                                    <  1:756058/146‑1 (MQ=255)
                                        tGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGaaa                                                                                                             <  2:760717/149‑1 (MQ=255)
                                             cgcgTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCTGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAcgc                                                                                                        <  1:624097/149‑1 (MQ=255)
                                                 tAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCt                                                                                                    >  2:819171/1‑149 (MQ=255)
                                                      tGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATGCCGGTTCTAAAATGGACGAAGTTGTCTACGAGGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTGAGATCGCTGAAAAACGCGTATTccag                                                                                               >  2:680962/1‑147 (MQ=255)
                                                      tGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGCGTTTAAAGGTACAGGCACCCTGGACCTGCCCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCg                                                                                               >  2:350110/1‑149 (MQ=255)
                                                          agagGGCGGCAGCCTGACCATTATCGCGACGGCGATTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTa                                                                                           >  2:441170/1‑149 (MQ=255)
                                                               gcggcgGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGAc                                                                                      <  2:101098/149‑1 (MQ=255)
                                                               gcggcAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGAc                                                                                      >  2:805053/1‑149 (MQ=255)
                                                                 ggcggcCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTa                                                                                    <  1:391150/149‑1 (MQ=255)
                                                                    aGCTTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACaa                                                                                 <  2:217370/149‑1 (MQ=255)
                                                                      ccTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACACGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAAcc                                                                               >  2:588769/1‑149 (MQ=255)
                                                                            cATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGAAAAAAAAAAAAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTg                                                                         <  1:482390/149‑1 (MQ=255)
                                                                                                 cggtACCGGTACTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTAGAGGCAACATGGAACTGCACACCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTACCGACTACAACCGTTCTGGTACCCGTAAAGACGAgctgc                                                    >  1:799803/4‑149 (MQ=255)
                                                                                                 cGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTAGAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAACAAGAgctgc                                                    >  1:722370/1‑149 (MQ=255)
                                                                                                  gATACCGGTTCTAAAATGGGCGAAGGTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctgct                                                   <  1:39747/149‑1 (MQ=255)
                                                                                                  gATACCGGTTCTAAAATGGACGAAGTTATATACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAgctgct                                                   <  1:552081/149‑1 (MQ=255)
                                                                                                       cGGTTCTAAAATGGACGAAGTTATCTAAGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGa                                              <  2:478541/149‑1 (MQ=255)
                                                                                                                   gggcGAAGTTATCTACGCAGAGTTTAAAGGTACAGGCAACATGGAGCTGCACCTCTCTCGTAAGAACACAGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAAc                                  <  1:805053/146‑1 (MQ=255)
                                                                                                                    gACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACt                                 >  1:475651/1‑149 (MQ=255)
                                                                                                                     aCGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTg                                >  2:496730/1‑149 (MQ=255)
                                                                                                                         aGTTATCTACGAAGAGTTTAAAGGTAGAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGGCTTCCCGGCTAGCGACTACAACCGTTCTGGTACCCGTAAAGAACAGCTGCTCACGACGCAGGAAGAACGACAGa                            >  2:557662/1‑149 (MQ=255)
                                                                                                                               cTACGAAGCTTTTAAATGCACAGGCAACATACCACCGCACCTCCCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtgt                      <  2:162457/149‑1 (MQ=255)
                                                                                                                                tACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtgtg                     <  1:105561/149‑1 (MQ=255)
                                                                                                                                 aCGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGGACTGCAGAAAATGTgg                    >  1:37717/1‑149 (MQ=255)
                                                                                                                                 aCGAAGAGCTTAAAGGTAGAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCCGGAAGAACTGCAGAAAACGTgg                    >  1:211090/1‑149 (MQ=255)
                                                                                                                                    aaGAGTTTAAAGGTAGAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAAGCGTCTTCCCGGCTAGCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGGAGAACTGCAGAAAATGTGGATc                 >  2:334949/1‑149 (MQ=255)
                                                                                                                                         tttAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgcg            >  1:408765/1‑149 (MQ=255)
                                                                                                                                           tAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAAAGCGTCTTACAGGCTATCGAATACAACTGTCCCGTGACCCGTAAAGAAGAACTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCa          >  2:854050/1‑149 (MQ=255)
                                                                                                                                                   cAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCAtt  <  2:475651/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGATGCCAACGCCCTGCATCGTCCGAAACGCTTCTTTGGTGCGGCGCGTAACGTGGAAGAGGGCGGCAGCCTGACCATTATCGCGACGGCGCTTATCGATACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATT  >  CP009273/3960642‑3960937

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