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,1570AT95.5% 79.2 / ‑3.3 22I148N (ATC→AAC) yddAABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (8/13);  total (8/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 = 7.04e-01

CCCGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAA  >  NZ_CP009273/1573070‑1573233
                                                                                       |                                                                            
cccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTc                                                                            <  2:583731/90‑1 (MQ=255)
                 aaTAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTc                                                           <  2:238668/90‑1 (MQ=255)
                 aaTAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTc                                                           <  1:465871/90‑1 (MQ=255)
                    aaCGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGc                                                        <  1:460147/90‑1 (MQ=255)
                        gTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATa                                                    <  1:439649/90‑1 (MQ=255)
                          aaaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACg                                                  >  1:356964/1‑90 (MQ=255)
                          aaaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACg                                                  >  2:173260/1‑90 (MQ=255)
                           aaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                 <  1:122516/90‑1 (MQ=255)
                                gATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTc                                                           <  1:23457/75‑1 (MQ=255)
                                gATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTc                                                           >  2:23457/1‑75 (MQ=255)
                                 aTCAGCATGCTAAGCGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCa                                           <  1:716276/90‑1 (MQ=255)
                                       aTGCTAAGTGACTGGAAGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTa                                     >  1:153835/1‑90 (MQ=255)
                                          cTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCt                                  >  2:23914/1‑90 (MQ=255)
                                            aaGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGt                                >  2:697968/1‑90 (MQ=255)
                                                        gAAGCCAACTGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                    <  1:569627/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                    <  1:111166/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                    <  2:461391/90‑1 (MQ=255)
                                                              aaaTGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGt                                >  2:178847/1‑72 (MQ=255)
                                                              aaaTGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGt                                <  1:178847/72‑1 (MQ=255)
                                                                       aCTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGt     <  2:631499/90‑1 (MQ=255)
                                                                         tGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGa   <  2:356964/90‑1 (MQ=255)
                                                                          gAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGaa  >  1:460210/1‑90 (MQ=255)
                                                                                       |                                                                            
CCCGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAA  >  NZ_CP009273/1573070‑1573233

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