Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,273,299 A→T intergenic (+227/‑115) nrfG → / → gltP heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,273,2990AT98.1% 1107.5 / 7.1 267intergenic (+227/‑115)nrfG/gltPheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base T (142/120);  total (145/122)
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.96e-01

Only 100 of 267 total aligned reads displayed.

CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCG  >  REL606/4273226‑4273371
                                                                         |                                                                        
cgTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                        <  2:6474698/76‑1 (MQ=255)
   tattatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                     >  1:2536908/1‑76 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATac                                                                    >  2:1573460/1‑75 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                    >  1:3529461/1‑76 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                   <  2:5442525/76‑1 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                   <  2:4671368/76‑1 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                   >  2:3968280/1‑76 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                   <  1:5394696/76‑1 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                   <  1:2755129/76‑1 (MQ=255)
      tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                  <  2:6316155/76‑1 (MQ=255)
         gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGtt                                                               >  1:1038385/1‑76 (MQ=255)
         gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGtt                                                               <  2:5572891/76‑1 (MQ=255)
           aTTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaa                                                             >  1:2588692/1‑76 (MQ=255)
            ttAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaa                                                            <  1:5703778/76‑1 (MQ=255)
            ttAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaa                                                            >  1:1673538/1‑76 (MQ=255)
            ttAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaa                                                            <  1:1350269/76‑1 (MQ=255)
             tAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaa                                                           <  2:1333649/76‑1 (MQ=255)
             tAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaa                                                           >  1:370379/1‑76 (MQ=255)
              aaTGCGGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          <  1:2027884/76‑1 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          >  2:5979755/1‑76 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          >  2:5374520/1‑76 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          >  2:3598254/1‑76 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          >  1:3467496/1‑76 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          <  2:2138788/76‑1 (MQ=255)
               aTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaa                                                         <  2:2349592/76‑1 (MQ=255)
                tGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaaa                                                        >  2:5178797/1‑76 (MQ=255)
                 gCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAAt                                                       >  2:2879351/1‑76 (MQ=255)
                  cAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATc                                                      <  1:3144341/76‑1 (MQ=255)
                   aGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCt                                                     <  2:6705442/76‑1 (MQ=255)
                     cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTAc                                                   <  1:7630764/76‑1 (MQ=255)
                      gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACa                                                  >  1:1531840/1‑76 (MQ=255)
                      gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACa                                                  <  1:6131688/76‑1 (MQ=255)
                       aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaa                                                 >  2:480969/1‑76 (MQ=255)
                       aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaa                                                 >  2:5073743/1‑76 (MQ=255)
                       aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaa                                                 <  1:2881829/76‑1 (MQ=255)
                        aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                                <  1:5730075/76‑1 (MQ=255)
                        aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACAGTTCATAAGCGTTAAAAAAATCTACaaa                                                >  2:3326035/1‑76 (MQ=255)
                          aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc                                              >  1:4592499/1‑76 (MQ=255)
                          aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc                                              <  1:7469232/76‑1 (MQ=255)
                           gCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCa                                             <  2:5815453/76‑1 (MQ=255)
                             tctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAAc                                           >  1:5877760/1‑76 (MQ=255)
                             tctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAAc                                           >  2:7071399/1‑76 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          >  1:4677070/1‑76 (MQ=255)
                                gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                        >  2:4242881/1‑76 (MQ=255)
                                gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                        <  1:4966099/76‑1 (MQ=255)
                                gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                        >  2:6135822/1‑76 (MQ=255)
                                 ttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCaa                                       <  1:7232541/76‑1 (MQ=255)
                                 ttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCaa                                       >  1:7667772/1‑76 (MQ=255)
                                   gttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAaca                                     <  1:1545226/76‑1 (MQ=255)
                                    ttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                    >  1:446874/1‑76 (MQ=255)
                                     tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                   <  1:2169708/76‑1 (MQ=255)
                                     tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                   <  1:3342592/76‑1 (MQ=255)
                                     tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAAAGCAAcaca                                   >  2:6991262/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                  >  2:6274097/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                  >  1:4146363/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                  >  2:1697302/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                  >  2:7425203/1‑76 (MQ=255)
                                       aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                 >  1:5176210/1‑76 (MQ=255)
                                       aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                 >  2:6015103/1‑76 (MQ=255)
                                       aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                 <  2:5100118/76‑1 (MQ=255)
                                          gggtGCGCAACATACCGCGCACATGATACTGTTGATAAGCGTTAAAAAAATATAGTCACCAACGTACCACACTTCa                              >  2:708514/4‑76 (MQ=255)
                                          ggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCa                              >  1:1298779/1‑76 (MQ=255)
                                          ggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCa                              >  1:3429880/1‑76 (MQ=255)
                                          ggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCa                              >  1:3874633/1‑76 (MQ=255)
                                          ggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCa                              <  1:6589134/76‑1 (MQ=255)
                                           ggtGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                             <  2:1934797/74‑1 (MQ=255)
                                           gTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                             <  2:815165/76‑1 (MQ=255)
                                             tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGc                           >  2:3653862/1‑76 (MQ=255)
                                              gcgcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGcc                          >  2:3546306/1‑76 (MQ=255)
                                                gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                        <  1:7178884/76‑1 (MQ=255)
                                                gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                        >  2:4532714/1‑76 (MQ=255)
                                                gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                        >  1:4763371/1‑76 (MQ=255)
                                                gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                        <  1:1903151/76‑1 (MQ=255)
                                                 cAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTg                       <  2:1151681/76‑1 (MQ=255)
                                                 cAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTg                       <  1:5541117/76‑1 (MQ=255)
                                                 cAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTg                       >  1:3218617/1‑76 (MQ=255)
                                                  aaCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTgg                      <  2:4950356/76‑1 (MQ=255)
                                                    cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    >  2:7564113/1‑76 (MQ=255)
                                                    cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    >  2:7288860/1‑76 (MQ=255)
                                                    cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    >  2:864020/1‑76 (MQ=255)
                                                    cATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    >  2:4572754/1‑76 (MQ=255)
                                                     aTACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAg                   >  1:84770/1‑76 (MQ=255)
                                                       aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa                 >  1:1145404/1‑76 (MQ=255)
                                                       aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa                 >  1:2816949/1‑76 (MQ=255)
                                                        ccGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt                <  2:5617845/76‑1 (MQ=255)
                                                          gcgcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGt              <  1:1163540/76‑1 (MQ=255)
                                                           cgcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTc             <  1:693397/76‑1 (MQ=255)
                                                            gcAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCa            >  2:4824448/1‑76 (MQ=255)
                                                               aaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGt         >  1:1453278/1‑76 (MQ=255)
                                                                aTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGtt        <  1:623272/76‑1 (MQ=255)
                                                                aTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGtt        <  1:3669492/76‑1 (MQ=255)
                                                                  gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct      >  2:1222954/1‑76 (MQ=255)
                                                                  gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct      >  1:6621163/1‑76 (MQ=255)
                                                                  gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct      >  1:4276357/1‑76 (MQ=255)
                                                                   aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc     >  1:6882804/1‑76 (MQ=255)
                                                                   aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc     >  2:5505766/1‑76 (MQ=255)
                                                                    tACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg    <  2:6907588/76‑1 (MQ=255)
                                                                    tACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg    >  2:2636840/1‑76 (MQ=255)
                                                                    tACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg    >  2:4196410/1‑76 (MQ=255)
                                                                     aCTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcgc   >  1:6506795/1‑76 (MQ=255)
                                                                         |                                                                        
CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCG  >  REL606/4273226‑4273371

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