Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 854,225 C→A A312S (GCC→TCC)  ybiW ← predicted pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606854,2250CA97.5% 122.9 / ‑2.2 40A312S (GCC→TCC) ybiWpredicted pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (22/17);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

TTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATAC  >  REL606/854153‑854296
                                                                        |                                                                       
ttGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAgtg                                                                      <  2:596079/76‑1 (MQ=255)
  gAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGtt                                                                    >  1:956328/1‑76 (MQ=255)
     tgtgAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTctc                                                                 <  2:844403/76‑3 (MQ=255)
         aGCCGGAGGAGATCTAGTACACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGCTGGAGTGTTCGCGATc                                                             <  1:17731/76‑1 (MQ=255)
         aGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATc                                                             >  1:1351998/1‑76 (MQ=255)
          gCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATcc                                                            >  1:1120763/1‑76 (MQ=255)
            cGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAg                                                          >  1:912029/1‑76 (MQ=255)
             ggAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGc                                                         >  1:1004962/1‑76 (MQ=255)
                 cGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCt                                                     >  2:567855/1‑76 (MQ=255)
                  ggATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTg                                                    <  2:546146/76‑1 (MQ=255)
                  ggATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTg                                                    <  1:1225556/76‑1 (MQ=255)
                   gATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTgg                                                   >  2:1481695/1‑76 (MQ=255)
                     tCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGtt                                                 <  1:1159024/76‑1 (MQ=255)
                       ttGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGa                                               >  2:378732/1‑76 (MQ=255)
                       ttGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGa                                               >  2:210078/1‑76 (MQ=255)
                        tGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAg                                              >  1:1215736/1‑76 (MQ=255)
                        tGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAg                                              >  2:377755/1‑76 (MQ=255)
                        tGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGACGGAGTGTTAGCGATCCAGCGTCTGGTCGAg                                              >  2:105023/1‑76 (MQ=255)
                          ttCACTTCCAGCAGNTTCAGCCAGCAGCTATGCAGCATCTCGGTGGAGTGTTCGCGATCCTGCCTCTGGTTGAGtt                                            <  1:1477488/76‑1 (MQ=255)
                                 ccAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTc                                     <  1:678133/76‑1 (MQ=255)
                                          tCAGCCAGCAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATa                            >  1:670903/1‑76 (MQ=255)
                                           cagccagcAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAg                           >  2:1223501/1‑76 (MQ=255)
                                             gccagcAGCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTa                         <  1:1091509/76‑1 (MQ=255)
                                               cagcagCTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACg                       <  2:1303800/76‑1 (MQ=255)
                                                  cagcTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGt                    >  1:361897/1‑76 (MQ=255)
                                                    gcTATGCAGCATCTCGATGNANTGTTCGCNATCCAGCGTCTGGGTGAGTTCNACNTCACGGCGATAGTACGGGTag                  <  2:248447/76‑1 (MQ=255)
                                                     cTATGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTaga                 >  1:1025747/1‑76 (MQ=255)
                                                       aTGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGgggaggga               >  2:479775/1‑76 (MQ=255)
                                                       aTGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGgagaga               >  2:479754/1‑76 (MQ=255)
                                                       aTGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAgaga               >  2:515232/1‑76 (MQ=255)
                                                       aTGCAGCATCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAgaga               >  2:961931/1‑76 (MQ=255)
                                                          cagcaTCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGCAGAGATAc            >  2:79089/1‑76 (MQ=255)
                                                           agcaTCTCGATGGAGTGTTCGCGATTCACCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACt           <  1:260449/76‑1 (MQ=255)
                                                            gcaTCTCTGTGGTGCGTTCTCGATCTAGCTTCTGGTTGAATTCAACGTCGCGGCGATAGTACGGGTATAGATATTg          <  2:34703/76‑1 (MQ=255)
                                                             caTCTCGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTgg         >  2:338375/1‑76 (MQ=255)
                                                               tctcGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTc       <  1:1159133/76‑1 (MQ=255)
                                                                 tcGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCa     <  1:222565/76‑1 (MQ=255)
                                                                  cGATGGAGTGTTCGCTATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCAt    <  1:973130/76‑1 (MQ=255)
                                                                  cGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCAt    >  2:1350299/1‑76 (MQ=255)
                                                                  cGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCAt    <  1:953783/76‑1 (MQ=255)
                                                                  cGATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCAt    <  2:769406/76‑1 (MQ=255)
                                                                   gATGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATa   <  2:1423040/76‑1 (MQ=255)
                                                                    aTGGAGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATAc  <  1:1623194/76‑1 (MQ=255)
                                                                        |                                                                       
TTGAGTGTGAGCCGGAGCGGATCTTGTTCACTTCCAGCAGTTTCAGCCAGCAGCTATGCAGCATCTCGATGGCGTGTTCGCGATCCAGCGTCTGGTTGAGTTCAACGTCGCGGCGATAGTACGGGTAGAGATACTGGTCCATAC  >  REL606/854153‑854296

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: