Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,823,628 G→A 12.6% intergenic (+373/+585) nrfG → / ← pheU heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,823,6280GA12.6% 289.9 / 26.3 103intergenic (+373/+585)nrfG/pheUheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base G (61/29);  new base A (6/7);  total (67/36)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 110 total aligned reads displayed.

CGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATG  >  minE/2823563‑2823694
                                                                 |                                                                  
cGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCgatgat                                                               <  1:3150336/71‑1 (MQ=255)
        ggCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCg                                                                    <  1:467427/58‑1 (MQ=255)
          cgcTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCcac                                                      >  1:3714782/1‑70 (MQ=255)
          cgcTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCcac                                                      >  1:2420438/1‑70 (MQ=255)
              ccATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAAcc                                                  >  1:3938496/1‑70 (MQ=255)
               cATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCg                                                                    <  1:2255143/51‑1 (MQ=255)
               cATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTa                                                >  1:2394670/1‑71 (MQ=255)
                 ttGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATaa                                                             >  1:54233/1‑56 (MQ=255)
                 ttGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAtt                                              >  1:2463536/1‑71 (MQ=255)
                 ttGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAtt                                              >  1:853448/1‑71 (MQ=255)
                 ttGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAtt                                              >  1:4075233/1‑71 (MQ=255)
                  tGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAttt                                             >  1:230719/1‑71 (MQ=255)
                  tGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAtt                                              >  1:1520810/1‑70 (MQ=255)
                   gAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATaa                                                             >  1:812392/1‑54 (MQ=255)
                   gAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAtttt                                            >  1:2013470/1‑71 (MQ=255)
                   gAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAttt                                             >  1:2165621/1‑70 (MQ=255)
                    aGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAttttt                                           <  1:3830551/71‑1 (MQ=255)
                      gAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTAc                                         >  1:1926336/1‑71 (MQ=255)
                         gTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCAt                                      <  1:2640135/71‑1 (MQ=255)
                          tcattcatAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCAt                                      >  1:2241750/1‑70 (MQ=255)
                           cattcatAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAt                                               >  1:913245/1‑60 (MQ=255)
                             ttcatAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAAtt                                                           <  1:2092902/46‑1 (MQ=255)
                               catAATAGAAACTCATTCGAAAAGGGAATGATGCgat                                                                  >  1:1155329/1‑37 (MQ=255)
                               catAATAGAAACTCATTCGAAAAGGGAATGATGCgat                                                                  >  1:2899427/1‑37 (MQ=255)
                               catAATAGAAACTCATTCGAAAAGGGAATGATGCga                                                                   >  1:4040529/1‑36 (MQ=255)
                               catAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 <  1:1685532/70‑1 (MQ=255)
                               catAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 <  1:3386912/70‑1 (MQ=255)
                               catAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 <  1:623330/70‑1 (MQ=255)
                               catAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 <  1:1392420/70‑1 (MQ=255)
                                ataataGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCta                                   >  1:4006042/1‑67 (MQ=255)
                                ataataGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAg                                <  1:3692205/70‑1 (MQ=255)
                                 taataGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 <  1:3215030/68‑1 (MQ=255)
                                  aataGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATg                             >  1:1247454/1‑71 (MQ=255)
                                   ataGAAACTCATTCGAAAAGGGAATGATGCgatgat                                                               <  1:1624188/36‑1 (MQ=255)
                                   ataGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGa                               >  1:2743956/1‑68 (MQ=255)
                                   ataGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATgg                            >  1:1322619/1‑71 (MQ=255)
                                   ataGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATgg                            >  1:2132685/1‑71 (MQ=255)
                                     aGAAACTCATTCGAAAAGGGAATGATGCGATGATaa                                                             >  1:2838755/1‑36 (MQ=255)
                                     aGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAtttt                                            >  1:199414/1‑53 (MQ=255)
                                     aGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCta                                   >  1:3073178/1‑62 (MQ=255)
                                     aGAAACTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGt                          >  1:2762139/1‑71 (MQ=255)
                                      gAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAt                                               >  1:2533659/1‑49 (MQ=255)
                                      gAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGa                               >  1:2829930/1‑65 (MQ=255)
                                       aaaCTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTa                                          <  1:3595352/53‑1 (MQ=255)
                                       aaaCTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGAt                              >  1:453655/1‑65 (MQ=255)
                                       aaaCTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGttt                        <  1:3300544/71‑1 (MQ=255)
                                        aaCTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGttta                       <  1:2680768/71‑1 (MQ=255)
                                         aCTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTAcc                                        >  1:2285278/1‑53 (MQ=255)
                                         aCTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTAATTTTACCATCTATAGATGGGtttat                      >  1:3574128/1‑71 (MQ=255)
                                          cTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 <  1:4126313/59‑1 (MQ=255)
                                          cTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGtttatt                     <  1:3113606/71‑1 (MQ=255)
                                          cTCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATc                                     <  1:2603234/55‑1 (MQ=255)
                                           tCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGAt                              <  1:2912635/61‑1 (MQ=255)
                                           tCATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGtttatt                     >  1:3366720/1‑70 (MQ=255)
                                            cATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 >  1:260142/1‑57 (MQ=255)
                                            cATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGtttattta                   >  1:2067535/1‑71 (MQ=255)
                                            cATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGtttatt                     >  1:2206066/1‑69 (MQ=255)
                                            cATTCGAAAAGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGtttattta                   >  1:1872913/1‑71 (MQ=255)
                                             aTTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTAc                                         >  1:3173241/1‑48 (MQ=255)
                                             aTTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCAt                                      >  1:4064047/1‑51 (MQ=255)
                                             aTTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 <  1:3840847/56‑1 (MQ=255)
                                             aTTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 <  1:2609374/56‑1 (MQ=255)
                                             aTTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGtttattta                   <  1:637196/70‑1 (MQ=255)
                                             aTTCGAAAAGGGAATGATGCAATGATAATTGCCACAAcc                                                  <  1:575915/39‑1 (MQ=255)
                                              ttCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACa                 >  1:3196579/1‑71 (MQ=255)
                                              ttCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACa                 >  1:2074747/1‑71 (MQ=255)
                                                cGAAAAGGGAATGATGCGATGATAATTGCCACAACCt                                                 >  1:2369353/1‑37 (MQ=255)
                                                cGAAAAGGGAATGATGCGATGATAATTGCCACAACCTa                                                >  1:507456/1‑38 (MQ=255)
                                                cGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAg                                >  1:684721/1‑54 (MQ=255)
                                                cGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACata               >  1:3879769/1‑71 (MQ=255)
                                                cGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACata               >  1:3570980/1‑71 (MQ=255)
                                                cGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACata               >  1:498682/1‑71 (MQ=255)
                                                 gAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACatat              <  1:2422398/71‑1 (MQ=255)
                                                  aaaaGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATc                                     <  1:1612339/47‑1 (MQ=255)
                                                   aaaGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACAtatta            >  1:903100/1‑71 (MQ=255)
                                                   aaaGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACAtatt             >  1:2392906/1‑70 (MQ=255)
                                                   aaaGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACAtatt             >  1:2397561/1‑70 (MQ=255)
                                                    aaGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGttt                        <  1:1354094/58‑1 (MQ=255)
                                                     aGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTATATattatt          <  1:905151/71‑1 (MQ=255)
                                                     aGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATattatt          <  1:2692151/71‑1 (MQ=255)
                                                     aGGGAATGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATattatt          <  1:1890416/71‑1 (MQ=255)
                                                      gggAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTg         <  1:1482875/71‑1 (MQ=255)
                                                       ggAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTgg        >  1:115361/1‑71 (MQ=255)
                                                         aaTGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACat                >  1:1075386/1‑61 (MQ=255)
                                                         aaTGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACat                >  1:1493213/1‑61 (MQ=255)
                                                         aaTGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTg      <  1:2570664/71‑1 (MQ=255)
                                                         aaTGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTg      <  1:1478107/71‑1 (MQ=255)
                                                         aaTGATGCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTg      <  1:1880529/71‑1 (MQ=255)
                                                          atgatgCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGtttattta                   >  1:3793422/1‑57 (MQ=255)
                                                           tgatgCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGtttatt                     <  1:1116996/54‑1 (MQ=255)
                                                            gatgCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGt                          >  1:2689650/1‑48 (MQ=255)
                                                            gatgCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGt       <  1:1822740/67‑1 (MQ=255)
                                                             atgCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 >  1:3353504/1‑40 (MQ=255)
                                                             atgCGATGATAATTGCCACAACCTATTTTTACCATCtata                                 >  1:1020016/1‑40 (MQ=255)
                                                             atgCGATGATAATTGCCACAACCTATTTTTACCATCTATAGa                               >  1:1969754/1‑42 (MQ=255)
                                                             atgCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACAtatta            >  1:2861199/1‑61 (MQ=255)
                                                             atgCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATg  >  1:3884283/1‑71 (MQ=255)
                                                              tgCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATattatt          >  1:2422957/1‑62 (MQ=255)
                                                              tgCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGaa    >  1:1324286/1‑68 (MQ=255)
                                                              tgCAATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATg  >  1:2195709/1‑70 (MQ=255)
                                                                 |                                                                  
CGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGCGATGATAATTGCCACAACCTATTTTTACCATCTATAGATGGGTTTATTTACATATTATTGGTGAATG  >  minE/2823563‑2823694

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