Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,072,858 G→A 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,8580GA97.4% 107.5 / NA 39A238A (GCC→GCT) yihU3‑sulfolactaldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/21);  total (17/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 = 1.00e+00

CAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGC  >  NC_000913/4072723‑4072997
                                                                                                                                       |                                                                                                                                           
cAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAgat                                                                                                                                      >  2:420020/1‑143 (MQ=255)
         ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGa                                                                                                                             >  1:648137/1‑143 (MQ=255)
          gACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGCTCGATCACGaa                                                                                                                            >  2:551739/1‑143 (MQ=255)
                    ggCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGaa                                                                                                                  <  1:396119/143‑1 (MQ=255)
                       gaccgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGa                                                                                                               >  1:513537/1‑143 (MQ=255)
                           gcgaccCGCTGCCCGCGCGTGGCTATAACCCTCCCGGTAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTTGCGACATCAAGGGCGATGCCAAGATCCTTCTGAGCAAGATCGATCATGAAGGCGGGAGAAAGAtctc                                                                                                           <  2:331955/143‑1 (MQ=255)
                            cgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAGGATCGATCATGAAGGCGGGAGGAAGATCTcc                                                                                                          >  2:377050/1‑143 (MQ=255)
                             gccccGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGGTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCg                                                                                                         <  1:387793/141‑1 (MQ=255)
                             gaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCg                                                                                                         <  1:75035/143‑1 (MQ=255)
                                ccGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTg                                                                                                      <  1:410368/143‑1 (MQ=255)
                                 cGCTGCGCGCGCGTCGCTATAAACCTCCCGTGAGGCGGCCGCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGa                                                                                                     <  2:88509/143‑1 (MQ=255)
                                 cGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGAGGGAGAAAGGTCTCCGCTGa                                                                                                     >  1:483687/1‑143 (MQ=255)
                                                  tataAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAGAGATCTCCGCTGAGGACTTTGTTTGCCCAg                                                                                    >  1:597471/1‑143 (MQ=255)
                                                        ccTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACCTCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgtg                                                                              <  2:340011/143‑1 (MQ=255)
                                                          tCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGGTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgtggt                                                                            <  1:143171/143‑1 (MQ=255)
                                                              gccaGGCGGCCCCCAGCGGCATTGGGACATGCAGCGGGTGGGCGACATCAAGGGCGATGCCACGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAg                                                                        <  1:199099/140‑1 (MQ=255)
                                                                   gcggcCCACAGCTGTGTTGGCACATGCAGCGGGGTGGCGACATCAAGGGCGATGCCAAGATGCTTATGAGGAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGAcc                                                                   <  2:424795/143‑1 (MQ=255)
                                                                        ccccAGCGGCATTGGCAGCTTCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTAc                                                              <  2:603894/143‑1 (MQ=255)
                                                                          ccAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACCTTGTTTGGCCAGGAAGAGGAGAAGGGACCATTACcg                                                            >  2:407734/1‑143 (MQ=255)
                                                                              gggCATTGGCACATGCCGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcgg                                                        <  2:617530/142‑1 (MQ=255)
                                                                                gCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcggcg                                                      >  1:498961/1‑143 (MQ=255)
                                                                                      gCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGCCGCCGGTaaag                                                <  1:535610/143‑4 (MQ=255)
                                                                                           tGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtcat                                           <  1:420020/143‑1 (MQ=255)
                                                                                                cTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACtt                                      >  1:354745/1‑143 (MQ=255)
                                                                                                 tGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACttt                                     <  2:513537/143‑1 (MQ=255)
                                                                                                        cGACATCAAGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCAGGAGAAAGATCTCCGATGAGGACTTTGTCTGGCCAGGAAGCGGTGAAACCAGCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCa                              >  2:136721/1‑143 (MQ=255)
                                                                                                               aaGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGa                       >  1:109900/1‑143 (MQ=255)
                                                                                                                aGGGCGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGaa                      >  1:197508/1‑143 (MQ=255)
                                                                                                                    cGATGCCAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGCGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTCCCGGCGGCGGTACCGCTCATCGCATTGGCGGCAACATCGAAGGGa                  >  2:252196/1‑143 (MQ=255)
                                                                                                                         ccAAGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCTCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCCGCGGTACCGCTCATACTTTTGCCGGCGGAATCGCATGGAAGAtt             >  1:148411/1‑143 (MQ=255)
                                                                                                                           gtgtgccTGATGGGGACGATCGATACTGAAGGCGGGAGAGAGATCTCCGCTGAGGACTTTGTTTGGCCAGGCATTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCa           <  2:48196/138‑1 (MQ=255)
                                                                                                                           aaGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGACGGGAGGATGCa           >  2:638551/1‑143 (MQ=255)
                                                                                                                           aaGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCa           <  1:282504/143‑1 (MQ=255)
                                                                                                                           aaGATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCa           <  1:580222/143‑1 (MQ=255)
                                                                                                                           aaGATCCCTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCa           <  1:283541/143‑1 (MQ=255)
                                                                                                                             gATCCTTATGAGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAgg         <  1:125749/143‑1 (MQ=255)
                                                                                                                             gATCCTTATGAGCAAGAGCGATCATGAAGGCGAGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGAAAGTGGTGAAGTGACCTTTACCCGCAGCGGTACCGCGCATCACCTGGACGGCTACATCGAAGAGACGAttgaag         >  2:267565/1‑139 (MQ=255)
                                                                                                                                 cctATGAGCAAGGTCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCtt     <  2:253023/141‑1 (MQ=255)
                                                                                                                                    aTGAGCCAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGc  <  1:224561/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACGGCAACATCGAAGGGAAGATTCAGGGCTTCGC  >  NC_000913/4072723‑4072997

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