Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 676,403 C→T A49T (GCC→ACC)  ybeQ ← Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913676,4030CT98.0% 163.1 / ‑4.1 50A49T (GCC→ACC) ybeQSel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/22);  total (28/22)
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.90e-01

TTATGACGGCTTTTATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGCTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAGCACGTTCGA  >  NC_000913/676309‑676498
                                                                                              |                                                                                               
ttATGACGGCTTTTATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGtttttt                                                                                            <  1:412833/100‑1 (MQ=255)
 tATGACGGCTTTTATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTtc                                                                                           <  2:1281341/100‑1 (MQ=255)
  aTGACGGCTTTTATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTATAAGGTTTTTtcg                                                                                          <  1:858735/100‑1 (MQ=255)
    gACGGCTTTTATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTtcgtc                                                                                        >  2:1228258/1‑100 (MQ=255)
    gACGGCTTTTATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTtcgtc                                                                                        >  1:285133/1‑100 (MQ=255)
        gCTTTTATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTcc                                                                                    >  1:194926/1‑100 (MQ=255)
         cTTTTATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCg                                                                                   <  1:378841/100‑1 (MQ=255)
            ttATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGgag                                                                                >  1:786152/1‑100 (MQ=255)
                ccTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATc                                                                            >  1:852489/1‑100 (MQ=255)
                 cTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCa                                                                           >  2:255833/1‑100 (MQ=255)
                   cGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAAt                                                                         <  2:260667/100‑1 (MQ=255)
                     gTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTg                                                                       >  1:51835/1‑100 (MQ=255)
                       ataCTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTgcg                                                                     <  1:13397/100‑1 (MQ=255)
                        taCTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTgcgc                                                                    <  1:1268649/100‑1 (MQ=255)
                          cTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAAACAATTGCGCTa                                                                  >  1:464970/1‑100 (MQ=255)
                             ttGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATcg                                                               <  2:568561/100‑1 (MQ=255)
                              tGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATcgc                                                              <  2:115812/100‑1 (MQ=255)
                               gCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATcgcg                                                             >  1:1148389/1‑100 (MQ=255)
                                cccTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATcgc                                                              >  2:469480/1‑98 (MQ=255)
                                cccTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGa                                                            >  2:539506/1‑100 (MQ=255)
                                 ccTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGAt                                                           <  1:137709/100‑1 (MQ=255)
                                 ccTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGAt                                                           >  2:44297/1‑100 (MQ=255)
                                 ccTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGAt                                                           <  2:595199/100‑1 (MQ=255)
                                 ccTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGAt                                                           <  2:595231/100‑1 (MQ=255)
                                  cTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGAt                                                           <  2:175316/99‑1 (MQ=255)
                                      ggAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTataa                                                      >  2:991805/1‑100 (MQ=255)
                                          tACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGCTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATaaa                                                  >  1:419307/1‑100 (MQ=255)
                                            cTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAAc                                                <  1:105494/100‑1 (MQ=255)
                                                   tCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCcaacaa                                         <  1:1228269/100‑1 (MQ=255)
                                                    cacaATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAat                                        >  1:377498/1‑100 (MQ=255)
                                                     acaATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAACACCCAACAata                                       <  2:1220989/100‑1 (MQ=255)
                                                     acaATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAata                                       <  1:1024186/100‑1 (MQ=255)
                                                     acaATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAata                                       <  2:117048/100‑1 (MQ=255)
                                                      caATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAatat                                      >  1:1212373/1‑100 (MQ=255)
                                                          gACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTg                                  >  2:216176/1‑100 (MQ=255)
                                                             cTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGACTCAATTGCGATATCGCGATTATAATAAAACCCAAAAATATATTgcgc                               >  2:429203/1‑100 (MQ=255)
                                                                    aGCGGCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCaa                        >  1:1262241/1‑100 (MQ=255)
                                                                       ggCCAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCg                     >  1:1034034/1‑100 (MQ=255)
                                                                         ccAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGcc                   >  2:197821/1‑100 (MQ=255)
                                                                          cAGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCt                  >  2:891463/1‑100 (MQ=255)
                                                                           aGCTTCAGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCtt                 >  1:390991/1‑100 (MQ=255)
                                                                                 aGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAg           <  1:1063986/100‑1 (MQ=255)
                                                                                 aGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAg           >  1:27517/1‑100 (MQ=255)
                                                                                 aGCCAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAg           <  2:1004548/100‑1 (MQ=255)
                                                                                  gCCAGTAAAAGGTTTTTTGGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAGc          >  1:736575/1‑100 (MQ=255)
                                                                                  gCCAGTAAAAGGTTTTTTCGTCGTCCGTAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAGc          >  2:526312/1‑100 (MQ=255)
                                                                                   ccAGTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAGCa         <  2:1212362/100‑1 (MQ=255)
                                                                                      gTAAAAGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAGCACGt      >  1:745968/1‑100 (MQ=255)
                                                                                        aaaaGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAGCACGTTc    <  1:854355/100‑1 (MQ=255)
                                                                                         aaaGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAGCACGTTCg   <  1:1029577/100‑1 (MQ=255)
                                                                                          aaGGTTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAGCACGTTCGa  >  2:684306/1‑100 (MQ=255)
                                                                                              |                                                                                               
TTATGACGGCTTTTATCCTCGGTATACTTTTGCCCTAAGGAATACTGTGCTTCACAATGACCTTGCTCAGCGGCCAGCTTCAGCCAGTAAAAGGCTTTTTCGTCGTCCGGAGAATCAATTGCGCTATCGCGATTATAATAAAACCCAACAATATATTGAGCCTCGCAATCGCCTTTTTCAGCACGTTCGA  >  NC_000913/676309‑676498

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