Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,175,524 G→T V83F (GTC→TTC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,175,5240GT95.7% 132.0 / NA 46V83F (GTC→TTC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (11/33);  total (12/34)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCG  >  CP009273/4175382‑4175670
                                                                                                                                              |                                                                                                                                                  
cGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTGCAAGCGCCTGAAACACCGTGGCTTCATCt                                                                                                                                              <  1:659771/149‑1 (MQ=255)
  gTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTc                                                                                                                                                  >  2:296944/1‑143 (MQ=255)
      aGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCt                                                                                                                                              <  2:506969/143‑1 (MQ=255)
       gTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGGGGCGTCAACTGTgagga                                                                                                                                       >  1:256486/1‑147 (MQ=255)
        ttAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCtgt                                                                                                                                            <  2:311197/143‑1 (MQ=255)
        ttAAAAAGCCGGAAACCATACACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCtgt                                                                                                                                            <  2:287916/143‑1 (MQ=255)
         tAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCtgtg                                                                                                                                           <  2:42021/143‑1 (MQ=255)
          aaaaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAAGCTGTAGACCTGTTCTGCGAGCATAGCTTGGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGGGAAGTg                                                                                                                                    >  1:188547/1‑149 (MQ=255)
            aaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTgcg                                                                                                                                  >  1:523395/1‑149 (MQ=255)
                  ggAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACAACGTGGCTTAAGCTGTGAGGAGTgcg                                                                                                                                  >  2:18359/1‑143 (MQ=255)
                    aaaCCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGaa                                                                                                                          <  1:538710/149‑1 (MQ=255)
                    aaaCCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGaa                                                                                                                          <  1:684084/149‑1 (MQ=255)
                         aTCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGa                                                                                                                           <  2:451010/143‑1 (MQ=255)
                          tCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGaa                                                                                                                          <  2:237316/143‑1 (MQ=255)
                          tCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGAcc                                                                                                                    <  1:710904/149‑1 (MQ=255)
                              cTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGGGCCGGTGCGGTAAGTACAAGCGCCTGAAACACCGTGTCTTCGTCTGTGAGAAGTGCGGCGTTGAAGTTACCCAGa                                                                                                                >  1:31348/1‑149 (MQ=255)
                                   gTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTaaa                                                                                                           <  1:603622/149‑1 (MQ=255)
                                   gTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCAGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTaaa                                                                                                           <  1:11512/149‑1 (MQ=255)
                                    tACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAg                                                                                                          >  1:477391/1‑149 (MQ=255)
                                              ccAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGa                                                                                                <  1:677178/149‑1 (MQ=255)
                                               cAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGTCTGTGCGGTGAGTGCAAGCGCCTGAAACGCCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACga                                                                                                     >  2:28634/1‑142 (MQ=255)
                                               cAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACCAAAGTACGc                                                                                                     >  2:472399/1‑143 (MQ=255)
                                                     gTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAg                                                                                               <  2:487917/143‑1 (MQ=255)
                                                      tGACGGCCTTTTTTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATgg                                                                                        <  1:80125/149‑1 (MQ=255)
                                                         cGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTCAAGTACGCCGTGAGCGTATGGGcc                                                                                     <  1:506296/149‑1 (MQ=255)
                                                              tgtgcTGCGCCCGTATCGTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGc                                                                                      <  2:564855/139‑1 (MQ=255)
                                                               tttCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGcc                                                                                     <  2:563439/143‑1 (MQ=255)
                                                                   tGCGCCCTTGTCGTTGGGGCGGGAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCTTCTCTTCATCTCTTAGAAGTTCCGCGTTGAATTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACAt                                                                                 <  2:272183/143‑1 (MQ=255)
                                                                     cgcCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGCATGGGCCACATCGAACTgg                                                                         <  1:333904/149‑1 (MQ=255)
                                                                         tgTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTTGGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACt                                                                           <  2:399414/142‑1 (MQ=255)
                                                                           tATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTgg                                                                         <  2:585473/143‑1 (MQ=255)
                                                                                     ccGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGca                                                         >  1:769463/1‑149 (MQ=255)
                                                                                              gATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACAt                                                      <  2:406071/143‑1 (MQ=255)
                                                                                               aTTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTc                                               >  1:211467/1‑149 (MQ=255)
                                                                                               aTTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGCGGCTTATTCGGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTc                                               <  1:23145/149‑1 (MQ=255)
                                                                                                     aGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCATATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGaaa                                         <  1:703661/149‑1 (MQ=255)
                                                                                                      gTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGGGAGCGTATGGGCCACATCGAACTGGCTTCCCGCACTCCCCCGCTCTGGTTCCTGaaac                                        >  1:715775/1‑148 (MQ=255)
                                                                                                       tgccttcacGGTAAGCACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTCACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCt                                             <  2:50770/135‑1 (MQ=255)
                                                                                                          cTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAACTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGaa                                          <  2:572301/143‑1 (MQ=255)
                                                                                                                  tAAGTACAAGCGCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTccc                            <  1:352086/149‑1 (MQ=255)
                                                                                                                        cAAGCCCCTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCg                      <  1:455509/149‑1 (MQ=255)
                                                                                                                              ccTGAAACACCGTGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTctgc                <  1:650922/149‑1 (MQ=255)
                                                                                                                              ccTGAAACACCGGGGCTTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCATACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCg                      <  2:465074/143‑1 (MQ=255)
                                                                                                                                tGAAACACCGTGGCTTCATCTGTGAGACGTGCGGCGTTGAAGTGCCCCAGCCTAAAGTACGCCGTGATCGGAGGGGCCACATGGAACTGGCTGCCCCGACGGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGt                    <  2:168426/143‑1 (MQ=255)
                                                                                                                                           ggCTTCATCTGTGAGAAGTGCGGCGTGGAAGTGACCCACACTAATGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGa         <  2:173570/143‑1 (MQ=255)
                                                                                                                                            gCGTCATCTGTAAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGCCTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCg  <  1:435273/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CGGTGAAGTTAAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCG  >  CP009273/4175382‑4175670

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: