Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,186,577 G→C D410H (GAT→CAT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,5770GC98.2% 199.7 / ‑5.9 57D410H (GAT→CAT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (29/27);  minor base T (0/1);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCA  >  NC_000913/4186434‑4186717
                                                                                                                                               |                                                                                                                                            
gTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATc                                                                                                                                         <  1:133247/149‑1 (MQ=255)
                    ccGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCcgcg                                                                                                                      >  8:91500/1‑148 (MQ=255)
                         gaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGa                                                                                                                     >  3:104834/1‑144 (MQ=255)
                         gaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGa                                                                                                                     <  4:104834/144‑1 (MQ=255)
                                  aCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAg                                                                                                                               >  5:109293/1‑125 (MQ=255)
                                  aCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAg                                                                                                                               <  6:109293/125‑1 (MQ=255)
                                      gAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGaa                                                                                                   >  4:141727/1‑149 (MQ=255)
                                               aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTg                                                                                                                                      >  4:169678/1‑105 (MQ=255)
                                               aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTg                                                                                                                                      <  3:169678/105‑1 (MQ=255)
                                               aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                                                          >  1:98762/1‑149 (MQ=255)
                                                 aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGc                                                                                             >  5:95311/1‑144 (MQ=255)
                                                 aCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGc                                                                                             <  6:95311/144‑1 (MQ=255)
                                                    gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgc                                                                                                       >  5:100864/1‑131 (MQ=255)
                                                    gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTActgc                                                                                                       <  6:100864/131‑1 (MQ=255)
                                                    gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGActc                                                                                     <  7:91500/149‑1 (MQ=255)
                                                    gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGActc                                                                                     <  1:44372/149‑1 (MQ=255)
                                                         tCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCAc                                                                                >  2:44821/1‑149 (MQ=255)
                                                          cTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTgg                                                                                                                                     <  2:56747/95‑1 (MQ=255)
                                                          cTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTgg                                                                                                                                     >  1:56747/1‑95 (MQ=255)
                                                           tACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGa                                                                                                                                    >  7:78806/1‑95 (MQ=255)
                                                           tACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGa                                                                                                                                    <  8:78806/95‑1 (MQ=255)
                                                             cGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCAcc                                                                               <  1:44821/146‑1 (MQ=255)
                                                                  aGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                                        >  7:70507/1‑148 (MQ=255)
                                                                   gCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                         <  8:99490/146‑1 (MQ=255)
                                                                   gCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                         >  7:99490/1‑146 (MQ=255)
                                                                        aaCTGCGTGGTCTTGCTACGACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAg                                                                 <  3:141727/149‑1 (MQ=255)
                                                                          cTGCGTTGGATTGCTACTACCATTTAATCTGCGTAGAAAATGGTTGTGCGCGAAGAAGCTGTCGTTTGGTATATCCTGGAAGAAGTTATCCGCGAACACCCGGTAGTGGTTAAACTTGCACCGACTCTTCAGCGTCTgg                                                                         <  8:121046/139‑1 (MQ=255)
                                                                          cTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTTGCATATCCAGTACGAAGTTATCCGCGAACACCCCTTACGGCTGAACCGTTCACCGACTCTGCACCGTCggg                                                                         >  7:121046/1‑139 (MQ=255)
                                                                                 gTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGaa                                                        <  2:98762/149‑1 (MQ=255)
                                                                                  tCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                                             >  5:154322/1‑143 (MQ=255)
                                                                                  tCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                                             <  6:154322/143‑1 (MQ=255)
                                                                                  tCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                       <  8:70507/149‑1 (MQ=255)
                                                                                      gCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCgg                                                    >  4:84746/1‑148 (MQ=255)
                                                                                        tACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                                                 >  6:9750/1‑149 (MQ=255)
                                                                                         accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACAcc                                                                                                                >  8:11632/1‑85 (MQ=255)
                                                                                         accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACAcc                                                                                                                <  7:11632/85‑1 (MQ=255)
                                                                                          ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                                                 >  5:139907/1‑147 (MQ=255)
                                                                                          ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc                                                 <  6:139907/147‑1 (MQ=255)
                                                                                           caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                         <  5:31845/122‑1 (MQ=255)
                                                                                           caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                         >  6:31845/1‑122 (MQ=255)
                                                                                             ccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                   <  4:92325/126‑1 (MQ=255)
                                                                                             ccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                   >  3:92325/1‑126 (MQ=255)
                                                                                               aTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATAACATGGACGAAGTTATCCGCGAACACCCGGTACAGCAGAACCGTGCACCGACTCAGCACCCTCTGGGTATCCAGGCATTTGAAAAGGAACTGATCg                                           >  8:108767/1‑148 (MQ=255)
                                                                                                     gCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTaa                                    <  5:9750/149‑1 (MQ=255)
                                                                                                      cTGCGAAGAAAATGGTTGCGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGAAGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                          >  2:121058/1‑142 (MQ=255)
                                                                                                      cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa                                          <  1:121058/142‑1 (MQ=255)
                                                                                                        gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                               >  6:112114/1‑135 (MQ=255)
                                                                                                        gCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                               <  5:112114/135‑1 (MQ=255)
                                                                                                              aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCt                                                                           <  3:94710/101‑1 (MQ=255)
                                                                                                              aaaaTGGTTGAGCGAGAAGAAGCTGTCGTTTGGCATATCCTGGACTAAGTTATCCTCGAACACACGGTACTGCTGAACCGTGCACCGACTCTGCACCGACt                                                                           >  4:94710/1‑101 (MQ=255)
                                                                                                                aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTc                                                                            <  3:73882/98‑1 (MQ=255)
                                                                                                                aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTc                                                                            >  4:73882/1‑98 (MQ=255)
                                                                                                                 aTGGTTGAGCGCGAAGAAGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                        >  7:57954/1‑149 (MQ=255)
                                                                                                                             gaagaaGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCa                                                                                 <  2:63330/80‑1 (MQ=255)
                                                                                                                             gaagaaGCTGTCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCa                                                                                 >  1:63330/1‑80 (MQ=255)
                                                                                                                                       tCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGc   <  8:57954/148‑1 (MQ=255)
                                                                                                                                       tCGTTTGGCATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCa  >  7:99162/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                            
GTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCA  >  NC_000913/4186434‑4186717

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