Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_017635 2,860,853 C→A R442S (AGG→AGT qseE ← two component system sensor histidine kinase QseE/GlrK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0176352,860,8530CA96.8% 113.5 / NA 31R442S (AGG→AGTqseEtwo component system sensor histidine kinase QseE/GlrK
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (18/12);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCC  >  NC_017635/2860806‑2860894
                                               |                                         
ctctCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACt                                         <  2:352231/50‑1 (MQ=255)
ctctCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACt                                         <  1:1240905/50‑1 (MQ=255)
  atcgtcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTg                                        <  1:1207884/48‑1 (MQ=255)
    cgtcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGc                                      >  2:1051176/1‑49 (MQ=255)
     gtcgACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGc                                      <  1:795586/48‑1 (MQ=255)
        gACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATg                                  <  2:182942/49‑1 (MQ=255)
        gACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATg                                  <  1:1506563/49‑1 (MQ=255)
         aCCAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCt                                >  2:1404348/1‑50 (MQ=255)
          ccAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCt                                >  1:2035649/1‑49 (MQ=255)
           cAGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTa                              <  1:875416/50‑1 (MQ=255)
            aGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTaa                             >  1:1172961/1‑50 (MQ=255)
            aGATACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTaa                             >  2:2039489/1‑50 (MQ=255)
               tACAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATcc                          <  1:1149090/50‑1 (MQ=255)
                 cAGTTCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATccc                         <  2:1961698/49‑1 (MQ=255)
                     tCCCCTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGAcc                    >  1:817530/1‑50 (MQ=255)
                      ccccTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCg                   >  1:642684/1‑50 (MQ=255)
                        ccTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCt                 <  1:1485786/50‑1 (MQ=255)
                        ccTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCt                 >  1:106749/1‑50 (MQ=255)
                         cTTGCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTg                >  1:1816308/1‑50 (MQ=255)
                          ttGCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGc               >  1:884006/1‑50 (MQ=255)
                            gCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGccc             >  2:1442823/1‑50 (MQ=255)
                            gCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGccc             >  1:1208717/1‑50 (MQ=255)
                            gCATACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGccc             >  1:121866/1‑50 (MQ=255)
                              aTACGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCtt           <  1:1178461/50‑1 (MQ=255)
                                aCGGCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCa         >  2:1056971/1‑50 (MQ=255)
                                  ggCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCAcc       <  2:1964380/50‑1 (MQ=255)
                                  ggCGAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCAcc       >  1:1235430/1‑50 (MQ=255)
                                     gAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcc    >  1:546541/1‑50 (MQ=255)
                                     gAATACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcc    >  2:222211/1‑50 (MQ=255)
                                      aaTACAATCACTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGccc   >  2:48908/1‑50 (MQ=255)
                                       aTACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGcccc  <  2:790650/50‑1 (MQ=255)
                                               |                                         
CTCTCGTCGACCAGATACAGTTCCCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCC  >  NC_017635/2860806‑2860894

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