Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,573,157 A→T I148N (ATC→AAC)  yddA ← ABC transporter ATP‑binding protein/permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,573,1570AT93.3% 47.7 / ‑3.9 16I148N (ATC→AAC) yddAABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (5/9);  minor base C (0/1);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

AATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAAGTAGT  >  NZ_CP009273/1573087‑1573238
                                                                      |                                                                                 
aaTAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTc                                                                <  1:42969/90‑1 (MQ=255)
         aaaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACg                                                       >  2:553793/1‑90 (MQ=255)
          aaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                      <  1:89894/90‑1 (MQ=255)
          aaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                      <  2:576017/90‑1 (MQ=255)
              tGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGt                                                  <  2:423707/90‑1 (MQ=255)
               gATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTc                                                                                 <  1:362571/58‑1 (MQ=255)
               gATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTc                                                 >  1:134934/1‑90 (MQ=255)
                      aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTa                                          >  1:222703/1‑90 (MQ=255)
                      aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTa                                          >  1:475807/1‑90 (MQ=255)
                      aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTa                                          >  2:340434/1‑90 (MQ=255)
                        gCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGCTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATc                                        <  2:11680/90‑1 (MQ=255)
                                       gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  2:200631/90‑1 (MQ=255)
                                       gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  2:371129/90‑1 (MQ=255)
                                                      aCTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGt          <  2:519510/90‑1 (MQ=255)
                                                         gAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGaa       >  2:174451/1‑90 (MQ=255)
                                                              ttttGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAagtagt  <  2:428910/90‑1 (MQ=255)
                                                                      |                                                                                 
AATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAAGTAGT  >  NZ_CP009273/1573087‑1573238

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