Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,786,523 C→T 10.5% R197H (CGC→CAC)  ppsA ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,786,5230CT10.5% 106.2 / 3.8 38R197H (CGC→CAC) ppsAphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base C (22/12);  new base T (2/2);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

ATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACGTTGAG  >  NC_000913/1786381‑1786671
                                                                                                                                              |                                                                                                                                                    
aTGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCa                                                                                                                                                >  1:154475/1‑149 (MQ=255)
      gTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGtt                                                                                                                                          >  2:144267/1‑149 (MQ=255)
         aaaCTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGaa                                                                                                                                       <  2:65440/149‑1 (MQ=255)
                         aCCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACACCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAgcg                                                                                                                        <  3:21063/148‑1 (MQ=255)
                                  tGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCacacaccg                                                                                                              >  6:49782/1‑145 (MQ=255)
                                       gACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGCGCAGAGCGCCCCCTCCGTTGCACACCGGCGGAGACCGCCACACCACgggtct                                                                                                         >  6:135970/1‑145 (MQ=255)
                                           aTCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACcgg                                                                                                                                 >  1:167449/1‑121 (MQ=255)
                                           aTCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACcgg                                                                                                                                 <  2:167449/121‑1 (MQ=255)
                                                aCCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAAcc                                                                                                  <  2:154475/147‑1 (MQ=255)
                                                 ccAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTg                                                                                               <  1:144267/149‑1 (MQ=255)
                                                         ccATTCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATAACGCCAGATGATGCTAGTTCAGAGCGCACCATCCGTTGAACACTGGCGGAGAGCGCCACACCCCGGTGATCGTCACACTGGTGCacac                                                                                       >  2:94596/1‑149 (MQ=255)
                                                              cGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATaa                                                                                  <  4:97520/149‑1 (MQ=255)
                                                                gAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATaa                                                                                  >  2:102597/1‑147 (MQ=255)
                                                                gAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAaga                                                                                <  1:57274/149‑1 (MQ=255)
                                                                 aaGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGTGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATaa                                                                                  >  8:31758/1‑146 (MQ=255)
                                                                       tAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGgcg                                                                         >  8:21204/1‑149 (MQ=255)
                                                                             ccaccTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGACAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCAGGTGCACACGATAAGAGATGGCGCGATCg                                                                   >  7:91545/1‑149 (MQ=255)
                                                                              caccTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGTGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAAGCCTGGTGCACACGATAaga                                                                                <  4:163177/135‑1 (MQ=255)
                                                                              caccTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGTGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAaga                                                                                >  3:163177/1‑135 (MQ=255)
                                                                              caccTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGt                                                                  >  4:127017/1‑149 (MQ=255)
                                                                                    gTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGTGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAAc                                                             <  7:31758/148‑1 (MQ=255)
                                                                                     tCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTg                                                                                                           >  5:197832/1‑101 (MQ=255)
                                                                                     tCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTg                                                                                                           <  6:197832/101‑1 (MQ=255)
                                                                                     tCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTg                                                                                                           <  6:194391/101‑1 (MQ=255)
                                                                                     tCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTg                                                                                                           >  5:194391/1‑101 (MQ=255)
                                                                                                   ggTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACAt                                             >  7:47436/1‑149 (MQ=255)
                                                                                                   ggTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACAACGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACAt                                             <  5:49782/149‑1 (MQ=255)
                                                                                                             gagaACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGc                                   >  6:19853/1‑149 (MQ=255)
                                                                                                                  cATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCAc                                >  3:101972/1‑147 (MQ=255)
                                                                                                                   aTCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCAc                                >  6:71305/1‑146 (MQ=255)
                                                                                                                       cGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACg                         <  5:180778/149‑1 (MQ=255)
                                                                                                                          cAGATGATGCGAGGTCAGAGCGCACAATCCGTTGAACACCGGCTGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGAGCGATAGTTAAACAGAGAAGCAAATACATGTATCACTGCCACGAGAACGGCg                      >  6:165576/1‑149 (MQ=255)
                                                                                                                            gatgagGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCATGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGt                     >  3:136668/1‑148 (MQ=255)
                                                                                                                                       gTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGaa         >  4:109693/1‑149 (MQ=255)
                                                                                                                                        tCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAAc        <  5:19853/149‑1 (MQ=255)
                                                                                                                                        tCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAAc        >  4:1692/1‑149 (MQ=255)
                                                                                                                                         cAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACg       >  2:130324/1‑149 (MQ=255)
                                                                                                                                            agCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATATGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACGtt     >  1:3264/1‑148 (MQ=255)
                                                                                                                                             gcgcACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTTGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACGTTGa   <  3:127017/149‑1 (MQ=255)
                                                                                                                                              cgcACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACGTTGAg  <  7:21204/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ATGCACGTAAAACTCATCCGGGTTAACCGCACCCTGCACGACCATCTCACCAAGGCCCCATGCGGAAGTGATAAACACCACCTGGTCAAAGCCGGATTCGGTATCAATGGAGAACATCACGCCAGATGATGCGAGGTCAGAGCGCACCATCCGTTGAACACCGGCGGAGAGCGCCACACCACGGTGATCGTAACCCTGGTGCACACGATAAGAGATGGCGCGATCGTTAAACAGAGAAGCAAATACATGTTTCACTGCCACGAGAACGGCGTCAAAACCCTGAACGTTGAG  >  NC_000913/1786381‑1786671

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: