Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,670,359 T→C G454G (GGA→GGG yphE ← sugar ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,670,3590TC97.9% 148.7 / ‑3.7 47G454G (GGA→GGGyphEsugar ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (31/15);  total (32/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAAC  >  NZ_CP009273/2670283‑2670443
                                                                            |                                                                                    
cgtgcTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         <  1:299223/90‑1 (MQ=255)
  tgcTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagccagc                                                                       <  2:171245/90‑1 (MQ=255)
          taaCAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGAc                                                               <  1:191910/90‑1 (MQ=255)
           aaCAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACa                                                              >  2:488513/1‑90 (MQ=255)
            aCAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACaa                                                             <  2:1338812/90‑1 (MQ=255)
             cAGGATGCGGTCACACACCAGCGGCGACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAAt                                                            <  1:1424061/90‑1 (MQ=255)
              aGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATa                                                           >  1:9731/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  2:200360/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         >  1:138703/1‑74 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         <  2:138703/74‑1 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  2:657214/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         <  2:339425/90‑1 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  2:1206538/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  2:1027500/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:1087651/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:903828/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:830694/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:560811/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:430071/1‑90 (MQ=255)
                gATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACg                                                         >  1:138138/1‑90 (MQ=255)
                  tGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGt                                                       <  1:657214/90‑1 (MQ=255)
                       tCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATc                                                  <  2:16750/90‑1 (MQ=255)
                         acacacCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATctg                                                >  1:1305914/1‑90 (MQ=255)
                              ccAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCCCGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGtt                                           <  1:1328240/90‑1 (MQ=255)
                              ccAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGtt                                           >  2:333988/1‑90 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAAt                                                            <  1:200016/71‑1 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAAt                                                            >  2:200016/1‑71 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg                                         >  1:425710/1‑90 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg                                         >  1:1352262/1‑90 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg                                         >  2:75466/1‑90 (MQ=255)
                                aGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTg                                         >  1:48643/1‑90 (MQ=255)
                                    gCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCtt                                     <  2:1244808/90‑1 (MQ=255)
                                    gCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCtt                                     >  1:521195/1‑90 (MQ=255)
                                      aaCTCCTCCACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCg                                   >  2:918350/1‑90 (MQ=255)
                                          cctccACTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGatat                               >  2:1213033/1‑90 (MQ=255)
                                               aCTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         >  1:287290/1‑43 (MQ=255)
                                               aCTTCACTGGAGATAAACACCACGCTTTTCCCTTCGGcagcca                                                                         <  2:287290/43‑1 (MQ=255)
                                                  tCACTGGAGATAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGAcgccg                       <  2:152946/90‑1 (MQ=255)
                                                          gaTAAACACCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTgg               >  2:745344/1‑90 (MQ=255)
                                                               acacCACGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCtcgt          >  1:1219895/1‑90 (MQ=255)
                                                                 accacGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCtcgtcg        >  2:1445683/1‑90 (MQ=255)
                                                                 accacGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCtcgtcg        >  2:626532/1‑90 (MQ=255)
                                                                    acGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGc     >  2:1376901/1‑90 (MQ=255)
                                                                    acGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGc     >  1:882047/1‑90 (MQ=255)
                                                                    acGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCATCGATATCGACGCCGCGCGTTGGCTCGTCGAGc     >  2:1195516/1‑90 (MQ=255)
                                                                     cGCTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGca    >  2:1297124/1‑90 (MQ=255)
                                                                       cTTTTCCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGcaac  <  2:1305914/90‑1 (MQ=255)
                                                                            |                                                                                    
CGTGCTGTAATAACAGGATGCGGTCACACACCAGCGGCAACTCCTCCACTTCACTGGAGATAAACACCACGCTTTTTCCTTCGGCAGCCAGCTCACGGACAATACGGTAAATCTGCTGTTTGGCTTCGATATCGACGCCGCGCGTTGGCTCGTCGAGCAAC  >  NZ_CP009273/2670283‑2670443

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