Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,189,512 +C coding (4164/4224 nt) rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,189,5111.C96.2% 152.2 / NA 54D1388A (GAC→GCC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base . (0/2);  major base C (29/21);  minor base G (0/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑A‑C‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATT  >  NC_000913/4189370‑4189646
                                                                                                                                               |                                                                                                                                         
gcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑A‑c                                                                                                                                          >  1:934731/1‑143 (MQ=255)
       aCTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCAGTTACGCGTACCACCAGGATCGCATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCt                                                                                                                                   >  1:930417/1‑143 (MQ=255)
            gcggcgCGAAAGGGCGCGTTAGCGTGGGTCGTCTGGTCCCGGCAGGTCCCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAg                                                                                                                              <  2:877950/143‑1 (MQ=255)
                ccTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCCGAAG‑ACC‑GCCTCTGCCAGccgg                                                                                                                          >  1:408492/1‑141 (MQ=255)
                ccTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTg                                                                                                                          >  1:754637/1‑143 (MQ=255)
                 acgAAGGAGTACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGGACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTgg                                                                                                                         <  1:717552/141‑1 (MQ=255)
                  tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTgcc                                                                                                                        >  1:523407/1‑141 (MQ=255)
                         gaACGTGACCGTGGGTCGGCTGAGCACGGCAGGTACCGGTTGGGCGGACCAGCAGGATCGTATGCGTCGCCGTGCGGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAActg                                                                                                                 <  2:676254/143‑1 (MQ=255)
                          accGTTCTCGTGGGTCGTCTGCTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGCAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAActgc                                                                                                                <  2:844996/141‑1 (MQ=255)
                              ttgtctggggTCGTCCGTCCCCGGCGGGTACCGGTTGCGCGTGCCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGaa                                                                                                            <  2:742852/136‑1 (MQ=255)
                              ttaccttcggTCGTATGACCCCGGCAGGTACCGTTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGaa                                                                                                            <  1:154867/135‑1 (MQ=255)
                                 tCGTGGGTCGTCTGAGCCCGGCAGGTACCGGTGACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGc                                                                                                         <  1:317313/143‑1 (MQ=255)
                                   gTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAg                                                                                                       <  2:592580/143‑1 (MQ=255)
                                   gTGGGGCGTCTGAGCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGGGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAg                                                                                                       <  1:696026/143‑1 (MQ=255)
                                    tGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGGCGTGCTGCGGGTGAAGCTGCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGACCTGCTGCAAGCAgg                                                                                                      >  1:867546/1‑143 (MQ=255)
                                         cgtgtcaTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTgg                                                                                                 <  1:434674/137‑1 (MQ=255)
                                           tcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc                                                                                               >  1:266537/1‑143 (MQ=255)
                                           tcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCACAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc                                                                                               >  1:474696/1‑143 (MQ=255)
                                            tgggtCCCGGCGGGTAGGGGTTAGGCGTCCCGCCAGGAGCGTGTGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGGAGCGG‑A‑CGGCATGGATAATCTTGGCAGAACTGCGGGACGCAGGTCTGGGCg                                                                                              <  2:202157/139‑1 (MQ=255)
                                            cTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGGGACTGCAGGAG‑ACC‑TCATCTGCCAGCCTGGCAGAACTGCTGAACGCACGTcccggcc                                                                                              >  2:152078/1‑137 (MQ=255)
                                               aTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGCCAAGCTCCGGCTGCACGGCAGGTGACTCCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCACGCCTGGGCGCtt                                                                                           >  2:373575/1‑143 (MQ=255)
                                                   ccgtacgTACCGGGTAGGCGTACCACCAGGGTCGTATGCGGCGCCGTGCTGCGGGTGGAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa                                                                                       <  2:628931/137‑1 (MQ=255)
                                                         gtgcccggtccgcgTACCTCCTGCGTCGTATGCGTCGCCGTGCTGCAGGTGAAGCTCCGGCTGCCCCGCCGGTGCCTGCAGAAG‑AGC‑GCATCTGCCCGCCTGGCAGAACTGCTGCACGCAGGTCTGGGCGGTTCTGATAACGa                                                                                 <  2:294727/133‑1 (MQ=255)
                                                          gtccGGCTGCGCGTACCACCAGGATCGTAGGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAg                                                                                <  2:833611/141‑1 (MQ=255)
                                                               gTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGACGCTCCGGCTGCACCGCAGCTGACTGCCGAAG‑ACC‑TCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTCCTAATAACGAGtcagc                                                                           >  2:35801/1‑139 (MQ=255)
                                                               gTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGGTCCGGCTGCACCGCAGGTGACTGCCGAAG‑ACC‑CCATCTGCCAGCCTGGCAGAACTGCTCAACGCAGGTCTGGGCGGTTCTGATactgcggaatc                                                                           >  1:797954/1‑133 (MQ=255)
                                                               gTTACGCGTACCACCAGGATCGTAGGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATc                                                                           <  1:914332/143‑1 (MQ=255)
                                                                ttACGCGTACCACCAGGATCGTATGCGTCGCCGTGATGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCGGCCAGCCTGGCAGGACAGCTGAACGCAGGGCTGGGCGGCACTGAAAACGAGAAATCg                                                                          >  2:865385/1‑143 (MQ=255)
                                                                    gcgTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGCCTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTAGACGCAGGTCTGGGCGGATCTGCTGACTAATAATCGTTac                                                                      >  2:507972/1‑142 (MQ=255)
                                                                           accaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCa‑a                                                              >  1:236391/1‑143 (MQ=255)
                                                                             caGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGCGTAATCGTTAATCCGCA‑CAt                                                            >  1:373928/1‑143 (MQ=255)
                                                                             caGGAGCGCATGCGGCGCCGTGCGGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AAt                                                            <  2:590805/143‑1 (MQ=255)
                                                                                gATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAAc                                                         >  1:101824/1‑143 (MQ=255)
                                                                                gATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCAGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCCGGCAGAACCGCTGCACGCAGGTCTGGGCGGTTCCGATACAGAGTAATCGGTAACCCGCA‑aacaac                                                         >  2:401477/1‑143 (MQ=255)
                                                                                         cctCGCCCTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACGGCAGAAG‑ACC‑GCATCTGCCAGCCGGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAAcc                                                <  1:699563/141‑1 (MQ=255)
                                                                                          gTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGCAG‑ACC‑GCATCTGCCAACCTGGCAGAACTGCTGAACGCAGGCCTGGGCGGTTCTGATAACGCGGCAACGTTAATCCCCA‑AATAAGGTAAAGAccc                                               >  2:321679/1‑143 (MQ=255)
                                                                                          gTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATCAGGAGTAATCGTTCATGCGCA‑AATAACGTAAAAAccc                                               >  2:880488/1‑143 (MQ=255)
                                                                                                 tgctgcGGGTGAAGCACCGGGTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTcgg                                        <  1:101173/143‑1 (MQ=255)
                                                                                                  gctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTcggc                                       >  1:483874/1‑143 (MQ=255)
                                                                                                  gcggcGGGTGAAGCTCCGGCGGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCAGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTcggc                                       <  1:483554/140‑1 (MQ=255)
                                                                                                   ccgcggggggagggccGGCTGCACCGCAGGTGACGGCAGAAG‑AGC‑GCAGCTGCGCCTCACCCAAAACAGCAAAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTcggcg                                      <  1:213069/129‑1 (MQ=255)
                                                                                                           gAAGCTCCGGCTGGACGGCGGGTGACGGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGtttttt                              <  2:101781/143‑1 (MQ=255)
                                                                                                                tCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgg                         >  1:153421/1‑143 (MQ=255)
                                                                                                                 ccGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAAGCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATggg                        >  1:147790/1‑143 (MQ=255)
                                                                                                                     ctcgatcGCAGGTTACTGGAGAAGCC‑C‑GCATCTGCCCGACTGGCATAACTGCTGAACGCAGGTGTGGGCGGTTCTGATAACGAGTAATCGTGAATCCGTA‑AATAACGTGAAAACCCTCTTCGGCGGGTTTTTTTATGGGGGGa                    <  2:70845/137‑1 (MQ=255)
                                                                                                                     cTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGa                    <  2:930417/143‑1 (MQ=255)
                                                                                                                        cACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGATTTTTTAGGGGGGGAGtt                 >  1:825290/1‑143 (MQ=255)
                                                                                                                          ccGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa               >  2:762546/1‑143 (MQ=255)
                                                                                                                          ccGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa               >  1:151801/1‑143 (MQ=255)
                                                                                                                            gCAGGTGGCTGCATAAG‑ACC‑GCATCTGCCAGGCGGGTCGAACTGCTGAACGCAGGTCTGGGCGGTTTTGATAACGAGTAATCGTTAATTCGGA‑AATAACGTAATAACGCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAgg             <  2:182005/143‑1 (MQ=255)
                                                                                                                            gCAGGTGAGTGCAGAAG‑ACC‑GCATCTGCCAGCCGGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAgg             <  1:190881/143‑1 (MQ=255)
                                                                                                                                gTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGATTAGGGaaa         >  1:522879/1‑143 (MQ=255)
                                                                                                                                     tGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCa    >  1:454528/1‑143 (MQ=255)
                                                                                                                                     tGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGAAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCC‑AATAACGTAAAAACACGCTTCGGCGGGATTTTCTACGGGGGGAGGTTCGGGAAAGGGCa    >  2:653954/1‑143 (MQ=255)
                                                                                                                                       cAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCACCATTAC‑TACGAACCTGCTTCGGCGGGTTGTTTTATGGGGGGAGTTTAGGGAAAGAGCAtt  >  1:926296/1‑143 (MQ=255)
                                                                                                                                               |                                                                                                                                         
GCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑A‑C‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATT  >  NC_000913/4189370‑4189646

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