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.90.0% 24.5 / 0.6 11coding (68/936 nt)rfaSLPS core biosynthesis protein RfaS
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (5/4);  total (5/6)
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.01e-01

GCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCC  >  NZ_CP009273/3798328‑3798491
                                                                                 |                                                                                  
gCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtttttttca                                                                            <  1:114727/90‑3 (MQ=255)
gCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtttttttca                                                                            >  2:736399/1‑88 (MQ=255)
   ttATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGt                                                                        >  2:163774/1‑90 (MQ=255)
     atccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTAc                                                                      >  2:414312/1‑90 (MQ=255)
          tccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAAt                                                                 >  1:637665/1‑90 (MQ=255)
                                            atatTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAAtttctg                               <  1:37059/90‑4 (MQ=255)
                                              atTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaa                             >  2:206643/1‑90 (MQ=255)
                                                     tGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAaataaata                      <  1:157894/90‑1 (MQ=255)
                                                        tgttgtTATATTTTTTATCTTGTCCTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTc                   <  1:248931/90‑1 (MQ=255)
                                                         gttgttATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTc                   <  2:345948/90‑1 (MQ=255)
                                                                    ttttATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagtta       <  1:680654/90‑1 (MQ=255)
                                                                      ttATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagttagt     >  2:43722/1‑90 (MQ=255)
                                                                         tCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTcc  <  2:637665/90‑1 (MQ=255)
                                                                                 |                                                                                  
GCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCC  >  NZ_CP009273/3798328‑3798491

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