Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,072,858 G→A 42.5% A238A (GCC→GCT)  yihU ← 3‑sulfolactaldehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,072,8580GA42.5% 15.0 / 43.3 45A238A (GCC→GCT) yihU3‑sulfolactaldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (13/13);  new base A (8/11);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

TGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGCACA  >  NC_000913/4072718‑4073000
                                                                                                                                            |                                                                                                                                              
tGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGGTGCCAAGATCCTTATGGGc                                                                                                                                              >  2:367609/1‑143 (MQ=255)
       gtaccgcgcACCAGTCCTGGCGACCGCGACCCGCGGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACAGGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAgatcg                                                                                                                                       <  1:618261/138‑1 (MQ=255)
         acggCGGACGAGCCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGGGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGGTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAgatcgat                                                                                                                                     <  2:262144/141‑1 (MQ=255)
               gACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGaa                                                                                                                               <  1:722137/143‑1 (MQ=255)
                  cAGTCCTGGCGACCGGGACCCGCTGCGCGCGCCTGGCTATGAACCTCCCGTGGGGCGGCCCGCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGc                                                                                                                            <  2:127489/143‑1 (MQ=255)
                      ccgggCGCCCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGga                                                                                                                        <  2:636444/140‑1 (MQ=255)
                        tGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCAAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGgaga                                                                                                                      >  1:803029/1‑143 (MQ=255)
                         ggCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGaa                                                                                                                     >  1:822861/1‑143 (MQ=255)
                         ggCCACCGCGACCCGCTGCGCCCGCCTGGCTACAATCCTCCCGTGAGGCGGCCCCCAGCGGCATGGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGaa                                                                                                                     <  1:300251/143‑1 (MQ=255)
                           cgaccgcgaccCGCTGCGCGCGCCTGGCTATAAACGTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACGTGGCGCGGGTCGGCGACATCAAGGGAGAGGGCACGGTCCGGATGGGCAAGAACCACCCTGAAGGCggggacaag                                                                                                                   >  2:221680/1‑137 (MQ=255)
                           cgaccgcgaccCCCTGCGCGCGCCGGGCTACAACCCTCCAGTGCGATGTTGTCCAGCGACAATGGCAAAGGCTGCTTTTGGGCGACATCAAGGGCGAGGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAg                                                                                                                   <  2:185710/143‑1 (MQ=255)
                              ccgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGCTCCTTAGGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtc                                                                                                                >  1:683037/1‑143 (MQ=255)
                                gcggccCGCTGCGCGCGCCTGGTTATAATCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACAGGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTTTGAGCAAGATCGATCATGAAGGCGGGAGAAAGAtctc                                                                                                              <  2:14187/143‑1 (MQ=255)
                                 cgaccCGCTGCGCGCGACTGGCTATAAACCACCCGAGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTcc                                                                                                             <  1:433997/143‑1 (MQ=255)
                                             cgcgcCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGGTCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACttt                                                                                                 >  1:246259/1‑143 (MQ=255)
                                                gcCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgtt                                                                                              >  1:700266/1‑143 (MQ=255)
                                                gcCTGGCTATAAACCTCCCGTGAGGCGGCAACCAGCGGCGTTGGCACACGCAGCTAGCTGGCGGCAACACGGGCAATGCCAAGAGGCTTATGGCGACGCTCGATACGGAAGGCGGCAGAAAGGTCACGGCTGAGGACtatgtg                                                                                              >  2:694558/1‑142 (MQ=255)
                                                   tGGCTATAAACCTTCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTgg                                                                                           <  2:25401/143‑1 (MQ=255)
                                                            aCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGACCTCCGCTGAGGACTTTGTTTGGCCAGGAAgt                                                                                  <  2:290538/143‑1 (MQ=255)
                                                             ccTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGGGAAAGATCTCCGCTGAGGGCCTTGTTTGGCCAggcggtg                                                                                 >  2:336860/1‑138 (MQ=255)
                                                                 ccGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTGTGGCCAGGAAGTGGTGa                                                                             >  1:1172/1‑143 (MQ=255)
                                                                  ggTGAGGCGGCCCCCAGCGGCATTGGCGCATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGaa                                                                            <  2:405343/142‑1 (MQ=255)
                                                                  cGTGGGGCGGCCCCCGGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGaa                                                                            <  2:553000/143‑1 (MQ=255)
                                                                        gctgcCCCCAACGTCATTCGCACATTCAGCTGGTTCGCGACATCAAGGGCGGTGCCACGATCCATGTAAATGATGTCAAGCCTGACTGCCGGTGAAAGATCTTCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGAcc                                                                      <  2:311273/140‑1 (MQ=255)
                                                                             ccccAGCGGCATTGGCACAGGCAGCTGGGTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTAc                                                                 <  1:94008/143‑1 (MQ=255)
                                                                                 aGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTATGATTGTCCAGGAAGTGGTGAAGTGACCTTTACcgac                                                             >  1:345987/1‑141 (MQ=255)
                                                                                    ggCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcggc                                                          >  1:850344/1‑143 (MQ=255)
                                                                                         tGGCACAGGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTGTGGCCAGGAAGTTGTGAAGTGACCTTTACCGGCGGCGGTAc                                                     <  1:840026/143‑1 (MQ=255)
                                                                                         gcgaacaGGCTGCTCGTGGGCGACATCAAGGGCGGTGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGGGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTAc                                                     <  2:423736/139‑1 (MQ=255)
                                                                                           gCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCg                                                   <  1:821451/143‑1 (MQ=255)
                                                                                           gCACATGCAGCTGGTGGGCGCAATAAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTAACg                                                   <  1:311485/143‑1 (MQ=255)
                                                                                            cacaTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTCTACCGGCGGCGGTACCGc                                                  >  2:378491/1‑143 (MQ=255)
                                                                                             acaTGCAGCTGGTTGGCGACCTCAAGGGCGATGCCAAGGTCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCt                                                 <  1:537802/143‑1 (MQ=255)
                                                                                               aTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtca                                               >  1:35452/1‑143 (MQ=255)
                                                                                                 gCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTAACGGCGGCGGTACCGCtcatc                                             >  2:45347/1‑143 (MQ=255)
                                                                                                       ggTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTATGATTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGAACCGCTCATCACTTTg                                       >  1:3810/1‑143 (MQ=255)
                                                                                                             cGACATCAAGGGCGAGGCCAAGATCCTTACGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCGTTAGCGGCGGCGGGACCGATCATCACCTTGCCgcgc                                 >  2:539423/1‑140 (MQ=255)
                                                                                                              gACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTGTACCGGCGACGGCACCGCACATCACTTTGACGGCaa                                >  2:420433/1‑143 (MQ=255)
                                                                                                                     aGGGCGATGCCGAGATACGTATGGGCGAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGaa                         <  2:372019/143‑1 (MQ=255)
                                                                                                                          gctGCCAAGAGCCCTGTGAGCAAGGTCGATCATGGAGGCGGGAGGAAGGTCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACAGGCGGCGGTACCGCGCATCACTTTGACGGCAACATCGAAGGGaa                    <  2:329863/141‑1 (MQ=255)
                                                                                                                               cAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAACGATCTCCGCTGAGGACTTTTTCTGGCCAGGAGGTGGTGAAGTGGTCTTTACCGGCCGCCGTCCCGCTCTTCAGCTTGTCTGCAACATCGATGGCAAGATTc               >  2:574439/1‑143 (MQ=255)
                                                                                                                                  gATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTCTCCCGGCGGGGGTTCCCGGTATCACTTTGACGGCAACATCGAAGGGAAGATTCAgg            <  1:372757/143‑1 (MQ=255)
                                                                                                                                   aTCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAggg           >  1:638689/1‑143 (MQ=255)
                                                                                                                                    tCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGc          >  2:878706/1‑143 (MQ=255)
                                                                                                                                        tATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCg      <  1:304940/143‑1 (MQ=255)
                                                                                                                                            ggCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGCCTTGGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGCCAACAGCGAAGGGAAGACTCACGGCTTCGcagc  >  2:492951/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                              
TGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGCACA  >  NC_000913/4072718‑4073000

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