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,1570AT96.5% 106.1 / ‑3.4 29I148N (ATC→AAC) yddAABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (11/17);  total (11/18)
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.17e-01

CCCGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAAGTAGT  >  NZ_CP009273/1573070‑1573238
                                                                                       |                                                                                 
cccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTc                                                                                 <  2:535445/90‑1 (MQ=255)
 ccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCa                                                                                >  2:173703/1‑90 (MQ=255)
    cgcTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTa                                                                             <  2:256647/90‑1 (MQ=255)
       tttCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAaga                                                                          <  2:72135/90‑1 (MQ=255)
                 aaTAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTc                                                                <  1:64504/90‑1 (MQ=255)
                          aaaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACg                                                       >  1:435769/1‑90 (MQ=255)
                           aaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                      <  1:132973/90‑1 (MQ=255)
                           aaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                      <  2:253411/90‑1 (MQ=255)
                               tGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGt                                                  <  2:470859/90‑1 (MQ=255)
                                gATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTc                                                 >  2:541532/1‑90 (MQ=255)
                                  tCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAg                                               >  1:192452/1‑90 (MQ=255)
                                    aGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAggg                                             <  2:274737/90‑1 (MQ=255)
                                    aGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAggg                                             <  1:46810/90‑1 (MQ=255)
                                               tGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTaa                                                                            >  1:331096/1‑48 (MQ=255)
                                               tGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTaa                                                                            <  2:331096/48‑1 (MQ=255)
                                                    ggATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTAt                             >  1:413528/1‑90 (MQ=255)
                                                     gATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATg                            <  1:320650/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTTTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  2:551727/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATc                                        <  1:266226/75‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATc                                        >  2:266226/1‑75 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  1:146717/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  1:77597/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  1:23576/90‑1 (MQ=255)
                                                              aaaTGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTa                   >  2:45453/1‑90 (MQ=255)
                                                                     aGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGc            >  1:450383/1‑90 (MQ=255)
                                                                          gAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGaa       >  2:78220/1‑90 (MQ=255)
                                                                            gCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAagt     >  2:184969/1‑90 (MQ=255)
                                                                               ttttGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAagtagt  <  2:227701/90‑1 (MQ=255)
                                                                               ttttGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAagtagt  <  1:391062/90‑1 (MQ=255)
                                                                                       |                                                                                 
CCCGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAAGTAGT  >  NZ_CP009273/1573070‑1573238

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: