Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,740,354 G→A G326G (GGC→GGT hycE ← hydrogenase 3, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,740,3540GA98.2% 207.6 / ‑4.2 57G326G (GGC→GGThycEhydrogenase 3, large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (35/21);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GTTGCGCCAGTTGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGGCCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGGAAGAACTGC  >  REL606/2740281‑2740426
                                                                         |                                                                        
gTTGCGCCAGTTGACGGGTCTGGATCATTTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGAcc                                                                        <  2:800790/76‑1 (MQ=255)
    cgcCAGTTGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAg                                                                    >  1:1024494/1‑76 (MQ=255)
       cAGTTGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGtt                                                                 >  1:1514672/1‑76 (MQ=255)
           tggcGGGTCTGGATCTTGTCATCTTCCAGCAGACCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAc                                                             <  2:1029079/73‑1 (MQ=37)
           tGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAc                                                             >  1:1373067/1‑76 (MQ=255)
            gACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAcg                                                            >  1:1219284/1‑76 (MQ=255)
             aCGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAcgc                                                           >  1:717448/1‑76 (MQ=255)
              cGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTAcgcg                                                          >  1:1180987/1‑76 (MQ=255)
               gggTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACgcgc                                                         >  1:1168078/1‑76 (MQ=255)
                 gTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCAc                                                       >  1:1150484/1‑76 (MQ=255)
                 gTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCAc                                                       >  1:1215315/1‑76 (MQ=255)
                 gTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCAc                                                       >  1:1517414/1‑76 (MQ=255)
                 gTCTGGATCATGTCATCTTTCAGCAGACCGCGACAAACCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCAc                                                       <  1:143877/76‑1 (MQ=37)
                    tGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCgg                                                    >  1:1852017/1‑76 (MQ=255)
                     ggATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGt                                                   >  1:1353845/1‑76 (MQ=255)
                       aTCATTTCATCTTTCAGCAGATCGCGACAAACCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTaa                                                 <  1:1036373/76‑1 (MQ=25)
                       aTCATGTCATCTTTCAGCAGACCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTaa                                                 <  2:388984/76‑1 (MQ=38)
                         cATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAgg                                               >  2:1632731/1‑76 (MQ=255)
                           tGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGAt                                             <  1:687492/76‑1 (MQ=255)
                            gTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGAtc                                            >  1:1373715/1‑76 (MQ=255)
                             tCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGAtct                                           >  2:4676/1‑76 (MQ=255)
                               aTCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATctct                                         >  1:1084330/1‑76 (MQ=255)
                               aTCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATctct                                         >  1:255939/1‑76 (MQ=255)
                               aTCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATctct                                         >  1:329539/1‑76 (MQ=255)
                                tCTTTCAGCAGAGCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTg                                        <  2:78211/76‑1 (MQ=255)
                                 cTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGc                                       >  2:530282/1‑76 (MQ=255)
                                  tttCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGcc                                      >  2:1479751/1‑76 (MQ=255)
                                  tttCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGcc                                      >  1:613074/1‑76 (MQ=255)
                                   tccagcagACCGCGACGAACCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCa                                     <  1:1245341/74‑1 (MQ=37)
                                      agcagATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCAttt                                  <  2:1169844/76‑1 (MQ=255)
                                      agcagATCGCGACGAATCCCGCCGACCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCAttt                                  >  2:862778/1‑76 (MQ=255)
                                          gATCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCAt                              >  2:1219819/1‑76 (MQ=255)
                                           aTCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATg                             >  1:1604403/1‑76 (MQ=255)
                                            tCGCGACGAATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATgg                            >  1:1824277/1‑76 (MQ=255)
                                                gacgaATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGt                        <  1:296578/76‑1 (MQ=255)
                                                 ccgaACCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGtt                       <  1:1459331/75‑1 (MQ=255)
                                                   gaATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTc                     >  2:283393/1‑76 (MQ=255)
                                                   gaATCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTc                     >  1:595918/1‑76 (MQ=255)
                                                    aaTCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTca                    <  1:809834/76‑1 (MQ=255)
                                                    aaTCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTca                    >  1:1813570/1‑76 (MQ=255)
                                                    aaTCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTca                    <  2:1367437/76‑1 (MQ=255)
                                                     aTCCCGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcac                   >  2:1863784/1‑76 (MQ=255)
                                                       cccGCCGATCAGGTTCAGACCGTAGTTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgc                 >  1:1775420/1‑76 (MQ=255)
                                                       cccGCCGATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgc                 <  1:284990/76‑1 (MQ=255)
                                                       cccGCCGATCAGGTTCAGACAGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgc                 >  2:1523530/1‑76 (MQ=255)
                                                         cgccgATCAGGTTCAGACCGTAGGTTTTATGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgcac               <  1:1859415/76‑1 (MQ=255)
                                                         cgccgATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTcacgcac               <  1:1341269/76‑1 (MQ=255)
                                                            cgATCAGGTTCAGGCCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCAcgcg            <  1:985968/76‑1 (MQ=255)
                                                            cgATCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCAcgcg            >  1:1325189/1‑76 (MQ=255)
                                                              aTCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGga          >  2:578347/1‑76 (MQ=255)
                                                              aTCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGga          <  2:903647/76‑1 (MQ=255)
                                                              aTCAGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGga          <  1:1648835/76‑1 (MQ=255)
                                                                 aGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGgaaga       <  2:262793/76‑1 (MQ=255)
                                                                 aGGTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGgaaga       >  2:1815502/1‑76 (MQ=255)
                                                                  tgTTCAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGgaagaa      <  2:509803/75‑1 (MQ=255)
                                                                  ggTTCAGCCCGTCGGTTTTACGCGCCCCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGgaagaa      <  1:944333/76‑1 (MQ=255)
                                                                  ggTTCAGACCGTAGGTTTTGCGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCggaggac      >  1:593141/1‑75 (MQ=38)
                                                                      cAGACCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGGAAGAACTGc  <  1:505527/76‑1 (MQ=255)
                                                                         |                                                                        
GTTGCGCCAGTTGACGGGTCTGGATCATGTCATCTTTCAGCAGATCGCGACGAATCCCGCCGATCAGGTTCAGGCCGTAGGTTTTACGCGCACCGGTAAGGATCTCTGCCATTTTCATGGAGGTTTCACGCACGCGGAAGAACTGC  >  REL606/2740281‑2740426

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