Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 502,414 T→C D42G (GAT→GGT)  ybaK ← Cys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273502,4140TC98.2% 193.7 / ‑4.2 54D42G (GAT→GGT) ybaKCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (33/20);  total (34/20)
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.13e-01

TCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACG  >  NZ_CP009273/502327‑502498
                                                                                       |                                                                                    
tCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGAccc                                                                                    <  1:440058/90‑1 (MQ=255)
tCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGAccc                                                                                    <  2:1248817/90‑1 (MQ=255)
       gACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCa                                                                             >  1:490125/1‑90 (MQ=255)
       gACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCa                                                                             <  2:1744159/90‑1 (MQ=255)
           tGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAAcc                                                                         <  2:2762667/90‑1 (MQ=255)
            gCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCt                                                                        >  1:2853448/1‑90 (MQ=255)
            gCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCt                                                                        >  2:122695/1‑90 (MQ=255)
                            ccACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACtt                                                        >  1:2118129/1‑90 (MQ=255)
                            ccACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACtt                                                        >  2:1136634/1‑90 (MQ=255)
                            ccACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACtt                                                        >  1:2700805/1‑90 (MQ=255)
                             cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                                       >  1:2677130/1‑90 (MQ=255)
                             cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                                       >  1:2444281/1‑90 (MQ=255)
                             cACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTc                                                       >  1:242398/1‑90 (MQ=255)
                                ggCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATc                                                    <  2:1678458/90‑1 (MQ=255)
                                    aGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCa                                                >  2:1874627/1‑90 (MQ=255)
                                     ggTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCaa                                               >  2:2424756/1‑90 (MQ=255)
                                        gtTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAt                                            <  1:372019/90‑1 (MQ=255)
                                         tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                                           >  2:1671318/1‑90 (MQ=255)
                                         tttCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAAtt                                           >  2:1012265/1‑90 (MQ=255)
                                             atatCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                                       <  2:2060741/90‑1 (MQ=255)
                                             atatCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                                       <  2:1351015/90‑1 (MQ=255)
                                             atatCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                                       <  1:1373628/90‑1 (MQ=255)
                                                tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATTGCCAAAATTGGTTTCa                                    >  2:1330912/1‑90 (MQ=255)
                                                tCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCa                                    >  1:2121677/1‑90 (MQ=255)
                                                  aCCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGc                                  >  2:2703111/1‑90 (MQ=255)
                                                  aCCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGc                                  >  1:1598250/1‑90 (MQ=255)
                                                  aCCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGc                                  >  2:1364135/1‑90 (MQ=255)
                                                  aCCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGc                                  <  2:1695024/90‑1 (MQ=255)
                                                   ccGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGcc                                 <  2:737014/90‑1 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                                       >  1:678535/1‑82 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                                       >  1:3025432/1‑82 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                                       <  2:678535/82‑1 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGtt                                       <  2:3025432/82‑1 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCgg                               <  2:2791699/90‑1 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCgg                               >  1:2376545/1‑90 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCgg                               >  2:1660695/1‑90 (MQ=255)
                                                     gTTCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCgg                               >  1:2284124/1‑90 (MQ=255)
                                                       tCACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGAt                             >  1:452349/1‑90 (MQ=255)
                                                        cACTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATc                            <  2:1003999/90‑1 (MQ=255)
                                                         aCTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCg                           <  2:490125/90‑1 (MQ=255)
                                                         aCTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCg                           <  1:2768639/90‑1 (MQ=255)
                                                          cTGCCACCAGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGt                          >  2:200293/1‑90 (MQ=255)
                                                             ccaccaGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCt                       >  2:2350040/1‑90 (MQ=255)
                                                               accaGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCg                     >  1:2583096/1‑90 (MQ=255)
                                                                ccaGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGt                    <  2:452349/90‑1 (MQ=255)
                                                                ccaGCAGCGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGt                    <  2:205571/90‑1 (MQ=255)
                                                                      gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt              <  2:563059/90‑1 (MQ=255)
                                                                      gcGTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTAt              >  1:38196/1‑90 (MQ=255)
                                                                        gTTTTGTAGACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATgg            >  1:2012566/1‑90 (MQ=255)
                                                                                gACCTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGaaa    <  1:122695/90‑1 (MQ=255)
                                                                                 aCCTGACCCGGATTCAACCCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAAc   <  2:1237380/90‑1 (MQ=255)
                                                                                  ccTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACg  >  2:256649/1‑90 (MQ=255)
                                                                                  ccTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACg  >  2:1773781/1‑90 (MQ=255)
                                                                                  ccTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACg  >  1:2184264/1‑90 (MQ=255)
                                                                                  ccTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACg  >  1:2373634/1‑90 (MQ=255)
                                                                                  ccTGACCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACg  >  1:374332/1‑90 (MQ=255)
                                                                                       |                                                                                    
TCCAGTTGACCTGCGACCGGCGTAACGGCCACGGCAAGGTGTTTCATATCACCGTTCACTGCCACCAGCAGCGTTTTGTAGACCTGATCCGGATTCAAACCTAATTTTTTGACGACTTCATCGCCAAAATTGGTTTCAGCCGGATCGTGCTCGTAGGTATGGATTTGAAACG  >  NZ_CP009273/502327‑502498

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