Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,176,029 C→A P251Q (CCG→CAG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,176,0290CA93.9% 133.1 / 0.8 49P251Q (CCG→CAG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (21/25);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

AGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTA  >  CP009273/4175886‑4176169
                                                                                                                                               |                                                                                                                                            
agagCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCATTCAGTCTGGTAACACAAAAGTTGCTATGTTCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTgg                                                                                                                                         >  1:856065/1‑149 (MQ=255)
  tgCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGtt                                                                                                                                       <  1:768040/148‑1 (MQ=255)
    cTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCt                                                                                                                                           <  2:626257/143‑1 (MQ=255)
       aaCGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGATAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGt                                                                                                                                        >  2:500035/1‑143 (MQ=255)
         cGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGAGCTGCGTCAGCTGGTTc                                                                                                                                      >  2:704622/1‑143 (MQ=255)
            aaCCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGc                                                                                                                                   <  2:77069/143‑1 (MQ=255)
                aaCTCCGAAACCAAGCGTAAAAGGCTGACCAAGCGTATCAAACTGCTGGGAGCGTTCGTTCAGTCTGGTAACATACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGCTCTGCGTCAGCTGATTCCGCTGGa                                                                                                                               <  2:243185/143‑1 (MQ=255)
                     cGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGAtggtg                                                                                                                          >  2:267426/1‑143 (MQ=255)
                        aaCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGGGGTc                                                                                                                       >  2:183926/1‑143 (MQ=255)
                         aCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCGGTCGGGTAACAAACCCGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTcgc                                                                                                                <  1:213226/149‑1 (MQ=255)
                         aCCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTcgc                                                                                                                <  1:771113/149‑1 (MQ=255)
                           cAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGa                                                                                                              <  1:538433/149‑1 (MQ=255)
                              gCGTAAAAAGCTGACCAAGCGTATCAACCTGCTGGGAGCGTTCGGTCAGTCGGGTAACAAACCAGAGTGGATGTTGCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGGTTcg                                                                                                                 <  2:540865/143‑1 (MQ=255)
                                  aaaaaGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTgg                                                                                                       >  1:730368/1‑148 (MQ=255)
                                    aaaGCTGACCAAGGGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGAcc                                                                                                     >  1:628073/1‑149 (MQ=255)
                                      aGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCTTGTCCGTTTTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTg                                                                                                   <  1:82654/149‑1 (MQ=255)
                                          gACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGCCCTGAACg                                                                                               >  1:541786/1‑149 (MQ=255)
                                           aCCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACCGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGa                                                                                              <  1:279208/149‑1 (MQ=255)
                                                   tATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTAtc                                                                                      <  1:404258/149‑1 (MQ=255)
                                                    aTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAACTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTAtcg                                                                                     <  1:910115/149‑1 (MQ=255)
                                                       aaaCTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACAAGAGCGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCACCTGGTTCCGCTGGAGGGTGGTCGTTTCGCGCCTTCTGACCttcacgcgctctt                                                                                        >  2:164098/1‑131 (MQ=255)
                                                           tgctgGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTc                                                                              >  1:286156/1‑149 (MQ=255)
                                                              tgggatcttgcggtcAGTCTGGTAACAAACCCGAGTGGATGACCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCCGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCAtt                                                                           <  1:14758/136‑1 (MQ=255)
                                                              tgGAAGCGTTCGTCCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTcgc                                                                                 <  2:829807/143‑1 (MQ=255)
                                                                gAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGCTCTGTATCGTCGCGt                                                                               >  2:227018/1‑143 (MQ=255)
                                                                gAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTaa                                                                         <  1:858154/149‑1 (MQ=255)
                                                                gAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGACGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTaa                                                                         <  1:504622/149‑1 (MQ=255)
                                                                  aGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGACCGATCTGTATCGTCGCGTCATTAAcc                                                                       >  1:269493/1‑149 (MQ=255)
                                                                     gttcgttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTa                                                                          >  2:96915/1‑143 (MQ=255)
                                                                            cagttcggTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGCCAAGGACGTTCTGTATCGTCGCGTCATTAACCGTaa                                                                   <  2:381271/137‑1 (MQ=255)
                                                                               tCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTaacaa                                                                <  2:427859/143‑1 (MQ=255)
                                                                                  ggTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCg                                                             >  2:155264/1‑143 (MQ=255)
                                                                                    caacaacCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACACATTTTAGCGTCGCGTCATTAACCGTAACAACCGTc                                                           <  2:809270/142‑1 (MQ=255)
                                                                                     aaCAAACCAGAGTGGATGATCCGGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCCGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCt                                                          <  2:747621/143‑1 (MQ=255)
                                                                                          aCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAAc                                                     >  2:355762/1‑143 (MQ=255)
                                                                                            tagagTGGATGATCATGACCGTTCTGCCGGTCCTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTctgctg                                             <  1:296598/148‑1 (MQ=255)
                                                                                            cagggtgcactctCCTGCCCGTTCCGCCGGTACTGCCTCCCGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGt                                                   <  2:423906/131‑1 (MQ=255)
                                                                                                    atgatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGCTCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAGCGTCTGCTGGa                                           >  2:209745/1‑143 (MQ=255)
                                                                                                      gatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGTGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCATCTGAAACGTCTGCTGGATCTGGCTg                                   >  1:704277/1‑149 (MQ=255)
                                                                                                      gatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCGGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATc                                         >  2:362339/1‑143 (MQ=255)
                                                                                                      gatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCAGGTTCCGCTGGATGGTGGGCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTGACCGTAACAACCGTCTGAAAAGTCTGCTGAATCTAGCTg                                   >  1:533100/1‑149 (MQ=255)
                                                                                                       attcTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGCCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTgc                                  <  1:240626/146‑1 (MQ=255)
                                                                                                       atCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTACCAACCGTCTGAAACGGCTGCTGGATCt                                        >  2:331459/1‑143 (MQ=255)
                                                                                                             aCCGTTCTGCCGGCACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGa                            <  1:655916/149‑1 (MQ=255)
                                                                                                                       cGGTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATTGTCGTGTCATTTTCCGTAACTACCGGCTGAAACGTTTGCTGGATCTGGCTGCGCGGGTCTTCATcgtt                  >  1:266461/1‑148 (MQ=255)
                                                                                                                         gTACTGCCGCCAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcgtacg                >  1:148013/1‑149 (MQ=255)
                                                                                                                                   cAGATCTGCGTCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCAGGACATCATCGTACGTAACGaaaaa      <  1:460297/149‑1 (MQ=255)
                                                                                                                                      aTCTGCGCCAGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGt   <  1:250403/149‑1 (MQ=255)
                                                                                                                                             tCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTa  >  2:460459/1‑143 (MQ=255)
                                                                                                                                               |                                                                                                                                            
AGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTA  >  CP009273/4175886‑4176169

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