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,1570AT94.6% 132.6 / ‑0.5 37I148N (ATC→AAC) yddAABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (14/21);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CCCGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAAGTAGTAG  >  NZ_CP009273/1573070‑1573240
                                                                                       |                                                                                   
cccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTGAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTc                                                                                   <  2:69094/90‑1 (MQ=255)
cccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTc                                                                                   <  2:327188/90‑1 (MQ=255)
                  aTAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAAt                                                                        >  1:46930/1‑83 (MQ=255)
                  aTAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAAt                                                                        <  2:46930/83‑1 (MQ=255)
                  aTAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCa                                                                 >  2:123121/1‑90 (MQ=255)
                      cGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGa                                                             <  2:307262/90‑1 (MQ=255)
                          aaaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACg                                                         >  2:732159/1‑90 (MQ=255)
                           aaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                        <  2:696560/90‑1 (MQ=255)
                           aaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                        <  2:372473/90‑1 (MQ=255)
                               tGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGt                                                    <  1:245563/90‑1 (MQ=255)
                               tGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGt                                                    <  1:88232/90‑1 (MQ=255)
                                gATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTc                                                   >  1:669769/1‑90 (MQ=255)
                                gATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTc                                                   >  1:729670/1‑90 (MQ=255)
                                  tCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAg                                                 >  2:305558/1‑90 (MQ=255)
                                       aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                        <  2:713285/78‑1 (MQ=255)
                                       aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                        >  1:713285/1‑78 (MQ=255)
                                         gCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATc                                          <  1:234562/90‑1 (MQ=255)
                                             aGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTa                                      >  2:567360/1‑90 (MQ=255)
                                                    ggATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTAt                               >  2:839149/1‑90 (MQ=255)
                                                    ggATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTAt                               >  2:757424/1‑90 (MQ=255)
                                                      aTGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGc                             <  2:207498/90‑1 (MQ=255)
                                                      aTGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGc                             <  2:345032/90‑1 (MQ=255)
                                                      aTGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGc                             <  2:409784/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                           <  2:487177/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                           <  1:305558/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                           >  1:328462/1‑90 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                           <  2:643444/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                           <  1:837212/90‑1 (MQ=255)
                                                              aaaTGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTa                     >  1:699570/1‑90 (MQ=255)
                                                                  gAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATc                 <  2:632636/90‑1 (MQ=255)
                                                                    aaGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTg               <  1:717288/90‑1 (MQ=255)
                                                                         tGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGa          <  2:681599/90‑1 (MQ=255)
                                                                         tGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGa          >  2:815107/1‑90 (MQ=255)
                                                                         tGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGa          >  1:835375/1‑90 (MQ=255)
                                                                          gAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGaa         >  1:625248/1‑90 (MQ=255)
                                                                               ttttGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAagtagt    <  2:758929/90‑1 (MQ=255)
                                                                                 ttGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAagtagtag  >  1:724567/1‑90 (MQ=255)
                                                                                       |                                                                                   
CCCGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAAGTAGTAG  >  NZ_CP009273/1573070‑1573240

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