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.4% 202.0 / ‑2.5 56I148N (ATC→AAC) yddAABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base T (24/29);  minor base C (1/1);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAAGTAGT  >  NZ_CP009273/1573070‑1573238
                                                                                       |                                                                                 
cccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTc                                                                                 <  1:570534/90‑1 (MQ=255)
cccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTc                                                                                 <  1:34053/90‑1 (MQ=255)
cccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTc                                                                                 <  2:1405917/90‑1 (MQ=255)
cccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTc                                                                                 <  1:269512/90‑1 (MQ=255)
 ccGCGCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCa                                                                                <  1:462135/90‑1 (MQ=255)
  cgcgCTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAg                                                                               >  2:147831/1‑90 (MQ=255)
   gcgcTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGt                                                                              >  2:782908/1‑90 (MQ=255)
   gcgcTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGt                                                                              >  2:335814/1‑90 (MQ=255)
    cgcTTTCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTa                                                                             <  2:187828/90‑1 (MQ=255)
       tttCCCAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAaga                                                                          >  2:164511/1‑90 (MQ=255)
          cccAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGATCAGTAAGAGaa                                                                       <  2:41092/90‑1 (MQ=255)
            cAAAGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATg                                                                     <  1:164511/90‑1 (MQ=255)
               aGAATAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCt                                                                  <  2:1378905/90‑1 (MQ=255)
                 aaTAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTc                                                                <  2:1533043/90‑1 (MQ=255)
                 aaTAACGGTAAAGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTc                                                                <  2:1057988/90‑1 (MQ=255)
                          aaaGGTGATCAGCATGCTAAGTGACTGGTTGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACg                                                       >  1:399007/1‑90 (MQ=255)
                          aaaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGCTCAGTAAGAGAATGTCTTCAGCGATACg                                                       >  1:1021761/1‑90 (MQ=255)
                           aaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGt                                                      <  1:545921/90‑1 (MQ=255)
                           aaGGTGATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGCTCAGTAAGAGAATGTCTTCAGCGATACGt                                                      <  2:1021761/90‑1 (MQ=255)
                                gATCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTc                                                 >  1:1545449/1‑90 (MQ=255)
                                 aTCAGCATGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCa                                                <  2:1387497/90‑1 (MQ=255)
                                       aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAaga                                                                          <  2:967189/58‑1 (MQ=255)
                                       aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAaga                                                                          >  1:967189/1‑58 (MQ=255)
                                       aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTa                                          >  2:1527270/1‑90 (MQ=255)
                                       aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTa                                          >  1:373922/1‑90 (MQ=255)
                                       aTGCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTa                                          >  1:835601/1‑90 (MQ=255)
                                         gCTAAGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATc                                        <  2:1033521/90‑1 (MQ=255)
                                            aaGTGACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGt                                     >  2:639406/1‑90 (MQ=255)
                                                gACTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTAttt                                 >  1:934567/1‑90 (MQ=255)
                                                  cTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTAttttt                               >  1:957123/1‑90 (MQ=255)
                                                  cTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTAttttt                               >  2:635179/1‑90 (MQ=255)
                                                  cTGGATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTAttttt                               >  2:1410026/1‑90 (MQ=255)
                                                    ggATGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTAt                             >  2:1460869/1‑90 (MQ=255)
                                                      aTGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGc                           >  1:839483/1‑90 (MQ=255)
                                                      aTGAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGc                           >  1:1543379/1‑90 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATc                                        >  2:326406/1‑75 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATc                                        <  1:326406/75‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  2:95020/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  2:909233/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  2:399007/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  1:1081802/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  2:1545449/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  2:1037410/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  1:720755/90‑1 (MQ=255)
                                                        gAAGCCAAATGAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTc                         <  1:1375115/90‑1 (MQ=255)
                                                                  gAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATc               <  1:800059/90‑1 (MQ=255)
                                                                  gAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATc               <  2:835601/90‑1 (MQ=255)
                                                                  gAAAGACTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATc               <  1:1323497/90‑1 (MQ=255)
                                                                       aCTGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGt          <  1:267415/90‑1 (MQ=255)
                                                                       aCTGAGCGTCTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGt          <  1:1218245/90‑1 (MQ=255)
                                                                         tGAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGa        >  2:831097/1‑90 (MQ=255)
                                                                          gAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGaa       >  2:790165/1‑90 (MQ=255)
                                                                          gAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGaa       >  2:98690/1‑90 (MQ=255)
                                                                          gAGCGTTTTGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGAATTTTTATGCTCGCCGTAGATCTGCGTGaa       >  2:665280/1‑90 (MQ=255)
                                                                               ttttGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAagtagt  <  2:1246902/90‑1 (MQ=255)
                                                                               ttttGCTGTTCAGTAAGAGAATGTCTTCAGCGATACGTTGGTCAGGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAagtagt  >  1:895177/1‑90 (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: