Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,260,284 A→G D372G (GAC→GGC)  hemA → glutamyl‑tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,260,2840AG91.7% 64.8 / ‑0.3 24D372G (GAC→GGC) hemAglutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (9/13);  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 = 9.92e-01

AAACCAGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGC  >  NZ_CP009273/1260152‑1260418
                                                                                                                                    |                                                                                                                                      
aaaCCAGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCa                                                                                                                                  <  1:337959/139‑1 (MQ=255)
aaaCCAGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCa                                                                                                                                  <  1:449026/139‑1 (MQ=255)
     aGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCa                                                                                                                             <  2:422651/139‑1 (MQ=255)
     aGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGCTAGCGGCCCTTGAGCAGGGCGGCGGCGCgccagcca                                                                                                                             >  2:763681/1‑139 (MQ=255)
                    tGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTgg                                                                                                              >  1:842038/1‑139 (MQ=255)
                    tGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTgg                                                                                                              >  1:353852/1‑139 (MQ=255)
                       cTGCGGGCAACAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCTCAAGCCATTATGCAGGATCTGGCAt                                                                                                           <  2:240078/139‑1 (MQ=255)
                          cGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGGGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGa                                                                                                        <  1:80209/139‑1 (MQ=255)
                                                   ttCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATgc                                                                               >  1:784625/1‑139 (MQ=255)
                                                    tCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAg                                                                                                                                <  1:529676/89‑1 (MQ=255)
                                                    tCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAg                                                                                                                                >  2:529676/1‑89 (MQ=255)
                                                    tCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATgcg                                                                              >  2:148796/1‑139 (MQ=255)
                                                              cGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAAt                                                                    <  2:842038/139‑1 (MQ=255)
                                                              cGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAAt                                                                    <  2:784625/139‑1 (MQ=255)
                                                                     aGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTc                                                              >  2:266572/1‑138 (MQ=255)
                                                                           agCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGgc                                                       >  2:778356/1‑139 (MQ=255)
                                                                                 ttCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCg                                                 >  1:491789/1‑139 (MQ=255)
                                                                                   cgcgATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGt                                                <  2:695073/138‑1 (MQ=255)
                                                                                    gcgATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGa                                              <  1:182684/139‑1 (MQ=255)
                                                                                     cgATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGAc                                             >  1:498865/1‑139 (MQ=255)
                                                                                     cgATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGAc                                             <  2:498865/139‑1 (MQ=255)
                                                                                     cgATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGAc                                             <  2:353852/139‑1 (MQ=255)
                                                                                               aCCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATaacg                                   <  1:622799/139‑1 (MQ=255)
                                                                                                   ccAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGGCGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATaacgaacg                               <  1:266572/139‑1 (MQ=255)
                                                                                                                             ggcggcggcGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGa             <  2:530267/131‑1 (MQ=255)
                                                                                                                             ggcggcggcGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGa             >  1:530267/1‑131 (MQ=255)
                                                                                                                             ggcggcggcGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCg     >  1:688609/1‑139 (MQ=255)
                                                                                                                                ggcggcGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGc  <  2:844366/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
AAACCAGCGAATTTATGGCGTGGCTGCGAGCACAAAGCGCCAGCGAAACCATTCGCGAGTATCGCAGCCAGGCAGAGCAAGTTCGCGATGAGTTAACCGCCAAAGCGTTAGCGGCCCTTGAGCAGGGCGGCGACGCGCAAGCCATTATGCAGGATCTGGCATGGAAACTGACTAACCGCTTGATCCATGCGCCAACGAAATCACTTCAACAGGCCGCCCGTGACGGGGATAACGAACGCCTGAATATTCTGCGCGACAGCCTCGGGC  >  NZ_CP009273/1260152‑1260418

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