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,2990AT97.7% 509.3 / 0.8 129intergenic (+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 (2/1);  new base T (71/55);  total (73/56)
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.40e-01

Only 100 of 129 total aligned reads displayed.

CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCG  >  REL606/4273226‑4273371
                                                                         |                                                                        
cgTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                        >  2:1715190/1‑76 (MQ=255)
  ttattaTGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCAt                                                                      <  1:1107961/76‑1 (MQ=255)
   tattatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                     >  1:1775858/1‑76 (MQ=255)
   tattatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                     >  1:283669/1‑76 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                    >  1:148691/1‑76 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                    >  1:1238805/1‑76 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                    >  1:1943023/1‑76 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                    >  1:218402/1‑76 (MQ=255)
    attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATaa                                                                    >  1:1808491/1‑76 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                   <  2:1717289/76‑1 (MQ=255)
     ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                   <  1:1446729/76‑1 (MQ=255)
      tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                  <  2:710142/76‑1 (MQ=255)
      tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                  <  1:724457/76‑1 (MQ=255)
      tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                  >  1:1291838/1‑76 (MQ=255)
        tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt                                                                >  1:250039/1‑76 (MQ=255)
        tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt                                                                <  2:319393/76‑1 (MQ=255)
         gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGtt                                                               <  1:191019/76‑1 (MQ=255)
           aTTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaa                                                             <  2:1173471/76‑1 (MQ=255)
           aTTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaa                                                             <  2:219834/76‑1 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          <  2:893982/76‑1 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          >  1:1645111/1‑76 (MQ=255)
              aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                          >  1:1445141/1‑76 (MQ=255)
               aTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaa                                                         >  2:1555317/1‑76 (MQ=255)
                tGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaaa                                                        <  2:1437841/76‑1 (MQ=255)
                  cAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATc                                                      >  1:639849/1‑76 (MQ=255)
                   aGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCt                                                     <  2:538300/76‑1 (MQ=255)
                   aGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCt                                                     <  1:395786/76‑1 (MQ=255)
                   aGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCt                                                     >  1:1144137/1‑76 (MQ=255)
                    gCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTa                                                    <  2:1400366/76‑1 (MQ=255)
                    gCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTa                                                    <  2:1936010/76‑1 (MQ=255)
                      gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                   >  1:1544588/1‑75 (MQ=255)
                      gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                   >  1:332258/1‑75 (MQ=255)
                      gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                   <  2:1544588/75‑1 (MQ=255)
                      gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                   <  2:332258/75‑1 (MQ=255)
                        aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                                >  1:852008/1‑76 (MQ=255)
                        aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACaaa                                                >  1:524825/1‑76 (MQ=255)
                         aaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAc                                               <  2:240725/76‑1 (MQ=255)
                         aaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAc                                               >  2:125523/1‑76 (MQ=255)
                          aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc                                              <  1:1754588/76‑1 (MQ=255)
                           gCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCa                                             <  1:922131/76‑1 (MQ=255)
                            ctctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCaa                                            >  2:38522/1‑76 (MQ=255)
                            ctctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCaa                                            >  2:1254192/1‑76 (MQ=255)
                             tctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAAc                                           >  2:1333525/1‑76 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          >  1:1172449/1‑76 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          >  1:777950/1‑76 (MQ=255)
                              ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                          <  1:1380298/76‑1 (MQ=255)
                                gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                        <  2:1334934/76‑1 (MQ=255)
                                gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                        <  2:1692640/76‑1 (MQ=255)
                                gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                        >  2:1497747/1‑76 (MQ=255)
                                gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                        >  2:406035/1‑76 (MQ=255)
                                gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                        <  2:585718/76‑1 (MQ=255)
                                gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                        >  1:1543271/1‑76 (MQ=255)
                                 ttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCaa                                       >  1:1263765/1‑76 (MQ=255)
                                  tgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAac                                      >  1:631451/1‑76 (MQ=255)
                                     tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                   >  1:1455015/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                  >  2:241122/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                  >  2:426586/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                  >  1:442319/1‑76 (MQ=255)
                                      aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                  >  2:784481/1‑76 (MQ=255)
                                       aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                 >  2:338772/1‑76 (MQ=255)
                                       aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                 >  2:1711237/1‑76 (MQ=255)
                                        aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                                >  2:1354723/1‑76 (MQ=255)
                                        aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                                <  2:503967/76‑1 (MQ=255)
                                          ggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCa                              >  2:1195983/1‑76 (MQ=255)
                                            ttGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATg                            >  2:1409668/1‑76 (MQ=255)
                                            ttGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATg                            >  1:574040/1‑76 (MQ=255)
                                             tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGc                           <  1:1100454/76‑1 (MQ=255)
                                             tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGc                           <  1:1100467/76‑1 (MQ=255)
                                             tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGc                           >  2:1951602/1‑76 (MQ=255)
                                             tGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGc                           <  2:42389/76‑1 (MQ=255)
                                              gcgcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGcc                          <  1:979304/76‑1 (MQ=255)
                                                gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTGAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                        >  2:1800174/1‑76 (MQ=255)
                                                 cAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTg                       >  2:24598/1‑76 (MQ=255)
                                                 cAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTg                       <  2:1905306/76‑1 (MQ=255)
                                                 cAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTg                       <  2:1419980/76‑1 (MQ=255)
                                                    cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    <  1:1366458/76‑1 (MQ=255)
                                                    cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    >  1:1371789/1‑76 (MQ=255)
                                                    cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    >  1:1122895/1‑76 (MQ=255)
                                                    cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    <  2:668506/76‑1 (MQ=255)
                                                    cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCa                    >  1:1168471/1‑76 (MQ=255)
                                                     aTACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAg                   >  1:1557584/1‑76 (MQ=255)
                                                       aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa                 >  2:178140/1‑76 (MQ=255)
                                                        ccGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt                <  2:1280074/76‑1 (MQ=255)
                                                        ccGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt                <  2:1280046/76‑1 (MQ=255)
                                                         cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg               <  1:1765386/76‑1 (MQ=255)
                                                         cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg               <  1:1316691/76‑1 (MQ=255)
                                                            gcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCa            <  2:165470/76‑1 (MQ=255)
                                                             cAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCAc           <  2:495721/76‑1 (MQ=255)
                                                              aaaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACg          >  1:893421/1‑76 (MQ=255)
                                                               aaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGt         >  1:680472/1‑76 (MQ=255)
                                                               aaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGt         >  1:486755/1‑76 (MQ=255)
                                                                aTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGtt        >  1:195330/1‑76 (MQ=255)
                                                                aTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGtt        <  1:632633/76‑1 (MQ=255)
                                                                aTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGtt        >  1:681218/1‑76 (MQ=255)
                                                                  gATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtct      <  2:854221/76‑1 (MQ=255)
                                                                   aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc     >  2:1480950/1‑76 (MQ=255)
                                                                   aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc     <  2:1578591/76‑1 (MQ=255)
                                                                   aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc     <  2:186686/76‑1 (MQ=255)
                                                                    tACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg    >  2:353812/1‑76 (MQ=255)
                                                                      cTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcgcg  >  1:736273/1‑76 (MQ=255)
                                                                         |                                                                        
CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCG  >  REL606/4273226‑4273371

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