Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,948,346 G→A 21.1% A43V (GCG→GTG)  rho ← transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,948,3460GA21.1% 96.6 / 19.8 52A43V (GCG→GTG) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base G (23/18);  new base A (6/5);  total (29/23)
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.44e-01

GCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCAT  >  NC_002947/5948200‑5948493
                                                                                                                                                  |                                                                                                                                                     
gCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCa                                                                                                                                                     >  8:299029/1‑149 (MQ=255)
  aGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACATg                                                                                                                                                   >  5:352535/1‑149 (MQ=255)
    gTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGct                                                                                                                                                 >  8:172301/1‑149 (MQ=255)
           cggcggATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGACAGCTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCAACCGAGATCTCATCGCCGCTTTTCGCATGCTTCTTcag                                                                                                                                          >  5:77999/1‑149 (MQ=255)
             gcggATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTcagca                                                                                                                                        <  7:299029/149‑1 (MQ=255)
                gATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACATGCTTCTTCAGCAggg                                                                                                                                     <  6:352535/149‑1 (MQ=255)
                                 cATAGATATCGTCTTGGCCTGCCAGGTAGGACGCATGAGCCGAACGCAGGAAACCGAAACGATCCTGGAGAATCTCCAGCACTCCGTCACCCGAGATGTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGt                                                                                                                  <  6:172537/149‑1 (MQ=255)
                                    aGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACAt                                                                                                                                                    >  6:241057/1‑114 (MQ=255)
                                    aGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACAt                                                                                                                                                    <  5:241057/114‑1 (MQ=255)
                                    aGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTg                                                                                                               >  6:301951/1‑149 (MQ=255)
                                            tCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCt                                                                                                                    >  6:298025/1‑136 (MQ=255)
                                            tCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCt                                                                                                                    <  5:298025/136‑1 (MQ=255)
                                                ggccggccAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACgg                                                                                                      >  6:113561/1‑146 (MQ=255)
                                                ggccggccAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACgg                                                                                                      <  5:113561/146‑1 (MQ=255)
                                                 gccggccAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCa                                                                                                  >  2:154788/1‑149 (MQ=255)
                                                       cAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCCCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCAc                                                                                                                         >  3:122152/1‑120 (MQ=255)
                                                       cAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCCCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCAc                                                                                                                         <  4:122152/120‑1 (MQ=255)
                                                             ggACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGc                                                                                      <  1:154788/149‑1 (MQ=255)
                                                                       gCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGccc                                                                                    <  2:57010/141‑1 (MQ=255)
                                                                       gCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGccc                                                                                    >  1:57010/1‑141 (MQ=255)
                                                                       gCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTc                                                                            >  2:209414/1‑149 (MQ=255)
                                                                        ccGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCg                                                                           <  5:301951/149‑1 (MQ=255)
                                                                               cAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCAttt                                                                    <  1:209414/149‑1 (MQ=255)
                                                                                  gaaaccgaaaccATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTcc                                                                  >  4:78390/1‑148 (MQ=255)
                                                                                  gaaaccgaaaccATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTcc                                                                  >  5:244622/1‑148 (MQ=255)
                                                                                      ccgaaaccATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCaa                                                                >  7:218438/1‑146 (MQ=255)
                                                                                        gaaaccATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCaaaa                                                              <  8:218438/146‑1 (MQ=255)
                                                                                           accATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCgg                                                        <  7:172301/149‑1 (MQ=255)
                                                                                             cATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGt                                                       >  5:69777/1‑148 (MQ=255)
                                                                                                    gagaTTCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCt                                               >  6:237975/1‑149 (MQ=255)
                                                                                                    gagaATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAGAAAAAAGGAGTTTTTGGGCCAACGGGCAATTTTTTGGTTGCCAATATGTTGGGCAATTTACAAAATATGGTTAATggg                                                 >  1:282080/1‑149 (MQ=255)
                                                                                                        aTCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTg                                           <  6:244622/149‑1 (MQ=255)
                                                                                                         tctcCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGc                                          >  6:115915/1‑149 (MQ=255)
                                                                                                            ccAGCACGCCGTCACCCGAGATCTCTTCGTCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCtt                                        >  2:268495/1‑148 (MQ=255)
                                                                                                             cAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCgg                                                        >  4:209782/1‑131 (MQ=255)
                                                                                                             cAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCgg                                                        <  3:209782/131‑1 (MQ=255)
                                                                                                                  cGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTtca                                 <  5:237975/149‑1 (MQ=255)
                                                                                                                  cGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCaa                                                                <  4:75516/118‑1 (MQ=255)
                                                                                                                  cGCCGTCACCCGAGATCTCTTCGCCGCTTTTCACATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCaa                                                                >  3:75516/1‑118 (MQ=255)
                                                                                                                   gCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGtttt                                                                                            >  5:245712/1‑89 (MQ=255)
                                                                                                                   gCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGtttt                                                                                            <  6:245712/89‑1 (MQ=255)
                                                                                                                     cGTCAGCCGAGATCTCTTCGCCGCTGTTCGCATGCAGCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTtcagtc                              <  6:77999/149‑1 (MQ=255)
                                                                                                                      gTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTACAAAAGATCg                                                         >  6:269289/1‑121 (MQ=255)
                                                                                                                         aCCCGAGATCTCTTCGTCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAgg                           <  1:268495/148‑1 (MQ=255)
                                                                                                                                 tctcTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGa                                      >  4:105289/1‑129 (MQ=255)
                                                                                                                                 tctcTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGa                                      <  3:105289/129‑1 (MQ=255)
                                                                                                                                     ttCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAggg                 <  2:282080/146‑1 (MQ=255)
                                                                                                                                     ttCGCCGCTTTTCACATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGaa              >  1:188360/1‑149 (MQ=255)
                                                                                                                                         ccgcTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGAc          >  8:234798/1‑149 (MQ=255)
                                                                                                                                          cgcTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTc                                                                   <  2:3643/91‑1 (MQ=255)
                                                                                                                                          cgcTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTc                                                                   >  1:3643/1‑91 (MQ=255)
                                                                                                                                              tttCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCAATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTaa      >  7:104506/1‑148 (MQ=255)
                                                                                                                                                  gCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCAt  <  8:104506/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAA‑‑AATCACGTCCTGTTTGCGCGAACGGGCCATGTTTTCGATGCCCATCTGTTCGGCCATTTCCAAAAGATCGGTAATCGGCTTTTGCTTGAGTTCAGTCAGGTTCATAAGGGGAATGACGTAATCAT  >  NC_002947/5948200‑5948493

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