Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,170,226 A→T 27.1% intergenic (‑290/‑97) katG ← / → PP_3669 catalase‑peroxidase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,170,2260AT27.1% 34.7 / 11.0 22intergenic (‑290/‑97)katG/PP_3669catalase‑peroxidase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (8/8);  new base T (2/4);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGA  >  NC_002947/4170078‑4170375
                                                                                                                                                    |                                                                                                                                                     
cGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTTTATCGGCCTCATCGCGGTCaa                                                                                                                                                      <  1:259863/150‑1 (MQ=255)
      tGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCcgc                                                                                                                                                >  2:250892/1‑150 (MQ=255)
         tGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGAt                                                                                                                                             >  2:386409/1‑150 (MQ=255)
           aCGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGa                                                                                                                                           <  1:491985/150‑1 (MQ=255)
                 tGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCgg                                                                                                                                     >  2:486889/1‑150 (MQ=255)
                                     ggggCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggcacaa                                                                                                                 <  1:386409/150‑1 (MQ=255)
                                      gggCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATacaaggcacaag                                                                                                                <  1:250892/150‑1 (MQ=255)
                                            cGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATg                                                                                                                                            >  1:167910/1‑116 (MQ=255)
                                            cGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATg                                                                                                                                            <  2:167910/116‑1 (MQ=255)
                                                       aGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTc                                                                                               <  1:486889/150‑1 (MQ=255)
                                                                       agGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGTAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                           <  2:433636/90‑1 (MQ=255)
                                                                       agGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGTAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGa                                                                                                                                           >  1:433636/1‑90 (MQ=255)
                                                                            tAATAGAGTGGCTGTTAGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTg                                                                          >  2:246258/1‑150 (MQ=255)
                                                                             aaTAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                          <  1:215029/85‑1 (MQ=37)
                                                                             aaTAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACCGCGATGAg                                                                                                                                          >  2:215029/1‑85 (MQ=37)
                                                                                                  tGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAg                                                                     <  2:70581/133‑1 (MQ=255)
                                                                                                  tGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAg                                                                     >  1:70581/1‑133 (MQ=255)
                                                                                                               ctGGCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGAcc                                       <  1:246258/150‑1 (MQ=255)
                                                                                                                    aGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTAAGCGGGACGCCGAGCCCTCATGAACCGCAACGACCTgcgc                                  >  1:281804/1‑150 (MQ=255)
                                                                                                                                       ccTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGGCGACATGAACctgc               >  1:499687/1‑150 (MQ=255)
                                                                                                                                        cTCATCGCGGGCTTGACCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTg                                                                          <  2:220818/90‑1 (MQ=37)
                                                                                                                                         tCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgctg             <  2:222772/150‑1 (MQ=255)
                                                                                                                                           aTCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGt           <  2:226482/150‑1 (MQ=255)
                                                                                                                                                    aaGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGa  >  1:233502/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCGCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGCTGTTCGA  >  NC_002947/4170078‑4170375

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: