Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 800,807 G→T 100% A17S (GCG→TCG)  ybhH → putative isomerase YbhH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913800,8070GT80.6% 63.3 / 8.3 37A17S (GCG→TCG) ybhHputative isomerase YbhH
Reads supporting (aligned to +/- strand):  ref base G (1/6);  new base T (12/16);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTATTAATCGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCC  >  NC_000913/800667‑800941
                                                                                                                                            |                                                                                                                                      
ttattaATCGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAAGTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGctcg                                                                                                                                      <  1:668805/143‑1 (MQ=255)
 taataATTGTTGCGCTAAAGTCAGTGGAGATTCTTTTCATTAAAGATAACCAGCTCGGCATCAGCGGTGTTTTTAAATTGAGGGCATTATTATTAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGt                                                                                                                                     <  2:250542/140‑1 (MQ=255)
     aaccGTTGCGCTAAAGTAAATAGAGACGCTTTGCACTGAATGTAACCAAATACCAAAAAGATATGTATTTAAATAGAGAGCAGAATAAAAAAAAAAAGACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCt                                                                                                                                 <  2:1000604/140‑1 (MQ=255)
      aTCGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAGATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCCGGTTCCTg                                                                                                                                >  2:845353/1‑143 (MQ=255)
      aTCGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTg                                                                                                                                <  1:367062/143‑1 (MQ=255)
       tCGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCAAATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGt                                                                                                                               <  1:388400/143‑1 (MQ=255)
       gcGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTCTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGt                                                                                                                               <  1:409952/142‑1 (MQ=255)
          ttGAGCTAAAGTCATTAGAGATGCTTTGCACTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTAGAGAAGCTCGTTCATGTTaa                                                                                                                            >  2:153350/1‑142 (MQ=255)
                tAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAAc                                                                                                                      <  1:546277/143‑1 (MQ=255)
                 aaaGTCATGAGAGATGCTTTGCCCTTAATGCAACCATAGCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACa                                                                                                                     <  2:167398/143‑1 (MQ=255)
                   aGTCATTGGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACATCGAGGGGCTCGTTCCTGTTAGTGGAACAtt                                                                                                                   >  1:339529/1‑143 (MQ=255)
                      cATTAGAGATGCATGGCACTTAATGTAACCATAGCGCAATAAGTTATGTTTTTAAATCGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGCGGGGCTCGTTCCTGTTAGCGGAACATTTAc                                                                                                                <  1:770332/143‑1 (MQ=255)
                        ttAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGTTGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGCACATTTAccc                                                                                                              >  2:195629/1‑143 (MQ=255)
                                  tttGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCaaa                                                                                                    <  2:19884/143‑1 (MQ=255)
                                         ttAAGGTAACCATATCGAAATATGTTATGCTTTTAAATTGAGGGCATTATTATGAAAAAAATAGCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAgcg                                                                                             <  2:860749/143‑1 (MQ=255)
                                          ttaTGTCACCTTATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAgcgc                                                                                            <  2:632020/141‑1 (MQ=255)
                                             tGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGAt                                                                                         >  2:658572/1‑143 (MQ=255)
                                                        cGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGAt                                                                              <  1:823722/143‑1 (MQ=255)
                                                          cAATAAGCTATGCTTTTAAATTGCTGGCATTCTTATGCAAACAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATgg                                                                            <  1:296746/143‑1 (MQ=255)
                                                               aGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCCGCGCGATAAAATATTGGTGGCAAtt                                                                       >  1:875924/1‑143 (MQ=255)
                                                                 tgaTGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTAt                                                                     <  1:13566/141‑1 (MQ=255)
                                                                      tttttAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGGGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGGTCGATAAGATATTGATGGCAATTATGGGtt                                                                >  1:955419/1‑143 (MQ=255)
                                                                          tAAAGTGAGGGCATCACTACGAAACAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCgg                                                            <  2:161469/143‑1 (MQ=255)
                                                                             atggAGGGCATGATGAGGAAAAAGATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTaa                                                         <  2:636382/140‑1 (MQ=255)
                                                                                       ttattattaAAAAAGTACCCTGCGGGATGTTGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACAGTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGaa                                               <  2:11428/143‑1 (MQ=255)
                                                                                        tattatGAAACAAATACCCTGCGTGATGATGCGGGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGaaa                                              <  1:326604/143‑1 (MQ=255)
                                                                                          ttatGAAAAAAATACCCTGCGTGATGATGCGAGGTGGTAACTCGGGGAGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCACACGCAGCGCGATCAAATATTGATGGCAATACTGGGCTCCGGTAACGCTCTGGagcgg                                            >  2:828288/1‑139 (MQ=255)
                                                                                             tGAAAAATATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGGAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGAc                                         <  2:742697/143‑1 (MQ=255)
                                                                                             ggAAAATAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGAc                                         <  2:684527/142‑1 (MQ=255)
                                                                                               aaaaaaaTACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTCCCCGAAGCTCAAACGCAGCGCGATAAAATATTGATCGCAAATATCGGTTCCGGTAACGACCTGCAAATTAAccg                                       >  1:814888/1‑141 (MQ=255)
                                                                                               aaaaaaaTAACCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACAGTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGCTCTGGAAATTGCCgg                                       >  1:700808/1‑143 (MQ=255)
                                                                                                       aCCCTGCGTGATGATGCGAGGCGTAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTggcg                               <  2:892693/143‑1 (MQ=255)
                                                                                                              gTGATGATGCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTaa                        >  1:22676/1‑143 (MQ=255)
                                                                                                                  tgatgCGAGGTGGAACCTCGCGGGGCTCGTTCCTGTTAGCGGAACCTTTACCCGAAGATCAAACGAAGCGCAATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCa                    >  1:10646/1‑143 (MQ=255)
                                                                                                                     tgggAGGTGGACCCTCGGGGTGCCCGTTCCTGTTAGCGGAACATTTCCCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTg                 <  1:703761/140‑1 (MQ=255)
                                                                                                                      gCGAGGTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGa                >  2:907307/1‑143 (MQ=255)
                                                                                                                           gTGGAACCTCGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGAGCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGGAATTGACGGTATTGGCGGCGGTAATCCACTGAcgcgt           >  1:681000/1‑139 (MQ=255)
                                                                                                                                    cGAGGGGCTCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGcc  <  2:636882/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                      
TTATTAATCGTTGAGCTAAAGTCATTAGAGATGCTTTGCCCTTAATGTAACCATATCGCAATAAGTTATGTTTTTAAATTGAGGGCATTATTATGAAAAAAATACCCTGCGTGATGATGCGAGGTGGAACCTCGAGGGGCGCGTTCCTGTTAGCGGAACATTTACCCGAAGATCAAACGCAGCGCGATAAAATATTGATGGCAATTATGGGTTCCGGTAACGATCTGGAAATTGACGGTATTGGCGGCGGTAATCCACTGACCAGTAAAGTCGCC  >  NC_000913/800667‑800941

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