Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 329,082 C→T E85K (GAG→AAG)  betI ← DNA‑binding transcriptional repressor BetI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913329,0820CT97.7% 145.9 / ‑4.3 44E85K (GAG→AAG) betIDNA‑binding transcriptional repressor BetI
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/20);  total (24/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 = 9.91e-01

TAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTCTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGAAATAGTGGCTGATGATCCCCGTAGAAACGCC  >  NC_000913/328960‑329217
                                                                                                                          |                                                                                                                                       
tAAACGATAGAGCATCGGCTGAGGCGTACGGCTGGCCTAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGaa                                                                                                                        <  1:414980/140‑1 (MQ=255)
 aaaCGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAg                                                                                                                       <  2:75582/140‑1 (MQ=255)
  aaCGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTCCCGCCAACCATCACCTATAATCGCTACTTTACACTGCCCTGCGGAAgg                                                                                                                      >  2:414980/1‑139 (MQ=255)
       tAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGTCTGGAATAGCTGCTTTTCACTCCCCTGCGGAAGTGCATg                                                                                                                 >  2:273229/1‑140 (MQ=255)
       tAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATg                                                                                                                 >  1:772544/1‑140 (MQ=255)
       tAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCAAACAATCGACTGTAATAGCTGCTTTACACTGCCCTGCGGAAGTGCATg                                                                                                                 >  2:294136/1‑140 (MQ=255)
          agCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTaa                                                                                                              >  1:644700/1‑140 (MQ=255)
           gCATCGGCTGATGCATACTGCTCGCTCAGTACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAAt                                                                                                             <  1:276627/140‑1 (MQ=255)
            cATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATc                                                                                                            >  1:73158/1‑140 (MQ=255)
               cGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGAt                                                                                                         >  2:78120/1‑140 (MQ=255)
                           aCTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGt                                                                                             >  1:440860/1‑140 (MQ=255)
                               ctgGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACg                                                                                         >  2:115716/1‑140 (MQ=255)
                                 ggCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTCTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCa                                                                                       >  1:85097/1‑140 (MQ=255)
                                  gCCCAGAACGCCAGCGAGGCTTTCACCGCGGCACTGCACACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAg                                                                                      <  1:420983/140‑1 (MQ=255)
                                  gCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGGACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACcg                                                                                                >  2:320480/1‑130 (MQ=255)
                                  gCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACcg                                                                                                <  1:320480/130‑1 (MQ=255)
                                  gCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAg                                                                                      >  1:483527/1‑140 (MQ=255)
                                  gCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAg                                                                                      <  2:73158/140‑1 (MQ=255)
                                   cccAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGc                                                                                     >  1:84067/1‑140 (MQ=255)
                                   cccAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGc                                                                                     <  2:84067/140‑1 (MQ=255)
                                       gAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGAc                                                                                 >  1:142069/1‑140 (MQ=255)
                                       gAACGCCAGCCAGGCTTTCATCACCCCACCGCTCACCCGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGAc                                                                                 <  1:273229/140‑1 (MQ=255)
                                          cGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTgg                                                                              <  2:771265/140‑1 (MQ=255)
                                                 cagGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATc                                                                       <  1:851664/140‑1 (MQ=255)
                                                    gCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATcgcg                                                                    <  2:168391/140‑1 (MQ=255)
                                                      tttCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCa                                                                  <  2:142069/140‑1 (MQ=255)
                                                      tttCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCa                                                                  <  2:777575/140‑1 (MQ=255)
                                                          aTCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGt                                                              >  1:93096/1‑140 (MQ=255)
                                                                   cTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCcag                                                     <  2:644700/140‑1 (MQ=255)
                                                                              cGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCAtttt                                          >  2:846168/1‑140 (MQ=255)
                                                                                tttCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTg                                        >  2:217879/1‑140 (MQ=255)
                                                                                tttCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTg                                        >  2:669009/1‑140 (MQ=255)
                                                                                   cATCGAAGTTTCCGCCAACAATCGCCTGTAACCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTAAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTcc                                     >  1:767637/1‑140 (MQ=255)
                                                                                     tCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCt                                   >  2:631293/1‑140 (MQ=255)
                                                                                        aaGTTTCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGaa                                >  2:161595/1‑140 (MQ=255)
                                                                                           tttCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGAAATa                             <  1:115716/140‑1 (MQ=255)
                                                                                             tCCGCCAACAATCGGCTGTAATCGCTGATTTGCACTGCCCTGCGGCAGCGCATGTAATCGATTAAAAACCGCGTCACGCAGCTGACCGCTGATATCGCGCATGGTTGCTTCCAGCAGACCATTATTGTCCCTGCAATAgg                           >  2:113758/1‑139 (MQ=255)
                                                                                             tCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGAAATAGt                           <  2:93096/140‑1 (MQ=255)
                                                                                             tCCGCCAACAATCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGAAATAGt                           <  1:846168/140‑1 (MQ=255)
                                                                                               cGCCAACAATCGCCTGTAATCGATGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGAAATAGTgg                         <  1:294136/140‑1 (MQ=255)
                                                                                                       aTCGCCTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGAAATAGTGGCTgatgat                 <  2:483527/140‑1 (MQ=255)
                                                                                                           ccTGTAATCGCTGCTTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGAAATAGTGGCTGATGATcccc             <  1:217879/140‑1 (MQ=255)
                                                                                                                    gctgctTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGAAATAGTGGCTGATGATCCCCGTAGAAACg    >  1:620953/1‑140 (MQ=255)
                                                                                                                      tgctTTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCACACCATTTTGGTCACTGGAATAGTGGGTGATAATCCCCGTAGAAACGcc  >  2:236185/1‑140 (MQ=255)
                                                                                                                          |                                                                                                                                       
TAAACGATAGAGCATCGGCTGATGCATACTGCTGGCCCAGAACGCCAGCCAGGCTTTCATCGCCGCACTGCTCACCTGCGTTTCATCGAAGTTTCCGCCAACAATCGCCTGTAATCGCTGCTCTGCACTGCCCTGCGGAAGTGCATGTAATCGATTCAAAACCGCGTCACGCAGCTGACTGGTGATATCGCGCATGGTTGCTTCCAGCAGACCATTTTTGTCCCTGAAATAGTGGCTGATGATCCCCGTAGAAACGCC  >  NC_000913/328960‑329217

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