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)
             aaTTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTcc                                                                                  >  1:125761/1‑100 (MQ=255)
                    catcGTGATGATCGCGATAACGTTGGGCATTGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACa                                                                           <  2:972515/100‑1 (MQ=255)
                    catcGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACATGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACa                                                                           >  1:1249300/1‑100 (MQ=255)
                          gatgatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGa                                                                     >  2:1019718/1‑100 (MQ=255)
                           atgatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGAt                                                                    <  2:183352/100‑1 (MQ=255)
                           atgatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGAt                                                                    <  2:395062/100‑1 (MQ=255)
                            tgatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATa                                                                   >  1:1287854/1‑100 (MQ=255)
                             gatCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATag                                                                  >  1:389101/1‑100 (MQ=255)
                                 gcgATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCa                                                              >  2:1195264/1‑100 (MQ=255)
                                        cGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATc                                                       <  1:68027/100‑1 (MQ=255)
                                           tGGGGATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGc                                                     <  2:1146641/99‑1 (MQ=255)
                                           tGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGc                                                     <  2:2583/100‑1 (MQ=255)
                                             ggCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAg                                                   <  2:1083274/99‑1 (MQ=255)
                                             ggCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGt                                                  >  1:177269/1‑100 (MQ=255)
                                               cATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTa                                                <  1:459236/100‑1 (MQ=255)
                                                  aGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACg                                              >  1:219678/1‑100 (MQ=255)
                                                   natatACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACg                                              >  2:671678/2‑98 (MQ=255)
                                                   natatACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCa                                            >  2:671693/2‑100 (MQ=255)
                                                       taCCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTAc                                        <  2:1003021/100‑1 (MQ=255)
                                                        aCCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTAcc                                       <  1:72758/100‑1 (MQ=255)
                                                          cAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCAc                                     <  1:623325/100‑1 (MQ=255)
                                                          cAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCAc                                     >  2:273684/1‑100 (MQ=255)
                                                               gagaCGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTaaca                                >  1:17624/1‑100 (MQ=255)
                                                                agaCGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTaacaa                               >  2:1170061/1‑100 (MQ=255)
                                                                agaCGCACCTGTTCGTTGGGTTGGCCATGTGAGGTAAAAAGCACCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTaacaa                               >  1:882874/1‑100 (MQ=255)
                                                                  aCGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcg                             <  1:937333/100‑1 (MQ=255)
                                                                  aCGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcg                             <  2:1022368/100‑1 (MQ=255)
                                                                    gCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcgcg                           <  2:964685/100‑1 (MQ=255)
                                                                    gCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcgcg                           <  2:601985/100‑1 (MQ=255)
                                                                     cACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAAcgcgc                          >  1:244843/1‑100 (MQ=255)
                                                                     cACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGCACCACTAACAAcgcgc                          >  2:760994/1‑100 (MQ=255)
                                                                           tttgttGGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtg                    <  2:882867/100‑1 (MQ=255)
                                                                           tttgttGGGTTGGCCATGTGAAGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtg                    <  1:1081420/100‑1 (MQ=255)
                                                                            ttgttgGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAGCATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtga                   >  2:714383/1‑100 (MQ=255)
                                                                            ttgttgGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAGCATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtga                   >  2:714290/1‑100 (MQ=255)
                                                                            ttgttgGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtga                   <  1:1245957/100‑1 (MQ=255)
                                                                            ttgttgGGTTGGCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtga                   >  1:192702/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 <  1:1126535/100‑1 (MQ=255)
                                                                              gttgggttggCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 >  1:833547/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 <  2:790544/100‑1 (MQ=255)
                                                                              gttgggttggCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 >  2:676176/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 >  1:428157/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 >  1:689087/1‑100 (MQ=255)
                                                                              gttgggttggCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGtgatg                 <  2:256090/100‑1 (MQ=255)
                                                                                 gggttggCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATgaga              <  1:1126000/100‑1 (MQ=255)
                                                                                    ttggCCATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAAt           <  2:556367/100‑1 (MQ=255)
                                                                                         cATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAgcgc      <  2:1271523/100‑1 (MQ=255)
                                                                                         cATGTGAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAgcgc      <  1:747407/100‑1 (MQ=255)
                                                                                           tgtgAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAgcgcgc    <  2:1287843/100‑1 (MQ=255)
                                                                                             tgAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAGCGCGCCa  <  1:734510/100‑1 (MQ=255)
                                                                                             tgAGGTAAAAAGCCCCCTCCATCAACATATTGATAGAGCAATGGATCGCAGTTACGCAGTACCACTAACAACGCGCGTAGTGATGAGAAATAGCGCGCCa  <  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: