Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,394,663 T→C G456G (GGA→GGG)  abgT ← p‑aminobenzoyl‑glutamate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,394,6630TC93.3% 37.2 / ‑3.6 15G456G (GGA→GGG) abgTp‑aminobenzoyl‑glutamate transporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (10/4);  total (11/4)
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.18e-01

TAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGGCAATACGGATGAGTCGGCAATACGAAAGAGGATTTGCG  >  NZ_CP009273/1394601‑1394742
                                                              |                                                                               
tAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGAt                                                      >  1:539032/1‑90 (MQ=255)
tAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGAt                                                      >  2:42386/1‑90 (MQ=255)
  gggCAAGGCTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATAc                                                    >  1:90738/1‑90 (MQ=255)
     cAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCgg                                                 >  1:247499/1‑90 (MQ=255)
       aGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGgcg                                               >  2:268106/1‑90 (MQ=255)
          cTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTa                                            <  1:42386/90‑1 (MQ=255)
            aaCGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTaaa                                          >  1:416569/1‑90 (MQ=255)
            aaCGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTaaa                                          >  2:259394/1‑90 (MQ=255)
                  tAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGGCAAt                                    >  1:578184/1‑90 (MQ=255)
                              aGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGGCAATACGGATGAGTCg                        >  1:372261/1‑90 (MQ=255)
                                     cgTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGGCAATACGGATGAGTCGGCAATAc                 <  1:259394/90‑1 (MQ=255)
                                     cgTCTGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGGCAATACGGATGAGTCGGCAATAc                 <  1:333703/90‑1 (MQ=255)
                                         tGGTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGGCAATACGGATGAGTCGGCAATACGaaa             >  2:65236/1‑90 (MQ=255)
                                          ggTTTGTAGCGTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGGCAATACGGATGAGTCGGCAATACGAAag            <  1:201388/90‑1 (MQ=255)
                                                    gTTGCAGGAACCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGGCAATACGGATGAGTCGGCAATACGAAAGAGGATTTgcg  >  2:521673/1‑90 (MQ=255)
                                                              |                                                                               
TAGGGCAAGACTAACGAATAGTAAGTACCCAGTTTCGCGTCTGGTTTGTAGCGTTGCAGGAATCCAAGAAACAGTGGAACAAAAGGAGATACCGGCGCTAAAGGCAATACGGATGAGTCGGCAATACGAAAGAGGATTTGCG  >  NZ_CP009273/1394601‑1394742

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