Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,447,483 G→T intergenic (‑10/+25) aroC ← / ← prmB chorismate synthase/50S ribosomal subunit protein L3 N(5)‑glutamine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,447,4830GT96.9% 86.7 / ‑6.0 33intergenic (‑10/+25)aroC/prmBchorismate synthase/50S ribosomal subunit protein L3 N(5)‑glutamine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (18/13);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAAC  >  NC_000913/2447343‑2447615
                                                                                                                                            |                                                                                                                                    
tCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCttt                                                                                                                                    >  2:189104/1‑143 (MQ=255)
                    cGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGttt                                                                                                                <  2:479660/143‑1 (MQ=255)
                       cAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGCCCCATTGTGTTGCCTGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACt                                                                                                             <  1:676422/143‑1 (MQ=255)
                            gAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGACAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATAGTTGTGTTTGCGTGTTTACTTAATc                                                                                                        >  2:576915/1‑143 (MQ=255)
                            gAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATc                                                                                                        >  1:217351/1‑143 (MQ=255)
                            gAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGAGTATATAATAAACAGATCACCTTTATTGTTGTGTTTGCGTGTTTACTTAATc                                                                                                        >  2:378168/1‑143 (MQ=255)
                               tGCCTGGCCGAACACCAGCGACGATGCGGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCttt                                                                                                     <  2:91223/143‑1 (MQ=255)
                                 ccTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGCTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTCCTTAATCTTtat                                                                                                   >  2:356774/1‑143 (MQ=255)
                                     gCGGAACACCATCGACGATGCAGCCCAGCGCCAGCCCTTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAAt                                                                                               <  2:85053/143‑1 (MQ=255)
                                           cacCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGACTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCGCGGCTCCTTTATTGTGGTGTTTGCGTGTTCCGTTACTCATTTTAGATCGCGaa                                                                                         >  2:334794/1‑143 (MQ=255)
                                            acCATCGACGATGCAGCCGAGCGCCAGTCCGTGCGATTCGCCGTAGGTGGTTACGCGTAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCTTTTATTTTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGaaa                                                                                        <  2:703576/143‑1 (MQ=255)
                                              cATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTCCGCGAGAGAGTTGTCCAATTATGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGATTGCGTGTTTACTGAATCTTTATAATTCGCGAAATg                                                                                      >  1:240840/1‑143 (MQ=255)
                                                          aGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGGGTTTGCGTGTGTACTTAATCTTTATAAATCGCGGAGTGTTCTCGTGCGGc                                                                          >  1:464271/1‑143 (MQ=255)
                                                              gAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGCGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTCATCTTTCTAAATCGCGACATGTTCTCGCTCGGCAATa                                                                      >  1:16397/1‑143 (MQ=255)
                                                                    cAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGTTTTTTTTTTTTTGTGTTTGTGTGTTTACTTAAGCTTTATAAAGCGGGAAATGTTCTCGTGCGGCAATAAgctgc                                                                <  1:161889/143‑1 (MQ=255)
                                                                              taTTCGCCGAAGCTGGGTACGCGAAAGCGTTGTCCAAGTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGa                                                      <  2:442113/142‑1 (MQ=255)
                                                                                ttCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGc                                                    <  2:217351/143‑1 (MQ=255)
                                                                                  cGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCAt                                                  <  1:189104/143‑1 (MQ=255)
                                                                                    ccGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGCCTCGGGTATCGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAGATGTTGTCGTGCGGCAATAAGCTGCTCTTCGGCGAGcgcac                                                >  2:864477/1‑139 (MQ=255)
                                                                                     cGAAGGGGGTTACGCGAAAGAGTTGTGCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGGTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCGCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATaaa                                               <  2:804582/143‑1 (MQ=255)
                                                                                      gAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAac                                              >  1:839595/1‑143 (MQ=255)
                                                                                      gAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAac                                              >  1:838340/1‑143 (MQ=255)
                                                                                                   cgAAAGAGTTGTCCAATTGTGTTCCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATcgccgc                                 <  1:95249/143‑1 (MQ=255)
                                                                                                           ttGTCCAATTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATcc                         >  1:745043/1‑142 (MQ=255)
                                                                                                             gTCCAAGTGGGTTTCCAGCGATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCaaa                       <  2:240554/143‑1 (MQ=255)
                                                                                                               ccACTTGTGTTTCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGACATGTTCCCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACt                     <  1:773549/143‑1 (MQ=255)
                                                                                                                     gtgtTTCCAGCCATCACGCCTCCCTCATTGTTGTGTTTGCGTGTTTACTTAATCTTTATCAATCGCGCAATGTTCTCGTGCGTCAACAACCTGCTCTTTGGTGAGCACAAACACACCATCGCCGCCGCGATCAAACCccagcc               >  2:300119/1‑143 (MQ=255)
                                                                                                                        tttCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTCCTTAATCTTTATAAATCGCGAAATGTTCTCGTCCGGCAATAAGCTGCTCTTTGGTGAGCATAAACCCCCCATCGCCGCCGTTCTCACACTCCACCCAgg            >  1:723224/1‑143 (MQ=255)
                                                                                                                        tttCCAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTCACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATACACACACCATCGCCGCCGCTACCAGACTCCCGCCAgg            >  2:269451/1‑143 (MQ=255)
                                                                                                                            cAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTCTCAAACTCCAGCCAGGTGaa        >  1:697933/1‑143 (MQ=255)
                                                                                                                            cAGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGaa        >  1:497443/1‑143 (MQ=255)
                                                                                                                             aGCCATCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGGGAAc       >  1:594076/1‑143 (MQ=255)
                                                                                                                               ccATCACGGCTCCTTTATTGTTAGGTGGGCGTGGTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACgg     >  2:307398/1‑143 (MQ=255)
                                                                                                                                  tCACGGCTCCTTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATCTTCTCGTGGGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCGTGCAGGTGGACGGAAc  >  1:55772/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                    
TCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGGTGAGCATAAACACACCATCGCCGCCGTTATCAAACTCCAGCCAGGTGAACGGAAC  >  NC_000913/2447343‑2447615

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: