Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,273,721 A→G E49G (GAA→GGA)  agaV → PTS N‑acetylgalactosamine transporter subunit IIB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,273,7210AG91.3% 58.3 / ‑0.2 23E49G (GAA→GGA) agaVPTS N‑acetylgalactosamine transporter subunit IIB
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (8/13);  total (9/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 = 9.04e-01

CAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGA  >  NZ_CP009273/3273636‑3273803
                                                                                     |                                                                                  
cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                <  2:596037/90‑1 (MQ=255)
cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                >  2:390473/1‑90 (MQ=255)
       tCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTgg                                                                         >  2:662977/1‑90 (MQ=255)
          gATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAg                                                                      >  2:691026/1‑90 (MQ=255)
               gCGGGGGCAAATCTGGTGCTGGTAG‑CAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGa                                                                <  2:244619/90‑1 (MQ=255)
                 gggggCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGAt                                                               <  1:405854/90‑1 (MQ=255)
                   gggCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCg                                                             >  2:486757/1‑90 (MQ=255)
                        aaTCTGTTGCT‑GTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                       <  1:512948/90‑1 (MQ=255)
                         aTCTGGTGGTGGTAGCCAACGATGAGGTTGCCGAAGGTCCGGTACAACAAAACCTTATTGGAATAGTTCTGGCAGAAGGGATCGCCGTAc                                                       <  2:263063/90‑1 (MQ=255)
                         aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                       <  1:76196/90‑1 (MQ=255)
                         aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                       <  2:705431/90‑1 (MQ=255)
                         aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                       <  1:586907/90‑1 (MQ=255)
                         aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                       <  2:199975/90‑1 (MQ=255)
                           cTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGt                                                     <  1:149956/90‑1 (MQ=255)
                              gTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTc                                                  >  1:712981/1‑90 (MQ=255)
                                  tGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGa                                              >  1:543626/1‑90 (MQ=255)
                                    gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                            <  1:702526/90‑1 (MQ=255)
                                    gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                            <  2:564843/90‑1 (MQ=255)
                                          aaCGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCaa                                      >  2:212374/1‑90 (MQ=255)
                                           aCGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCaaa                                     >  2:580403/1‑90 (MQ=255)
                                                      gCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGac                          >  2:496833/1‑90 (MQ=255)
                                                       ccGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGaca                         <  2:230658/90‑1 (MQ=255)
                                                                              cTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGa  <  2:712981/90‑1 (MQ=255)
                                                                                     |                                                                                  
CAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGA  >  NZ_CP009273/3273636‑3273803

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