Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 990,162 T→G E111A (GAA→GCA)  ycbC ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606990,1620TG97.7% 174.0 / ‑3.2 43E111A (GAA→GCA) ycbCconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (21/21);  total (22/21)
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.29e-01

TATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCAC  >  REL606/990091‑990234
                                                                       |                                                                        
tATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCtgcatt                                                                      >  1:1027766/1‑76 (MQ=255)
 aTTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCtgcattc                                                                     >  1:1725667/1‑76 (MQ=255)
 aTTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCtgcattc                                                                     <  2:1451005/76‑1 (MQ=255)
 aTTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCtgcattc                                                                     >  1:1704035/1‑76 (MQ=255)
  ttGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCtgcattca                                                                    <  1:1697147/76‑1 (MQ=255)
   tGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAg                                                                   >  2:1027070/1‑76 (MQ=255)
     gTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGc                                                                 <  2:1123982/76‑1 (MQ=255)
      tttttGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCg                                                                >  2:1711680/1‑76 (MQ=255)
        tttGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTg                                                              <  2:1478888/76‑1 (MQ=255)
        tttGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTg                                                              <  2:1310611/76‑1 (MQ=255)
         ttGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTgg                                                             >  2:907360/1‑76 (MQ=255)
          tGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGc                                                            >  1:943912/1‑76 (MQ=255)
             tACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGa                                                         >  2:1733052/1‑76 (MQ=255)
              aCGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGAc                                                        >  1:478411/1‑76 (MQ=255)
                gCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTg                                                      <  2:799355/76‑1 (MQ=255)
                      cGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTg                                                <  1:11292/76‑1 (MQ=255)
                      cGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTg                                                >  2:1554607/1‑76 (MQ=255)
                       gTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGa                                               >  1:385829/1‑76 (MQ=255)
                          aaGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTa                                            >  1:258862/1‑76 (MQ=255)
                             aTCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAAt                                         <  1:239597/76‑1 (MQ=255)
                              tCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAAtt                                        <  1:1706532/76‑1 (MQ=255)
                              tCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAAtt                                        <  1:790102/76‑1 (MQ=255)
                               cAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTa                                       <  2:1745333/76‑1 (MQ=255)
                               cAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTa                                       <  2:502925/76‑1 (MQ=255)
                                tgTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTag                                      >  2:1877813/2‑76 (MQ=255)
                                aGTTTTGATCCCGGATTTGCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTag                                      <  1:875313/76‑1 (MQ=255)
                                   tttGATCCCGGATTTTCGCGCCATAGGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGc                                   >  2:600016/1‑76 (MQ=255)
                                   tttGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGc                                   <  2:404503/76‑1 (MQ=255)
                                    ttGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCt                                  <  1:1793706/76‑1 (MQ=255)
                                     tGATCCCGGATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTc                                 <  1:429416/76‑1 (MQ=255)
                                            ggATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGccc                          >  1:648975/1‑76 (MQ=255)
                                            ggATTTTCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGccc                          >  1:1779376/1‑76 (MQ=255)
                                               ttttCGCGCCATAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACt                       >  2:650476/1‑76 (MQ=255)
                                                         aTAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCa             <  1:1343528/76‑1 (MQ=255)
                                                         aTAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCa             <  1:277413/76‑1 (MQ=255)
                                                         aTAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCa             <  1:408073/76‑1 (MQ=255)
                                                          tAAGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAg            <  1:1222769/76‑1 (MQ=255)
                                                           aaGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAgg           >  1:1671790/1‑76 (MQ=255)
                                                            aGCGAATACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGt          <  2:51872/76‑1 (MQ=255)
                                                                aaTACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGtatct      >  2:182605/1‑74 (MQ=255)
                                                                aaTACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATaa      <  2:1737265/76‑1 (MQ=255)
                                                                aaTACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATaa      >  2:1643313/1‑76 (MQ=255)
                                                                aaTACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCAAGGCATaa      >  1:216029/1‑76 (MQ=255)
                                                                 aTACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAac     >  2:1558547/1‑76 (MQ=255)
                                                                 aTACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAac     >  1:1793746/1‑76 (MQ=255)
                                                                  tACCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAacc    >  1:91912/1‑76 (MQ=255)
                                                                   aCCTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAacca   >  1:1746777/1‑76 (MQ=255)
                                                                    ccTGCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAaccac  <  2:1441213/76‑1 (MQ=255)
                                                                       |                                                                        
TATTGGTTTTTGCTACGCCTCCCGTGAAGATCAGTTTTGATCCCGGATTTTCGCGCCATAAGCGAATACCTTCATTCAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCAC  >  REL606/990091‑990234

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