Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,798,416 (T)8→7 coding (61/936 nt) rfaS ← LPS core biosynthesis protein RfaS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,798,4090T.80.0% 17.3 / 3.4 10coding (68/936 nt)rfaSLPS core biosynthesis protein RfaS
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (3/5);  total (5/5)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

CCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACAT  >  NZ_CP009273/3798320‑3798498
                                                                                         |                                                                                         
ccAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACt                                                                                           <  1:662968/90‑1 (MQ=255)
 cAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtt                                                                                          <  1:833620/90‑1 (MQ=255)
        gCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtttttttca                                                                                   >  1:742534/1‑88 (MQ=255)
            tATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGt                                                                               >  1:320222/1‑90 (MQ=255)
                                 tttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAAt                                                                        <  1:388356/75‑1 (MQ=255)
                                 tttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAAt                                                                        >  2:388356/1‑75 (MQ=255)
                                  ttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATaa                                                        <  2:621206/90‑1 (MQ=255)
                                                           tttGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaaataa                                >  1:842019/1‑90 (MQ=255)
                                                                tgttgtTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTc                          <  1:709312/90‑1 (MQ=255)
                                                                   tgttATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCata                       >  1:888308/1‑90 (MQ=255)
                                                                     ttATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCataat                     <  2:742534/90‑1 (MQ=255)
                                                                              ttATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagttagt            >  2:529488/1‑90 (MQ=255)
                                                                                aTCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTc          <  2:75527/90‑1 (MQ=255)
                                                                                         ttttttttCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACAt  <  2:965463/90‑1 (MQ=255)
                                                                                         |                                                                                         
CCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACAT  >  NZ_CP009273/3798320‑3798498

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