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,7210AG98.2% 172.6 / ‑4.1 54E49G (GAA→GGA) agaVPTS N‑acetylgalactosamine transporter subunit IIB
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (22/31);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

TTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAG  >  NZ_CP009273/3273634‑3273806
                                                                                       |                                                                                     
ttCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGgaa                                                                                     <  2:1529671/90‑1 (MQ=255)
ttCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGgaa                                                                                     <  2:1556973/90‑1 (MQ=255)
 tCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAAt                                                                                    <  1:101666/90‑1 (MQ=255)
 tCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAAt                                                                                    <  2:600810/90‑1 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   >  1:1164595/1‑90 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   <  1:997412/90‑1 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   <  1:713011/90‑1 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   <  1:471854/90‑1 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   >  2:791014/1‑90 (MQ=255)
  cAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATg                                                                                   <  1:1466363/90‑1 (MQ=255)
      gggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAc                                                                               <  2:1152989/90‑1 (MQ=255)
      gggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTAc                                                                               <  1:210938/90‑1 (MQ=255)
      gggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTAATGGGAATGGTAc                                                                               <  1:1494419/90‑1 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  2:131683/1‑90 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              <  2:218281/90‑1 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  2:676112/1‑90 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  1:355165/1‑90 (MQ=255)
       ggTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  1:409716/1‑90 (MQ=255)
            gATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAg                                                                         <  2:601635/90‑1 (MQ=255)
                 gCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAggg                                                                    >  2:738555/1‑90 (MQ=255)
                    ggggCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATc                                                                 <  1:243973/90‑1 (MQ=255)
                    ggggCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATc                                                                 <  1:1145616/90‑1 (MQ=255)
                     gggCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCg                                                                >  1:195907/1‑90 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  1:1360068/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  2:808210/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  1:330174/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  1:1229268/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  2:355165/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  2:334641/90‑1 (MQ=255)
                           aTCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  1:1507169/90‑1 (MQ=255)
                           aTCTGGGGCTGGTAGCCAACGATGAGGTTGCCGAAAGTCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTAc                                                          <  2:754954/90‑1 (MQ=255)
                            tCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACg                                                         <  2:809/90‑1 (MQ=255)
                                  gCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTg                                                   >  2:583750/1‑90 (MQ=255)
                                    tGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              <  2:1029251/61‑1 (MQ=255)
                                    tGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACt                                                                              >  1:1029251/1‑61 (MQ=255)
                                     ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAc                                                >  2:961996/1‑90 (MQ=255)
                                     ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAc                                                >  2:926468/1‑90 (MQ=255)
                                     ggTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGAc                                                >  1:947801/1‑90 (MQ=255)
                                      gTAGCCAACGATGAGGTTGTCGAAGATCCGGTACAACAAAACCTGATGGGAAAGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                               <  1:528335/90‑1 (MQ=255)
                                      gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                               <  2:197689/90‑1 (MQ=255)
                                      gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                               <  2:245303/90‑1 (MQ=255)
                                      gTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACg                                               <  1:647974/90‑1 (MQ=255)
                                       tAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGc                                              <  1:375732/90‑1 (MQ=255)
                                         gCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTg                                            >  1:1226235/1‑90 (MQ=255)
                                                       tGCCGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGa                              >  1:698318/1‑90 (MQ=255)
                                                         ccGAAGATCCGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGaca                            <  2:332063/90‑1 (MQ=255)
                                                                  cGGTACAACAAAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTTATc                   >  1:76965/1‑90 (MQ=255)
                                                                        aacaaAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGcgccg             >  1:387651/1‑90 (MQ=255)
                                                                         acaaAACCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGcgccgc            >  2:1563049/1‑90 (MQ=255)
                                                                              aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc       >  1:1552970/1‑90 (MQ=255)
                                                                              aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc       <  2:1244452/90‑1 (MQ=255)
                                                                              aCCTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc       >  1:1013607/1‑90 (MQ=255)
                                                                              aCCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATc       >  2:10671/1‑90 (MQ=255)
                                                                                cTGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGa     <  1:738555/90‑1 (MQ=255)
                                                                                 tGATGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGAc    >  2:1146515/1‑90 (MQ=255)
                                                                                   aTGGGAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAg  >  1:780117/1‑90 (MQ=255)
                                                                                       |                                                                                     
TTCAATGGGTCGGATTTGCGGGGGCAAATCTGGTGCTGGTAGCCAACGATGAGGTTGCCGAAGATCCGGTACAACAAAACCTGATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAG  >  NZ_CP009273/3273634‑3273806

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