Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,331,016 A→G Y310C (TAT→TGT)  acnA → aconitate hydratase AcnA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,331,0160AG76.5% 21.6 / 6.9 17Y310C (TAT→TGT) acnAaconitate hydratase AcnA
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (7/6);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

TGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGG  >  NZ_CP009273/1330929‑1331096
                                                                                       |                                                                                
tggGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAAtgtg                                                                                <  2:99145/90‑1 (MQ=255)
     aaaTTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATGTGGTGcc                                                                           <  2:389077/90‑1 (MQ=255)
                          gATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATc                                                      >  1:540658/1‑90 (MQ=255)
                           aTGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCg                                                     <  1:75967/90‑1 (MQ=255)
                             ggTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGAt                                                   <  2:635815/90‑1 (MQ=255)
                                       cACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAAccc                                         >  1:166944/1‑90 (MQ=255)
                                       cACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAAccc                                         >  1:482356/1‑90 (MQ=255)
                                        aCTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCt                                        >  1:124668/1‑90 (MQ=255)
                                              gTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTa                                  >  1:275754/1‑90 (MQ=255)
                                                                 aTTGCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGgcgc               <  2:144524/90‑1 (MQ=255)
                                                                 aTTGCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGgcgc               <  2:436687/90‑1 (MQ=255)
                                                                 aTTGCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGgcgc               <  2:536707/90‑1 (MQ=255)
                                                                 aTTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGgcgc               <  2:275754/90‑1 (MQ=255)
                                                                    gCCAATATGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGc            >  2:368516/1‑90 (MQ=255)
                                                                        atatGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAg        >  1:449931/1‑90 (MQ=255)
                                                                         tatGTCGCCAGAATGTGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGa       >  1:831438/1‑90 (MQ=255)
                                                                            gTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCa    >  2:235562/1‑90 (MQ=255)
                                                                              cGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAgg  >  1:559204/1‑90 (MQ=255)
                                                                                       |                                                                                
TGGGGAAATTCGTCGAATTTTATGGTGATGGTCTGGATTCACTACCGTTGGCGGATCGCGCCACCATTGCCAATATGTCGCCAGAATATGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGG  >  NZ_CP009273/1330929‑1331096

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