Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,936,315 G→A G196G (GGC→GGT torY ← TMAO reductase III (TorYZ), cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,936,3150GA99.5% 746.0 / ‑5.7 187G196G (GGC→GGTtorYTMAO reductase III (TorYZ), cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (101/85);  total (102/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 190 total aligned reads displayed.

CGGTGAGATCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTGGC  >  REL606/1936242‑1936387
                                                                         |                                                                        
cgGTGAGATCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCAcc                                                                        >  2:4719678/1‑76 (MQ=255)
 ggTGAGATCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCa                                                                       <  2:5491294/76‑1 (MQ=255)
  gtgtGATCGGTACCAGGCTTCCCGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCAt                                                                      <  2:1069689/76‑1 (MQ=255)
  gTGAGATCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCAt                                                                      >  2:5855909/1‑76 (MQ=255)
  gTGAGATCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCAt                                                                      <  2:5131671/76‑1 (MQ=255)
   tGAGATCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCAtt                                                                     >  2:193728/1‑76 (MQ=255)
    gagaTCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTa                                                                    <  1:5616918/76‑1 (MQ=255)
     agaTCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAg                                                                   <  2:4784491/76‑1 (MQ=255)
     agaTCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAg                                                                   <  2:4251522/76‑1 (MQ=255)
       aTCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTc                                                                 >  1:73407/1‑76 (MQ=255)
       aTCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTc                                                                 <  1:251526/76‑1 (MQ=255)
        tCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCa                                                                <  2:3091776/76‑1 (MQ=255)
        tCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCa                                                                <  1:4535798/76‑1 (MQ=255)
        tCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCa                                                                <  1:6022747/76‑1 (MQ=255)
          ggTACCAGGATTCCCGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACt                                                              <  1:3211530/76‑1 (MQ=255)
          ggTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACt                                                              <  2:5497245/76‑1 (MQ=255)
          ggTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACAATTAGTCAca                                                              >  2:1434383/1‑75 (MQ=255)
           gTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTg                                                             >  2:3777028/1‑76 (MQ=255)
           gTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTg                                                             >  1:1778202/1‑76 (MQ=255)
            tACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGa                                                            >  1:3286094/1‑76 (MQ=255)
            tACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGa                                                            <  1:6229909/76‑1 (MQ=255)
            tACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGa                                                            >  1:5192632/1‑76 (MQ=255)
            tACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGa                                                            >  1:388067/1‑76 (MQ=255)
             aCCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGaa                                                           >  2:676577/1‑76 (MQ=255)
             aCCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGaa                                                           <  2:5694476/76‑1 (MQ=255)
              ccAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAg                                                          >  2:2437464/1‑76 (MQ=255)
               gaGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGt                                                         <  2:3542980/75‑1 (MQ=255)
                aGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTa                                                        <  1:4180573/76‑1 (MQ=255)
                aGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTa                                                        >  1:948007/1‑76 (MQ=255)
                      cACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTc                                                  >  2:1353460/1‑76 (MQ=255)
                       aCGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCt                                                 <  1:3355269/76‑1 (MQ=255)
                       aCGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCt                                                 >  2:5308031/1‑76 (MQ=255)
                        cgttAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTg                                                <  1:1659663/73‑1 (MQ=255)
                        cGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTg                                                <  1:2232332/76‑1 (MQ=255)
                         ggTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGa                                               >  2:3463039/1‑76 (MQ=255)
                          gTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGAc                                              >  2:5905805/1‑76 (MQ=255)
                           tAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACt                                             >  1:4086472/1‑76 (MQ=255)
                            aGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACtt                                            <  1:1015550/76‑1 (MQ=255)
                             gCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACttt                                           >  1:3955735/1‑76 (MQ=255)
                              ccAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTc                                          >  2:73960/1‑76 (MQ=255)
                               cAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTcc                                         >  2:6140522/1‑76 (MQ=255)
                               cAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTcc                                         <  2:1217312/76‑1 (MQ=255)
                               cAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTcc                                         <  1:2500066/76‑1 (MQ=255)
                                aGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCa                                        <  1:1349399/76‑1 (MQ=255)
                                aGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCa                                        <  2:1485041/76‑1 (MQ=255)
                                aGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCa                                        >  2:2133010/1‑76 (MQ=255)
                                 gCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAg                                       >  2:2606183/1‑76 (MQ=255)
                                 gCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAg                                       <  2:4561144/76‑1 (MQ=255)
                                  cTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGc                                      <  1:225538/76‑1 (MQ=255)
                                    cGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTc                                    >  1:5815287/1‑76 (MQ=255)
                                    cGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTc                                    >  2:251361/1‑76 (MQ=255)
                                    cGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTc                                    >  2:4398874/1‑76 (MQ=255)
                                    cGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTc                                    <  2:2420476/76‑1 (MQ=255)
                                     gCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCa                                   <  1:2848418/76‑1 (MQ=255)
                                     gCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCa                                   >  2:2049081/1‑76 (MQ=255)
                                      ccTATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCAt                                  >  2:1720974/1‑76 (MQ=255)
                                        tATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATgc                                <  1:3492792/76‑1 (MQ=255)
                                        tATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATgc                                >  1:1891852/1‑76 (MQ=255)
                                        tATGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATgc                                >  2:4442307/1‑76 (MQ=255)
                                         aTGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATgcg                               <  1:5274742/76‑1 (MQ=255)
                                          tGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATgcgc                              >  1:5728471/1‑76 (MQ=255)
                                          tGCGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATgcgc                              >  1:3864601/1‑76 (MQ=255)
                                            cGAGAAGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGcgccg                            >  2:1585256/1‑76 (MQ=255)
                                             gagaAGTTTTGAAAGGATAAATATGTGCGCCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATgcgcgcc                           >  2:2981247/1‑73 (MQ=255)
                                                aaGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGtt                        <  1:4646785/76‑1 (MQ=255)
                                                aaGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGtt                        <  1:4838900/76‑1 (MQ=255)
                                                aaGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGtt                        <  1:634041/76‑1 (MQ=255)
                                                aaGTTTTGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGtt                        <  2:3343147/76‑1 (MQ=255)
                                                   ttttGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTc                     >  1:5568815/1‑76 (MQ=255)
                                                   ttttGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTc                     >  1:4330104/1‑76 (MQ=255)
                                                   ttttGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTc                     >  1:4222718/1‑76 (MQ=255)
                                                    tttGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCa                    >  1:3640499/1‑76 (MQ=255)
                                                     ttGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCAt                   >  2:5757263/1‑76 (MQ=255)
                                                     ttGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCAt                   <  2:5400974/76‑1 (MQ=255)
                                                     ttGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCAt                   >  1:2239456/1‑76 (MQ=255)
                                                     ttGAAAGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCAt                   >  1:5876431/1‑76 (MQ=255)
                                                        aaaGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATcaa                >  2:3049940/1‑74 (MQ=255)
                                                         aaGGATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCAt               <  2:5589304/76‑1 (MQ=255)
                                                           ggATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCAttt             >  2:6004839/1‑76 (MQ=255)
                                                            gATAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCAtttt            >  1:145573/1‑76 (MQ=255)
                                                             aTAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCAttttt           <  1:6163062/76‑1 (MQ=255)
                                                             aTAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCAttttt           <  2:3145754/76‑1 (MQ=255)
                                                             aTAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCAttttt           <  1:4575652/76‑1 (MQ=255)
                                                             aTAAATATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCAttttt           >  2:2172996/1‑76 (MQ=255)
                                                                aaTATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTAtt        >  2:428669/1‑76 (MQ=255)
                                                                aaTATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTAtt        >  1:365217/1‑76 (MQ=255)
                                                                aaTATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTAtt        >  2:6289334/1‑76 (MQ=255)
                                                                aaTATGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTAtt        >  1:4995007/1‑76 (MQ=255)
                                                                 atatGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTAttt       >  2:6074715/1‑76 (MQ=255)
                                                                  tatGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTc      <  2:3640550/76‑1 (MQ=255)
                                                                  tatGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTc      >  2:3268016/1‑76 (MQ=255)
                                                                  tatGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTc      >  2:1128961/1‑76 (MQ=255)
                                                                   atGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCt     >  2:5365956/1‑76 (MQ=255)
                                                                   atGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCt     >  1:5436215/1‑76 (MQ=255)
                                                                   atGTGCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCt     <  1:2205944/76‑1 (MQ=255)
                                                                    tgtgCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTg    <  1:2721261/76‑1 (MQ=255)
                                                                    tgtgCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTg    >  2:5604546/1‑76 (MQ=255)
                                                                    tgtgCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTg    <  2:3333179/76‑1 (MQ=255)
                                                                      tgCACCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTGGc  >  1:1022562/1‑76 (MQ=255)
                                                                      tgCACCATTAGTCACTGAAGTAGCGGTCTAACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTGGc  >  1:200276/1‑76 (MQ=255)
                                                                         |                                                                        
CGGTGAGATCGGTACCAGGATTCACGGTAGCCAGCTCGCCTATGCGAGAAGTTTTGAAAGGATAAATATGTGCGCCATTAGTCACTGAAGTAGCGGTCTGACTTTCCAGCTCATGCGCCGCGTTGTCATCCATTTTTATTTCTGGC  >  REL606/1936242‑1936387

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