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.95.2% 78.9 / ‑1.3 21coding (68/936 nt)rfaSLPS core biosynthesis protein RfaS
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (8/12);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

CAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCC  >  NZ_CP009273/3798321‑3798491
                                                                                        |                                                                                  
cAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtt                                                                                   >  2:1511386/1‑90 (MQ=255)
       gCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtttttttca                                                                            >  1:1786902/1‑88 (MQ=255)
       gCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtttttttca                                                                            >  1:891428/1‑88 (MQ=255)
        cATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCa                                                                           <  1:1276202/90‑1 (MQ=255)
            atccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTAc                                                                      >  2:3011260/1‑90 (MQ=255)
            atccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTAc                                                                      <  1:10518/90‑1 (MQ=255)
            atccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTAc                                                                      >  2:3022498/1‑90 (MQ=255)
                 tccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACt                                                                                    >  2:2381273/1‑72 (MQ=255)
                 tccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACt                                                                                    <  1:2381273/72‑1 (MQ=255)
                 tccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAAt                                                                 >  1:2256730/1‑90 (MQ=255)
                                      tattaATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGc                                            <  1:3022498/90‑1 (MQ=255)
                                          aaTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAAcc                                        <  2:2256730/90‑1 (MQ=255)
                                             ccAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCaa                                     <  2:788032/90‑1 (MQ=255)
                                               aGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAAtt                                   >  1:1145110/1‑90 (MQ=255)
                                                   atatTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTtata                               <  1:452190/90‑1 (MQ=255)
                                                     atTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaa                             <  2:151990/90‑1 (MQ=255)
                                                        tttttGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaaata                          >  1:1773073/1‑90 (MQ=255)
                                                          tttGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaaataaa                        >  2:1962617/1‑90 (MQ=255)
                                                          tttGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaaataaa                        >  1:2003962/1‑90 (MQ=255)
                                                            tGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAaataaata                      <  2:462326/90‑1 (MQ=255)
                                                              gtgtTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGt                    <  2:801527/90‑1 (MQ=255)
                                                               tgttgtTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTc                   <  2:1145110/90‑1 (MQ=255)
                                                                     tataTTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCataata             <  2:1773073/90‑1 (MQ=255)
                                                                      atatTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCAtaataa            <  2:317852/90‑1 (MQ=255)
                                                                          tttttATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagtt        >  2:1286089/1‑90 (MQ=255)
                                                                                tCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTcc  <  1:1962617/90‑1 (MQ=255)
                                                                                        |                                                                                  
CAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCC  >  NZ_CP009273/3798321‑3798491

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