Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,809,612 C→T R497C (CGC→TGC)  katE → catalase HPII

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,809,6120CT78.8% 72.8 / 11.2 33R497C (CGC→TGC) katEcatalase HPII
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base T (13/13);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGT  >  NZ_CP009273/1809478‑1809743
                                                                                                                                      |                                                                                                                                   
ggATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAg                                                                                                                                 >  1:628700/1‑139 (MQ=255)
        acTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGt                                                                                                                         <  2:457325/139‑1 (MQ=255)
          tAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGttt                                                                                                                       >  1:532426/1‑139 (MQ=255)
          tAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGGTCGCGAGCGCAGCCCATCGttt                                                                                                                       >  1:331710/1‑139 (MQ=255)
            aCCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTgg                                                                                                                     <  1:659096/139‑1 (MQ=255)
             cccGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGc                                                                                                                    <  2:828150/139‑1 (MQ=255)
               cGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGa                                                                                                                  <  2:532426/139‑1 (MQ=255)
                 gCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGCGCGCGTGGACGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGaa                                                                                                                 >  2:725083/1‑138 (MQ=255)
                         cgaaccgaacTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTc                                                                                                         <  1:711227/138‑1 (MQ=255)
                                                 cTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGt                                                                                <  2:91214/139‑1 (MQ=255)
                                                   ggCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGt                                                                                >  2:190942/1‑137 (MQ=255)
                                                      cgcgcgAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGAc                                                                           <  2:628700/139‑1 (MQ=255)
                                                          cgAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGc                                                                                                                    >  1:462880/1‑94 (MQ=255)
                                                          cgAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGc                                                                                                                    <  2:462880/94‑1 (MQ=255)
                                                                acCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGag                                                                 <  1:719059/139‑1 (MQ=255)
                                                                  cgccgGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGagca                                                               <  2:577047/139‑1 (MQ=255)
                                                                   gccgGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGagcag                                                              <  2:125547/139‑1 (MQ=255)
                                                                                   ggcggTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGcc                                                          <  1:225776/127‑1 (MQ=255)
                                                                                   ggcggTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGcc                                                          >  2:225776/1‑127 (MQ=255)
                                                                                   ggcggTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATg                                              >  2:432551/1‑139 (MQ=255)
                                                                                                 aTACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAg                                <  1:190942/139‑1 (MQ=255)
                                                                                                 aTACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGtt                                     <  1:657497/134‑1 (MQ=255)
                                                                                                 aTACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGtt                                     >  2:657497/1‑134 (MQ=255)
                                                                                                  tACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGt                               >  1:692/1‑139 (MQ=255)
                                                                                                       ggAGCGCGTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAg                           >  2:830260/1‑138 (MQ=255)
                                                                                                           cgcgTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAg                      >  2:675925/1‑139 (MQ=255)
                                                                                                           cgcgTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAg                      >  2:171548/1‑139 (MQ=255)
                                                                                                             cgTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGAc                                                                           >  2:373178/1‑84 (MQ=255)
                                                                                                             cgTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGAc                                                                           <  1:373178/84‑1 (MQ=255)
                                                                                                             cgTGGAAGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGtc                    >  1:553519/1‑139 (MQ=255)
                                                                                                             cgTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGtc                    >  2:215460/1‑139 (MQ=255)
                                                                                                                   aGGCAATAAAGTTCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGtcgttcgt              <  1:171548/139‑1 (MQ=255)
                                                                                                                               tCGCGAGTGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGt  >  2:288194/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGACGCCATTTGAGCAGCGCCATATTGTCGATGGTTTCAGTTTTGAGTTAAGCAAAGTCGTTCGTCCGTATATTCGT  >  NZ_CP009273/1809478‑1809743

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑