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,2990AT96.3% 213.8 / ‑1.1 54intergenic (+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 (18/34);  total (18/36)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGC  >  REL606/4273227‑4273370
                                                                        |                                                                       
gTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCa                                                                      <  1:1436231/76‑1 (MQ=255)
gTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCa                                                                      <  1:1385366/76‑1 (MQ=255)
  tattatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                    <  2:676898/76‑1 (MQ=255)
   attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                   <  1:1260007/76‑1 (MQ=255)
     tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                 >  1:228051/1‑76 (MQ=255)
     tatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGc                                                                 >  1:494692/1‑76 (MQ=255)
       tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt                                                               <  1:341313/76‑1 (MQ=255)
       tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt                                                               <  1:760388/76‑1 (MQ=255)
       tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt                                                               <  1:1301484/76‑1 (MQ=255)
        gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGtt                                                              <  2:1069699/76‑1 (MQ=255)
        gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGtt                                                              <  2:619962/76‑1 (MQ=255)
          aTTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaa                                                            >  1:206338/1‑76 (MQ=255)
          aTTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaa                                                            <  1:948761/76‑1 (MQ=255)
            tAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaa                                                          <  2:154609/76‑1 (MQ=255)
            tAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaa                                                          <  2:39103/76‑1 (MQ=255)
              aTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaan                                                        <  1:235457/76‑2 (MQ=255)
               tGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaaa                                                       >  1:1288162/1‑76 (MQ=255)
                gCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAAt                                                      <  2:823731/76‑1 (MQ=255)
                    cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                  <  2:1294737/76‑1 (MQ=255)
                    cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                  <  1:254643/76‑1 (MQ=255)
                    cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                  <  1:39070/76‑1 (MQ=255)
                    cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                  <  1:1124149/76‑1 (MQ=255)
                    cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                  <  1:1124132/76‑1 (MQ=255)
                     gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACa                                                 <  1:1034404/76‑1 (MQ=255)
                      aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaa                                                <  1:65974/76‑1 (MQ=255)
                       aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                               <  2:1021286/76‑1 (MQ=255)
                         aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc                                             <  1:1102898/76‑1 (MQ=255)
                         aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc                                             >  2:1069844/1‑76 (MQ=255)
                             ctGTTGTTAAAGGGTTGCGTAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                         >  2:620734/1‑76 (MQ=255)
                             ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                         >  2:274723/1‑76 (MQ=255)
                              tgttgtTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGc                                        >  1:89654/1‑76 (MQ=255)
                              tgttgtTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGc                                        <  1:1042770/76‑1 (MQ=255)
                               gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                       <  1:841691/76‑1 (MQ=255)
                                  gttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAaca                                    <  1:295795/76‑1 (MQ=255)
                                   ttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                   >  1:1471951/1‑76 (MQ=255)
                                   ttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                   >  2:1074124/1‑76 (MQ=255)
                                   ttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAacac                                   >  2:1146727/1‑76 (MQ=255)
                                    tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                  <  2:1236839/76‑1 (MQ=255)
                                    tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                  <  1:499345/76‑1 (MQ=255)
                                     aaaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACaa                                 >  2:1313041/1‑76 (MQ=255)
                                      aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                >  2:217767/1‑76 (MQ=255)
                                       aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                               >  2:1243467/1‑76 (MQ=255)
                                       aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                               <  2:266014/76‑1 (MQ=255)
                                         ggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCa                             >  2:922825/1‑76 (MQ=255)
                                          gTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                            <  1:409702/76‑1 (MQ=255)
                                                  aCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGc                    <  2:1446422/76‑1 (MQ=255)
                                                    aTACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAg                  <  2:903668/76‑1 (MQ=255)
                                                     tACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGt                 <  2:222785/76‑1 (MQ=255)
                                                      aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa                >  1:362006/1‑76 (MQ=255)
                                                      aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa                <  1:1157853/76‑1 (MQ=255)
                                                            cAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCAc          >  1:470714/1‑76 (MQ=255)
                                                             aaaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACg         <  1:1285968/76‑1 (MQ=255)
                                                                  aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc    <  2:1066368/76‑1 (MQ=255)
                                                                    aCTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcgc  >  1:964224/1‑76 (MQ=255)
                                                                        |                                                                       
GTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGC  >  REL606/4273227‑4273370

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