Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,819,176 C→A 21.0% intergenic (‑241/‑15) mltB ← / → srlA membrane‑bound lytic murein transglycosylase B/glucitol/sorbitol‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,819,1760CA21.0% 39.4 / 4.7 24intergenic (‑241/‑15)mltB/srlAmembrane‑bound lytic murein transglycosylase B/glucitol/sorbitol‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (7/11);  new base A (2/3);  total (10/14)
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

AACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCATTCTTCCGCTGTTGATTAGCCTGCTGGTTATCATGAACGCACTGA  >  CP009273/2819039‑2819323
                                                                                                                                         |                                                                                                                                                   
aactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAGTATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAACCTGATTAGATTAGGTTGCCGTTTGGTAATAAATCAATAAATCATGAAg                                                                                                                                                <  1:241701/143‑1 (MQ=255)
 actcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGaa                                                                                                                                         <  2:286947/149‑1 (MQ=255)
      ttCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTATTAAAAAAATAAATCTTAAaagaaaga                                                                                                                                          >  1:367588/1‑143 (MQ=255)
          ccTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACaa                                                                                                                                      <  1:161770/143‑1 (MQ=255)
           cgcGCTATCCTTTGATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAAt                                                                                                                                     <  1:6624/141‑1 (MQ=255)
           cTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATATTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCATGAAGGAGAGaaacat                                                                                                                                     >  1:378434/1‑139 (MQ=255)
             cGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGaaa                                                                                                                             >  2:148717/1‑149 (MQ=255)
                tATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGCTTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAg                                                                                                                                >  1:33665/1‑143 (MQ=255)
                tATCCTTTTATTCAAACTTTCAAATTAAAATAGTTATCTTTCGTTTTGCGATCAAAATAACAATTTTATGTCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCa                                                                                                                          >  2:407753/1‑149 (MQ=255)
                     ttttATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCATGAAGGAGAGAACAATGATAGAAACCATTACt                                                                                                                     >  2:245683/1‑149 (MQ=255)
                     ttttATTCAAACTTTCAAATTAAAATATTTATATTTCATTTTGCGATCAAAATAAAACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATACTGAAGGAGAGAACAATGATAGAAAcc                                                                                                                           >  1:297645/1‑143 (MQ=255)
                          ttCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTAc                                                                                                                      <  1:391010/143‑1 (MQ=255)
                                             atatTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCCGATGGACAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCggg                                                                                             <  2:409858/149‑1 (MQ=255)
                                                  tagcTTTCATTTTGTGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCgg                                                                                              <  1:215461/140‑1 (MQ=255)
                                                                            cacTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCg                                                              <  2:327469/149‑1 (MQ=255)
                                                                              cTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCCAAAGGGCGGAGAGGTGTTTACCggg                                                            <  2:243098/149‑1 (MQ=255)
                                                                                                tGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCAttctt                                          <  2:23924/149‑1 (MQ=255)
                                                                                                         tagGTTGCCGTTTGGTAATAAAACAATAAATCACGAAGGAGACAACAAGGATAGAACCCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCATTCTTCCg                                       <  1:126463/143‑1 (MQ=255)
                                                                                                             ttGCGGTTTGGTAATAAAACAATAAATCCTTAAGGCGAGAACAATTATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCATTCTTCCGCTGTTGATTa                             <  2:332772/149‑1 (MQ=255)
                                                                                                                  gTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCATTCTTCCGCTGTTGATTAGCctg                        >  2:290772/1‑149 (MQ=255)
                                                                                                                    ttGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTGCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCAGTCTTCCGCTGTGGAGTAg                            >  1:356264/1‑143 (MQ=255)
                                                                                                                     tGGTAATAAAACAATAAATCATGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCATTCTTCCGCTGTTGATTAGc                           <  1:259665/143‑1 (MQ=255)
                                                                                                                        ttataaAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCATTCTTCCGCTGTTGGTTAGCctt                        >  1:338459/3‑142 (MQ=255)
                                                                                                                           taaAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCATTCTTCCGCTGTTGATTAGCCTGCTGGTTATc               <  2:73491/149‑1 (MQ=255)
                                                                                                                                        ccTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCATTCTTCCGCTGTTGATTAGCCTGCTGGTTATCATGAACGCACTGa  >  2:128382/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
AACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGCAGAGTGGTTTATCGGGCTGTTCCAAAAGGGCGGAGAGGTGTTTACCGGGATGGTGACCGGCATTCTTCCGCTGTTGATTAGCCTGCTGGTTATCATGAACGCACTGA  >  CP009273/2819039‑2819323

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