Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,186,086 C→T P246L (CCG→CTG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,186,0860CT97.2% 108.6 / ‑5.6 36P246L (CCG→CTG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/16);  minor base G (0/1);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

ACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAAC  >  NC_000913_3_pae_pgi/4185946‑4186231
                                                                                                                                            |                                                                                                                                                 
aCAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCCGacc                                                                                                                                           >  2:533733/1‑147 (MQ=255)
 cAGCTGCGTGAAGAGCTGAACGAAACCAACGACGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCTTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCt                                                                                                                                          <  2:360643/149‑1 (MQ=255)
         tGAAGAGCTGAACGAAACCAACTCCGAAACCAGGCGAAAACAGCTGACCAAACGTATCAAACTGCTGGAAGCGTTCGTTCAGACTAGTAACAAACCAGAGTGGATGATCCTAACCTTTCTGCAGTAGCTGCTGCCAGGTCTGCGTCGGc                                                                                                                                  >  2:292166/1‑149 (MQ=255)
             gagCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGt                                                                                                                              <  1:84266/149‑1 (MQ=255)
                cTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTcc                                                                                                                           >  2:410365/1‑149 (MQ=255)
                  gAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGc                                                                                                                         >  1:178715/1‑149 (MQ=255)
                   aaCGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGGGTCCGCTGGTTCCGCt                                                                                                                        >  2:279975/1‑149 (MQ=255)
                              cTCAGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGTTGGTc                                                                                                             >  1:123923/1‑149 (MQ=255)
                               tCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGAt                                                                                                                    >  2:227342/1‑141 (MQ=255)
                               tCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGAt                                                                                                                    <  1:227342/141‑1 (MQ=255)
                                             tAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTg                                                                                              <  2:606404/149‑1 (MQ=255)
                                                 aaGCCGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCt                                                                                          <  2:461160/149‑1 (MQ=255)
                                                    ctgcgccaGCGTATCAAATTGCTGGATGCGGTCGTTCAGTTTGGTAACAACCCAGAGTGGATGATCCTGACCGTTCTGCCGCTACTGCTTCCAGATCTGCGTCCGCTGGTTCCGCTGTATGGTGGTCGTTTCGCGACTTCTGACCTGaa                                                                                       <  2:373158/142‑1 (MQ=255)
                                                     tGACCACGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAATCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAAc                                                                                      <  2:433895/149‑1 (MQ=255)
                                                       aCCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTCAACGa                                                                                    >  2:33732/1‑149 (MQ=255)
                                                        gcAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGAt                                                                                   <  2:450326/148‑1 (MQ=255)
                                                        ccAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGAt                                                                                   >  1:11721/1‑149 (MQ=255)
                                                          aaGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCt                                                                                 <  1:253645/149‑1 (MQ=255)
                                                                    aaCTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTcgc                                                                       <  2:178715/149‑1 (MQ=255)
                                                                          tgGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCAtt                                                                 >  1:345734/1‑149 (MQ=255)
                                                                          tgGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGGGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCAtt                                                                 >  2:68850/1‑149 (MQ=255)
                                                                             aaGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAAc                                                              >  1:534513/1‑149 (MQ=255)
                                                                                        cAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCg                                                   >  2:239256/1‑149 (MQ=255)
                                                                                             tGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGa                                              >  1:206784/1‑149 (MQ=255)
                                                                                                tAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAAc                                           <  2:55630/149‑1 (MQ=255)
                                                                                                  aCAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTcgc                                                                                                      <  2:110870/88‑1 (MQ=255)
                                                                                                  aCAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTcgc                                                                                                      >  1:110870/1‑88 (MQ=255)
                                                                                                      aCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTctgc                                     <  1:410365/149‑1 (MQ=255)
                                                                                                        cAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTctgctg                                   >  1:397966/1‑149 (MQ=255)
                                                                                                              ggctgatCCGGACCGTTCTGCCGGTAGTGCGGCCAACTCTGCGTCCGCTGCTTACGCTGGATGGTTGTCGTTTCGCGACTTCTGACCTGAACGATCTGAATCCTCGCGTCATTGACCGTAACAACCGTCTGAAACGTCTGCTGGATCTg                             <  2:594487/146‑1 (MQ=255)
                                                                                                               gatgatCCTGACCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTgg                            >  1:451507/1‑149 (MQ=255)
                                                                                                               gaggttCCGGACCGTTCTGCCGGTACTGCTTCCAGTTCGGCGTCCGCTGGTTCCGCGGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTATCAACCGTTTTAAGCGTCTGCTGGATCTgg                            <  2:598955/144‑1 (MQ=255)
                                                                                                                       ggACCGTTCTTGCGGTACTGCTACCAGATCTTCGGCCTCTCGTTCCGCTGCAGGCTGCTCTTGTCGCTACTTCTAACCTCAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACTTCTGCTGGATCTGGCTACGcct                    <  2:613894/148‑2 (MQ=255)
                                                                                                                         aCCGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATGTGTATCGTCGCGTCATTAACCATAACAAGCGTCTGAGACAGCTGCGGGATCTGGCTgcgcgcga                  >  1:605519/1‑145 (MQ=255)
                                                                                                                           cGTTCTGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGaag                >  1:122708/1‑147 (MQ=255)
                                                                                                                                tGCCGGTACTGCTGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACatcatc           <  1:51457/149‑1 (MQ=255)
                                                                                                                                         tgctgcCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAAc  >  1:383367/1‑149 (MQ=255)
                                                                                                                                         tgctgcCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAAc  >  2:385959/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
ACAGCTGCGTGAAGAGCTGAACGAAACCAACTCCGAAACCAAGCGTAAAAAGCTGACCAAGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAAC  >  NC_000913_3_pae_pgi/4185946‑4186231

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