Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,639,576 (A)6→5 coding (145/336 nt) uspB ← universal stress (ethanol tolerance) protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,639,5710A.96.1% 216.1 / 2.2 51coding (150/336 nt)uspBuniversal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (22/27);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

ACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAGCGCGCCA  >  NC_000913/3639472‑3639665
                                                                                                   |                                                                                              
acaGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTa                                                                                                >  1:2583/1‑100 (MQ=255)
 caGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTaa                                                                                               >  2:1126525/1‑100 (MQ=255)
   gcggcgAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTTGCCATGTGAGGTaaaa                                                                                             <  1:998877/100‑1 (MQ=255)
      gcgAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTaaaaagn                                                                                          <  1:537136/100‑3 (MQ=255)
      gcgAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTaaaaagc                                                                                          >  1:1203560/1‑98 (MQ=255)
       cgAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTaaaaagcc                                                                                         >  1:1271534/1‑97 (MQ=255)
        gAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTaaaaagccc                                                                                        <  1:144295/100‑5 (MQ=255)
             aaTTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTcc                                                                                  >  1:125761/1‑100 (MQ=255)
                    catcGTGATGATCGCGATAACGTTGGGCATTGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACa                                                                           <  2:972515/100‑1 (MQ=255)
                    catcGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACATGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACa                                                                           >  1:1249300/1‑100 (MQ=255)
                          gatgatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGa                                                                     >  2:1019718/1‑100 (MQ=255)
                           atgatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGAt                                                                    <  2:183352/100‑1 (MQ=255)
                           atgatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGAt                                                                    <  2:395062/100‑1 (MQ=255)
                            tgatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATa                                                                   >  1:1287854/1‑100 (MQ=255)
                             gatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATag                                                                  >  1:389101/1‑100 (MQ=255)
                                 gcgATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCa                                                              >  2:1195264/1‑100 (MQ=255)
                                        cGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATc                                                       <  1:68027/100‑1 (MQ=255)
                                           tGGGGATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGc                                                     <  2:1146641/99‑1 (MQ=255)
                                           tGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGc                                                     <  2:2583/100‑1 (MQ=255)
                                             ggCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAg                                                   <  2:1083274/99‑1 (MQ=255)
                                             ggCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGt                                                  >  1:177269/1‑100 (MQ=255)
                                               cATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTa                                                <  1:459236/100‑1 (MQ=255)
                                                  aGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACg                                              >  1:219678/1‑100 (MQ=255)
                                                   natatACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACg                                              >  2:671678/2‑98 (MQ=255)
                                                   natatACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCa                                            >  2:671693/2‑100 (MQ=255)
                                                       taCCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTAc                                        <  2:1003021/100‑1 (MQ=255)
                                                        aCCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTAcc                                       <  1:72758/100‑1 (MQ=255)
                                                          cAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCAc                                     <  1:623325/100‑1 (MQ=255)
                                                          cAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCAc                                     >  2:273684/1‑100 (MQ=255)
                                                               gagaCGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTaaca                                >  1:17624/1‑100 (MQ=255)
                                                                agaCGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTaacaa                               >  2:1170061/1‑100 (MQ=255)
                                                                agaCGCACCTGTTCGTTGGGTTGGCCATGTGAGGT‑AAAAAGCACCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTaacaa                               >  1:882874/1‑100 (MQ=255)
                                                                  aCGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcg                             <  1:937333/100‑1 (MQ=255)
                                                                  aCGCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcg                             <  2:1022368/100‑1 (MQ=255)
                                                                    gCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcgcg                           <  2:964685/100‑1 (MQ=255)
                                                                    gCACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcgcg                           <  2:601985/100‑1 (MQ=255)
                                                                     cACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcgcgc                          >  1:244843/1‑100 (MQ=255)
                                                                     cACCTGTTTGTTGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGCACCACTAACAAcgcgc                          >  2:760994/1‑100 (MQ=255)
                                                                           tttgttGGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtg                    <  2:882867/100‑1 (MQ=255)
                                                                           tttgttGGGTTGGCCATGTGAAGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtg                    <  1:1081420/100‑1 (MQ=255)
                                                                            ttgttgGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAGCATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtga                   >  2:714383/1‑100 (MQ=255)
                                                                            ttgttgGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAGCATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtga                   >  2:714290/1‑100 (MQ=255)
                                                                            ttgttgGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtga                   <  1:1245957/100‑1 (MQ=255)
                                                                            ttgttgGGTTGGCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtga                   >  1:192702/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 <  1:1126535/100‑1 (MQ=255)
                                                                              gttgggttggCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 >  1:833547/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 <  2:790544/100‑1 (MQ=255)
                                                                              gttgggttggCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 >  2:676176/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 >  1:428157/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 >  1:689087/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 <  2:256090/100‑1 (MQ=255)
                                                                                 gggttggCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATgaga              <  1:1126000/100‑1 (MQ=255)
                                                                                    ttggCCATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAAt           <  2:556367/100‑1 (MQ=255)
                                                                                         cATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAgcgc      <  2:1271523/100‑1 (MQ=255)
                                                                                         cATGTGAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAgcgc      <  1:747407/100‑1 (MQ=255)
                                                                                           tgtgAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAgcgcgc    <  2:1287843/100‑1 (MQ=255)
                                                                                             tgAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAGCGCGCCa  <  1:734510/100‑1 (MQ=255)
                                                                                             tgAGGT‑AAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAGCGCGCCa  <  2:173395/100‑1 (MQ=255)
                                                                                                   |                                                                                              
ACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAGCGCGCCA  >  NC_000913/3639472‑3639665

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: ‑