Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,364,611 A→T 11.8% V24D (GTC→GAC)  folP ← dihydropteroate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,364,6110AT11.8% 145.1 / 11.3 51V24D (GTC→GAC) folPdihydropteroate synthase
Reads supporting (aligned to +/- strand):  ref base A (20/25);  new base T (4/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

GCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGACGCTCATCCTGGGCCGCAGGCGTGCCCGAG  >  NC_002947/5364462‑5364759
                                                                                                                                                     |                                                                                                                                                    
gcgcACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGa                                                                                                                                                      <  1:1120745/150‑1 (MQ=255)
     cccGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATgg                                                                                                                                                 <  1:986511/150‑1 (MQ=255)
       cGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGTCATGGGt                                                                                                                                               >  2:75880/1‑150 (MQ=255)
             cgTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGTCATGGGTAttctc                                                                                                                                         <  1:75880/150‑5 (MQ=255)
              gTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGTCATGGGTAttctca                                                                                                                                        >  1:464295/1‑145 (MQ=255)
                ggACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACaa                                                                                                                                      >  1:29577/1‑150 (MQ=255)
                     cGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGCGAATACCCATGACATGGGTATGCGACAAATCaa                                                                                                                                 <  1:803359/150‑1 (MQ=255)
                     cGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCaa                                                                                                                                 <  1:545184/150‑1 (MQ=255)
                      gCCACAGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGCGACATTGATAATACCCATGACATGGCTATGCGACAAACCaaa                                                                                                                                >  2:1202490/1‑149 (MQ=255)
                         aCCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAAc                                                                                                                             >  2:978779/1‑150 (MQ=255)
                             atatCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCgg                                                                                                                         <  2:1205567/150‑1 (MQ=255)
                             atatCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCgg                                                                                                                         >  2:149966/1‑150 (MQ=255)
                                   aTCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCg                                                                                                                   >  2:881099/1‑150 (MQ=255)
                                      aaCGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCAACAAATCAGTAACCCGGCTGCCGCaa                                                                                                                >  2:69840/1‑150 (MQ=255)
                                               ccGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCgg                                                                                                       <  1:881099/150‑1 (MQ=255)
                                                    ggcAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCgggtagg                                                                                                  <  1:729163/150‑1 (MQ=255)
                                                    ggcAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCgggtagg                                                                                                  >  2:729163/1‑150 (MQ=255)
                                                          cATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACAt                                                                                                                                                   <  2:270938/95‑1 (MQ=255)
                                                          cATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTg                                                                                            <  1:923076/150‑1 (MQ=255)
                                                          cATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTg                                                                                            <  2:29577/150‑1 (MQ=255)
                                                          cATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATAGTGAGAATACCCATGACAt                                                                                                                                                   >  1:270938/1‑95 (MQ=255)
                                                              gCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACt                                                                                        <  2:228781/150‑1 (MQ=255)
                                                                ttCGGCGTGCCGCAGGGCCTCGCCGCGCTGACTGAAGCGACCGCCATTGGATAAGCAATCGGGGGTGATATGGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACTCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGa                                                                                      <  1:183175/150‑1 (MQ=255)
                                                                tgcGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGa                                                                                      >  1:578096/3‑150 (MQ=255)
                                                                         gcgcAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGa                                                                             <  1:1185317/150‑1 (MQ=255)
                                                                             aGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGcc                                                                         <  1:149966/150‑1 (MQ=255)
                                                                              gggCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGccc                                                                        <  1:252740/150‑1 (MQ=255)
                                                                                 ccTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCa                                                                                                                                  <  1:1107745/89‑1 (MQ=255)
                                                                                 ccTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCa                                                                                                                                  >  2:1107745/1‑89 (MQ=255)
                                                                                      tcgCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGtt                                                                <  1:68272/150‑1 (MQ=255)
                                                                                             gactgaAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGc                                                         <  1:1067013/150‑1 (MQ=255)
                                                                                                         cGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGAt                                             <  1:686650/150‑1 (MQ=255)
                                                                                                         cGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGAt                                             >  2:393996/1‑150 (MQ=255)
                                                                                                              tCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTgg                                          <  2:1150836/148‑1 (MQ=255)
                                                                                                              tCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACAGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTgg                                          >  1:1150836/1‑148 (MQ=255)
                                                                                                                 gagaAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCAc                                                                                         <  1:1116590/98‑1 (MQ=255)
                                                                                                                 gagaAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCAc                                                                                         >  2:1116590/1‑98 (MQ=255)
                                                                                                                      ggAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGa                                >  1:50243/1‑150 (MQ=255)
                                                                                                                       gAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGAc                               >  2:852077/1‑150 (MQ=255)
                                                                                                                         aTCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGACGc                             <  1:393996/150‑1 (MQ=255)
                                                                                                                              gggTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGACGCTCATc                        <  2:578096/150‑1 (MQ=255)
                                                                                                                                     atTGAGAATACCCATGACATGGGTATGCGACAAATCaa                                                                                                                                 >  2:1157037/1‑38 (MQ=255)
                                                                                                                                     atTGAGAATACCCATGACATGGGTATGCGACAAATCaa                                                                                                                                 <  1:1157037/38‑1 (MQ=255)
                                                                                                                                      ttGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGACGCTCATCCTGGGCCg                >  1:1101249/1‑150 (MQ=255)
                                                                                                                                        tcgcaTACCCATGTCATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTgg                                          >  1:1018560/5‑122 (MQ=255)
                                                                                                                                        tcgcaTACCCATGTCATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTgg                                          <  2:1018560/118‑1 (MQ=255)
                                                                                                                                          gcaTACCCATGTCATGGGTATGCGACAAAGCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTCGTGTTCAGCCGCTGGGCGGCCGATTGGCGATTCCGGACGCTAATCCTGGGCCGCAgg            >  1:324110/3‑150 (MQ=255)
                                                                                                                                                ccATGACAGGGGGATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGACGCTCATCCTGGGCCGCAGGCGTGcc      >  2:599403/1‑150 (MQ=255)
                                                                                                                                                 cATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAg                                                                    <  1:343704/87‑1 (MQ=255)
                                                                                                                                                 cATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAg                                                                    >  2:343704/1‑87 (MQ=255)
                                                                                                                                                   tGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGACGCTCATCCTGGGCCGCAGGCGTGCCCGa   <  1:1185772/150‑1 (MQ=255)
                                                                                                                                                    gACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGACGCTCATCCTGGGCCGCAGGCGTGCCCGAg  >  2:705304/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACCATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGGGTATGCGACAAATCAAGAACCCGGCTGCCGCAAGGCAACCGGGTAGGGTACTGCACTGAGCTCATAGATGCCCTTAGTGTTCAGCCGCTGGGCCGCCGATTGGCGATTCCGGACGCTCATCCTGGGCCGCAGGCGTGCCCGAG  >  NC_002947/5364462‑5364759

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