Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,819,095 T→C intergenic (‑160/‑96) mltB ← / → srlA lytic murein transglycosylase B/PTS glucitol/sorbitol transporter subunit IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,819,0950TC94.7% 53.4 / ‑3.7 19intergenic (‑160/‑96)mltB/srlAlytic murein transglycosylase B/PTS glucitol/sorbitol transporter subunit IIC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (10/8);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATA  >  NZ_CP009273/2819012‑2819171
                                                                                   |                                                                            
cTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCAtttt                                                                        >  2:235539/1‑90 (MQ=255)
    ccGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGa                                                                    <  1:613279/90‑1 (MQ=255)
      gATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATc                                                                  <  2:330866/90‑1 (MQ=255)
            gACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATa                                                            >  1:115905/1‑90 (MQ=255)
                ttACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAacac                                                        >  1:10031/1‑90 (MQ=255)
                          caactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCt                                              >  2:685436/1‑90 (MQ=255)
                            actcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCttt                                            >  1:713873/1‑90 (MQ=255)
                                 ttCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCGTTTTGCGATCAAAATAACACTTTTAAATCTTTCAATc                                       <  1:693124/90‑1 (MQ=255)
                                 ttCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATc                                       >  2:590009/1‑90 (MQ=255)
                                 ttCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATc                                       >  1:306814/1‑90 (MQ=255)
                                 ttCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATc                                       >  2:835131/1‑90 (MQ=255)
                                     ccTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTgat                                   <  1:426591/90‑1 (MQ=255)
                                      cTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTgatt                                  <  2:306814/90‑1 (MQ=255)
                                      cTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTgatt                                  >  1:377124/1‑90 (MQ=255)
                                                    aTTCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCg                    >  1:795494/1‑90 (MQ=255)
                                                      tCAAACTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGtt                  >  2:692800/1‑90 (MQ=255)
                                                        aaaCTTTCAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTg                <  1:647718/90‑1 (MQ=255)
                                                               cAAATTAAAATATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGtaataa         <  2:795494/90‑1 (MQ=255)
                                                                      aaaTATTTATCTTCCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATa  <  1:685436/90‑1 (MQ=255)
                                                                                   |                                                                            
CTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATA  >  NZ_CP009273/2819012‑2819171

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