Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,505,986 T→C intergenic (‑283/+7) BW25113_RS17530 ← / ← trpS GntR family transcriptional regulator YhfZ/tryptophan‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,505,9860TC93.6% 132.0 / 2.2 47intergenic (‑283/+7)BW25113_RS17530/trpSGntR family transcriptional regulator YhfZ/tryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (32/12);  total (34/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

TTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTCC  >  NZ_CP009273/3505900‑3506068
                                                                                      |                                                                                  
tttCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGt                                                                                 <  2:1434844/90‑1 (MQ=255)
  tCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGa                                                                               >  1:176891/1‑90 (MQ=255)
  tCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGa                                                                               >  1:671289/1‑90 (MQ=255)
     ttGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAAtt                                                                            >  2:1232325/1‑90 (MQ=255)
     ttGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAAtt                                                                            >  2:424874/1‑90 (MQ=255)
     ttGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAAtt                                                                            >  1:765741/1‑90 (MQ=255)
       gCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTAc                                                                          >  2:260893/1‑90 (MQ=255)
          atatGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGc                                                                       >  1:255854/1‑90 (MQ=255)
            atGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCtt                                                                     >  2:399930/1‑90 (MQ=255)
            atGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCtt                                                                     >  2:1580676/1‑90 (MQ=255)
                 cATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCca                                                                >  2:713656/1‑90 (MQ=255)
                    tACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACaa                                                             <  1:164257/90‑1 (MQ=255)
                      cTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGt                                                                                 <  1:1054788/68‑1 (MQ=255)
                      cTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGt                                                                                 >  2:1054788/1‑68 (MQ=255)
                      cTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACaaaa                                                           >  2:1389478/1‑90 (MQ=255)
                        ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAAcc                                                         >  2:905844/1‑90 (MQ=255)
                        ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAAcc                                                         >  2:598861/1‑90 (MQ=255)
                          gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCaa                                                       >  1:1357771/1‑90 (MQ=255)
                           tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAAt                                                      <  1:563888/90‑1 (MQ=255)
                           tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAAt                                                      <  2:362285/90‑1 (MQ=255)
                           tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAAt                                                      <  1:239555/90‑1 (MQ=255)
                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCg                                                    >  2:986009/1‑90 (MQ=255)
                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCg                                                    <  2:840929/90‑1 (MQ=255)
                               ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCt                                                  >  2:778111/1‑90 (MQ=255)
                               ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCt                                                  >  2:116402/1‑90 (MQ=255)
                                      aGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTaca                                           >  2:1514158/1‑90 (MQ=255)
                                       ggCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTacac                                          >  1:1077450/1‑90 (MQ=255)
                                        gCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACAcc                                         <  2:1007158/90‑1 (MQ=255)
                                         cGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCg                                        >  2:1420876/1‑90 (MQ=255)
                                           ttCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCt                                      >  2:1385454/1‑90 (MQ=255)
                                           ttCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCt                                      >  1:779560/1‑90 (MQ=255)
                                            tCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCtt                                     >  1:1187208/1‑90 (MQ=255)
                                                   gcATCCGGCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGt                              <  1:230699/90‑1 (MQ=255)
                                                         ggCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGa                        >  1:513729/1‑90 (MQ=255)
                                                          gCATGAACAAAGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGaa                       <  1:778111/90‑1 (MQ=255)
                                                             tGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCg                    >  2:810329/1‑90 (MQ=255)
                                                                   aaGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcg              >  2:1307669/1‑90 (MQ=255)
                                                                   aaGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcg              >  1:682256/1‑90 (MQ=255)
                                                                   aaGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcg              >  1:808590/1‑90 (MQ=255)
                                                                    aGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcgc             <  1:1499360/90‑1 (MQ=255)
                                                                    aGCGCAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcgc             <  2:699210/90‑1 (MQ=255)
                                                                        cAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACAGCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGc         >  2:1195462/1‑90 (MQ=255)
                                                                        cAATTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGc         >  1:318379/1‑90 (MQ=255)
                                                                          aTTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcgc             >  2:1421029/1‑84 (MQ=255)
                                                                          aTTTGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTgcgcgc             <  1:1421029/84‑1 (MQ=255)
                                                                             tGCCAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCttttt    >  2:730502/1‑90 (MQ=255)
                                                                               ccAGCAACAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTcc  >  2:992969/1‑90 (MQ=255)
                                                                                      |                                                                                  
TTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTCC  >  NZ_CP009273/3505900‑3506068

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