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,2990AT99.1% 923.9 / ‑3.1 216intergenic (+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 (0/2);  new base T (132/82);  total (132/84)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCG  >  REL606/4273226‑4273371
                                                                         |                                                                         
cgTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                         <  1:2512535/76‑1 (MQ=255)
 gTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCa                                                                        <  1:126548/76‑1 (MQ=255)
 gTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCa                                                                        <  2:4020133/76‑1 (MQ=255)
  ttattaTGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCAt                                                                       >  1:1398303/1‑76 (MQ=255)
  ttattaTGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCAt                                                                       >  2:5832295/1‑76 (MQ=255)
  ttattaTGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCAt                                                                       <  2:4774661/76‑1 (MQ=255)
  ttattaTGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCAt                                                                       >  2:3291789/1‑76 (MQ=255)
   tattatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                      >  2:6232108/1‑76 (MQ=255)
   tattatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                      <  2:5551276/76‑1 (MQ=255)
   tattatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                      >  2:2264534/1‑76 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                     >  1:5804669/1‑76 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                     <  1:947437/76‑1 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                     <  2:384224/76‑1 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                    <  1:3428524/76‑1 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                    <  2:4114561/76‑1 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                    >  2:4199202/1‑76 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                    <  2:6269709/76‑1 (MQ=255)
      tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                   >  2:2960673/1‑76 (MQ=255)
      tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                   >  1:4806033/1‑76 (MQ=255)
      tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                   >  1:1235943/1‑76 (MQ=255)
      tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                   >  1:2942218/1‑76 (MQ=255)
        tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt                                                                 >  1:6637552/1‑76 (MQ=255)
          cATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTa                                                              <  1:4818334/76‑1 (MQ=255)
          cATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTa                                                              <  2:6473345/76‑1 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          >  1:683119/1‑76 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          >  2:5269477/1‑76 (MQ=255)
                tGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaaa                                                        >  1:317981/1‑76 (MQ=255)
                  cAGCGAAAAGCTCTGTTGTTAATGGGTTGTGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATc                                                      >  2:4837982/1‑76 (MQ=255)
                   aGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCt                                                     >  1:567021/1‑76 (MQ=255)
                   aGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCt                                                     >  2:4654573/1‑76 (MQ=255)
                    gCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTa                                                    <  2:3066718/76‑1 (MQ=255)
                    gCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTa                                                    <  2:6810905/76‑1 (MQ=255)
                     cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                   <  1:3674177/76‑1 (MQ=255)
                     cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                   <  1:542835/76‑1 (MQ=255)
                      gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACa                                                  <  1:234621/76‑1 (MQ=255)
                       aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaa                                                 >  2:3687756/1‑76 (MQ=255)
                       aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaa                                                 <  1:4844732/76‑1 (MQ=255)
                        aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                                <  1:2014363/76‑1 (MQ=255)
                        aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                                >  1:608738/1‑76 (MQ=255)
                        aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                                >  1:5316038/1‑76 (MQ=255)
                         aaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAc                                               <  2:1387668/76‑1 (MQ=255)
                           gCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCa                                             <  2:167879/76‑1 (MQ=255)
                            ctctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCaa                                            <  2:5717038/76‑1 (MQ=255)
                            ctctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCaa                                            <  1:2806880/76‑1 (MQ=255)
                             tctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAAc                                           >  1:4489176/1‑76 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          <  1:71146/76‑1 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          >  1:2615468/1‑76 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          <  2:271472/76‑1 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          <  1:1850681/76‑1 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          >  1:5024191/1‑76 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          >  1:4079029/1‑76 (MQ=255)
                                   gttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAaca                                     >  1:7009691/1‑76 (MQ=255)
                                   gttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAaca                                     <  2:6670785/76‑1 (MQ=255)
                                     tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                   >  1:3089192/1‑76 (MQ=255)
                                     tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                   >  1:1279439/1‑76 (MQ=255)
                                     tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                   >  2:1779379/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAAATCTACAAACCAACGCAAcaca                                   >  2:3979078/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                  >  1:3009704/1‑76 (MQ=255)
                                       aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                 >  2:4942636/1‑76 (MQ=255)
                                         gggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTc                               <  1:4247403/76‑1 (MQ=255)
                                         gggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTc                               <  1:4152938/76‑1 (MQ=255)
                                         gggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACACAACACAATTc                               <  2:747786/76‑1 (MQ=255)
                                          ggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCa                              >  1:3129036/1‑76 (MQ=255)
                                           gTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                             >  2:2309567/1‑76 (MQ=255)
                                           gTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                             <  2:3678013/76‑1 (MQ=255)
                                           gTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                             >  2:2572252/1‑76 (MQ=255)
                                            ttGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATg                            >  2:1574570/1‑76 (MQ=255)
                                            ttGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATg                            >  2:4579022/1‑76 (MQ=255)
                                             tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGc                           >  2:6334705/1‑76 (MQ=255)
                                             tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGc                           >  2:1133405/1‑76 (MQ=255)
                                              gcgcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGcc                          >  2:2776140/1‑76 (MQ=255)
                                              gcgcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGcc                          >  1:2086292/1‑76 (MQ=255)
                                               cgcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGccc                         >  2:338968/1‑76 (MQ=255)
                                                gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATTTACAAACCAACGCAACACAATTCATGCCCt                        <  1:4408594/76‑1 (MQ=255)
                                                gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                        >  1:4847775/1‑76 (MQ=255)
                                                gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                        <  1:5771878/76‑1 (MQ=255)
                                                 cAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTg                       >  2:4881817/1‑76 (MQ=255)
                                                  aaCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTgg                      >  2:6065618/1‑76 (MQ=255)
                                                   aCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGc                     <  1:142587/76‑1 (MQ=255)
                                                    cATACCGCGCAAATGATACTGATCATAAGCGTTAAAACAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    <  1:1655213/76‑1 (MQ=255)
                                                     aTACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAg                   >  2:1517326/1‑76 (MQ=255)
                                                     aTACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAg                   >  2:6416837/1‑76 (MQ=255)
                                                      tACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGt                  >  1:1012329/1‑76 (MQ=255)
                                                       aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa                 >  1:4603306/1‑76 (MQ=255)
                                                       aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa                 >  1:457964/1‑76 (MQ=255)
                                                         cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg               <  1:6342751/76‑1 (MQ=255)
                                                         cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg               <  2:6941442/76‑1 (MQ=255)
                                                         cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg               <  2:6922133/76‑1 (MQ=255)
                                                         cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg               <  2:1080419/76‑1 (MQ=255)
                                                         cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg               <  2:1930353/76‑1 (MQ=255)
                                                          gcgcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGt              <  2:4096388/76‑1 (MQ=255)
                                                          gcgcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGt              >  2:3531885/1‑76 (MQ=255)
                                                            gcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCa            >  1:6530615/1‑76 (MQ=255)
                                                            gcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCa            >  1:3410099/1‑76 (MQ=255)
                                                             cAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTAAc           <  1:4349844/76‑1 (MQ=255)
                                                                  gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct      >  2:5759184/1‑76 (MQ=255)
                                                                  gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct      >  1:389175/1‑76 (MQ=255)
                                                                   aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc     >  1:3150715/1‑76 (MQ=255)
                                                                   aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc     >  2:5235157/1‑76 (MQ=255)
                                                                     aCTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcgc   >  1:5508141/1‑76 (MQ=255)
                                                                         |                                                                         
CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCG  >  REL606/4273226‑4273371

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