Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,948,251 C→T 44.5% 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,2510CT44.5% 11.8 / 49.7 46G75S (GGC→AGC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base C (14/11);  new base T (10/10);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATG  >  NC_002947/5948104‑5948399
                                                                                                                                                   |                                                                                                                                                    
ggCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggccg                                                                                                                                                     <  1:322164/149‑1 (MQ=255)
 gCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGggctgg                                                                                                                                                    <  1:328683/149‑1 (MQ=255)
   ggTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGGTCCCCTTGGGTCGGTGGGCGGATCTTGCCGAGGATGGGGTCGCCGGTACGCAGGTTGAAACGGCGGATGTGGGGGGGCGAGTGATAGATAGCGTGTGggccg                                                                                                                                                     >  2:207919/1‑146 (MQ=255)
             gATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACg                                                                                                                                        >  4:268501/1‑149 (MQ=255)
                ggTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGTCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCa                                                                                                                                      <  6:155436/148‑1 (MQ=255)
                        ccTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGaa                                                                                                                             <  5:38108/149‑1 (MQ=255)
                          ttCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACg                                                                                                                           >  5:182582/1‑149 (MQ=255)
                                    cGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGgaaacc                                                                                                                  >  6:169684/1‑148 (MQ=255)
                                             gCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATcc                                                                                                        <  1:141599/149‑1 (MQ=255)
                                             gCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACAGAAACCATcc                                                                                                        >  5:346787/1‑149 (MQ=255)
                                             gCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATAAGCCGAACGCAGGAAACCGAAACCATcc                                                                                                        >  7:65722/1‑149 (MQ=255)
                                                 cccTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATc                                                                                                         <  4:125395/144‑1 (MQ=255)
                                                 cccTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATc                                                                                                         >  3:125395/1‑144 (MQ=255)
                                                 cccTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTg                                                                                                      <  8:152417/147‑1 (MQ=255)
                                                   cttccttcGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCa                                                                                                                         <  2:146045/126‑1 (MQ=255)
                                                   cttccttcGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCa                                                                                                                         >  1:146045/1‑126 (MQ=255)
                                                      ccttcGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGaa                                                                                                 <  3:268501/147‑1 (MQ=255)
                                                           ggTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAg                                                                                          <  2:85042/149‑1 (MQ=255)
                                                                gCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGc                                                                                     >  4:260744/1‑149 (MQ=255)
                                                                         tccGACGATGGTGTCGCCGGTACGCAGGTTGAAACGTCGTAGCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGAAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCAccc                                                                             <  8:65722/147‑1 (MQ=255)
                                                                           cgacgaTGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCg                                                                                   <  4:108711/140‑1 (MQ=255)
                                                                           cgacgaTGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCg                                                                                   >  3:108711/1‑140 (MQ=255)
                                                                               gaTGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGAtct                                                                      >  4:104874/1‑149 (MQ=255)
                                                                                      tCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAg                                                                                          <  4:26801/122‑1 (MQ=255)
                                                                                      tCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAg                                                                                          >  3:26801/1‑122 (MQ=255)
                                                                                      tCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAAt                                                                                                >  8:33977/1‑116 (MQ=255)
                                                                                      tCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAAt                                                                                                <  7:33977/116‑1 (MQ=255)
                                                                                       cGCCGGTACGCAGGTTGACACGGCGGATCTGGCTGGGCGAGACATAGATATCGGCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTTGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTcgccg                                                              <  7:32134/149‑1 (MQ=255)
                                                                                         ccGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCt                                                            >  8:286470/1‑149 (MQ=255)
                                                                                          cGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAGCCGAACCCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCtt                                                           >  6:124656/1‑149 (MQ=255)
                                                                                                  aGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATg                                                   <  6:182582/149‑1 (MQ=255)
                                                                                                    gTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGct                                                 >  6:211209/1‑149 (MQ=255)
                                                                                                     ttGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGctt                                                <  7:286470/149‑1 (MQ=255)
                                                                                                        aaaCGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGcttc                                               >  6:153435/1‑147 (MQ=255)
                                                                                                        aaaCGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGcttc                                               <  5:153435/147‑1 (MQ=255)
                                                                                                         aaCGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGACACCGAAACCAGCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGccctgctg                                               >  2:274141/1‑139 (MQ=255)
                                                                                                           cggcggATCTGGCTGGGCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGcttctt                                             >  7:317078/1‑146 (MQ=255)
                                                                                                                     ggctggGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGaa                                 >  2:294351/1‑148 (MQ=255)
                                                                                                                         gggCGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCa                            <  3:260744/149‑1 (MQ=255)
                                                                                                                         gggCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCa                            >  6:351997/1‑149 (MQ=255)
                                                                                                                          ggCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCAc                           <  5:351997/149‑1 (MQ=255)
                                                                                                                            cGAGACATAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGt                         <  3:104874/149‑1 (MQ=255)
                                                                                                                               gaCATAGATATCGTCTGGGCCGGCCAGGCAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCt                      <  5:211209/149‑1 (MQ=255)
                                                                                                                                 cATAGATATCGTCTGGGCGGCCCAGGTCGCACGCATCCGCCGACCGCAGGAAACCGAACCCATCCTGGAGAATCTCCAGCACGCCGTCCCCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCTAAAATCACGTCCTGt                    <  6:359426/149‑1 (MQ=255)
                                                                                                                                  aTAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGtt                   >  5:359426/1‑149 (MQ=255)
                                                                                                                                  aTAGATATCGCCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCATGGAGAATCTCCAGCACGCCGTCACCAGAGAGCTCTTCGCCGCTTTTCGCATGGTGATTCAGCAGGGCGAAAATCACGTCCTGtt                   >  8:133673/1‑149 (MQ=255)
                                                                                                                                   tAGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGttt                  <  5:169684/149‑1 (MQ=255)
                                                                                                                                    aGATATCGTCTGGGCTGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTg                 >  1:92199/1‑149 (MQ=255)
                                                                                                                                    aGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTg                 >  2:118968/1‑149 (MQ=255)
                                                                                                                                                   cggccAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATg  >  2:316584/1‑149 (MQ=255)
                                                                                                                                                   cggccAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATg  >  4:155336/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCGGTCGAAGTTGATGGTATCGACCTTCAGCAGGGCGAAGTAACGCTCCCCTTCCTTCGGTGGGCGGATCTTGCCGACGATGGTGTCGCCGGTACGCAGGTTGAAACGGCGGATCTGGCTGGGCGAGACATAGATATCGTCTGGGCCGGCCAGGTAGGACGCATCAGCCGAACGCAGGAAACCGAAACCATCCTGGAGAATCTCCAGCACGCCGTCACCCGAGATCTCTTCGCCGCTTTTCGCATGCTTCTTCAGCAGGGCGAAAATCACGTCCTGTTTGCGCGAACGGGCCATG  >  NC_002947/5948104‑5948399

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