Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,823,561 A→C 60.6% intergenic (+306/+652) 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
*Exported2,823,5610AC60.6% 139.7 / 100.2 109intergenic (+306/+652)nrfG/pheUheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (15/28);  new base C (26/40);  total (41/68)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

TTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGC  >  Exported/2823492‑2823627
                                                                     |                                                                  
  aaaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCAc                                                                  <  1:1570985/70‑1 (MQ=255)
   aaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACgg                                                                <  1:1036278/71‑1 (MQ=255)
   aaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACgg                                                                <  1:196280/71‑1 (MQ=255)
   aaaaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACgg                                                                <  1:1953000/71‑1 (MQ=255)
     aaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACggg                                                               >  1:426767/1‑70 (MQ=255)
     aaaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACgggg                                                              <  1:1523505/71‑1 (MQ=255)
      aaaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGa                                                             <  1:459520/71‑1 (MQ=255)
      aaaTCTACAAACCAACGCAACACAATTCATGCCCAGGCAGTATGTCACGTTCTCGCGTTTCTGCACggg                                                               >  1:3275463/1‑69 (MQ=255)
       aaTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGa                                                             <  1:2489265/70‑1 (MQ=255)
        aTCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAAc                                                           >  1:2913080/1‑71 (MQ=255)
           tACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCAc                                                                  <  1:1174270/61‑1 (MQ=255)
           tACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACggg                                                               <  1:2206791/64‑1 (MQ=255)
             cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACgg                                                         <  1:3434228/68‑1 (MQ=255)
             cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACgg                                                         <  1:2557499/68‑1 (MQ=255)
             cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACgg                                                         <  1:3130077/68‑1 (MQ=255)
             cAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGgcgc                                                      <  1:3179479/71‑1 (MQ=255)
                aCCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTc                                                    >  1:473221/1‑70 (MQ=255)
                aCCAACGCAACACAATTCATGCACTGGCAGTATGTCACGTTCTCGCGTTTCTCaa                                                                   >  1:898567/1‑55 (MQ=255)
                   aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                                   >  1:3590782/1‑52 (MQ=255)
                   aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                                   >  1:1056391/1‑52 (MQ=255)
                   aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCAtt                                                <  1:1928345/71‑1 (MQ=255)
                   aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCAtt                                                <  1:2664479/71‑1 (MQ=255)
                   aaCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCAtt                                                <  1:3331176/71‑1 (MQ=255)
                    aCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                                   >  1:3149285/1‑51 (MQ=255)
                    aCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                                   >  1:916323/1‑51 (MQ=255)
                    aCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                                   >  1:1544231/1‑51 (MQ=255)
                    aCGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTg                                               <  1:580328/71‑1 (MQ=255)
                      gCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAg                                             <  1:3108039/71‑1 (MQ=255)
                      gCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACg                                                                 >  1:798001/1‑51 (MQ=255)
                        aaCACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCt                                                     >  1:1356321/1‑61 (MQ=255)
                        aaCACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCt                                                     >  1:2136624/1‑61 (MQ=255)
                        aaCACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAgg                                            >  1:3133376/1‑70 (MQ=255)
                         acacAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACgg                                                                <  1:3683298/49‑1 (MQ=255)
                           acaATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACgggg                                                              >  1:3270017/1‑49 (MQ=255)
                           acaATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACgggg                                                              >  1:985826/1‑49 (MQ=255)
                           acaATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACggg                                                               >  1:1032228/1‑48 (MQ=255)
                             aaTTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGa                                           <  1:2094606/66‑1 (MQ=255)
                             aaTTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGa                                           <  1:2002410/66‑1 (MQ=255)
                              aTTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACgg                                                                <  1:876964/44‑1 (MQ=255)
                               ttCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcatt                                    >  1:2125330/1‑71 (MQ=255)
                                tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                                   >  1:2145697/1‑39 (MQ=255)
                                tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                                   >  1:268267/1‑39 (MQ=255)
                                tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                                   >  1:687826/1‑39 (MQ=255)
                                tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGaa                                                                   >  1:1297122/1‑39 (MQ=255)
                                tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAgg                                            <  1:845450/62‑1 (MQ=255)
                                tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGATCGGGGAACGGCGCTCCATTGa                                              <  1:1304714/60‑1 (MQ=255)
                                tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGa                                              <  1:2405933/60‑1 (MQ=255)
                                tCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGa                                              <  1:2730958/60‑1 (MQ=255)
                                  aTGCCCTGGCAGTATGTCACGTTCTCGGGTTTCTGCACGGGGAACGGCGCTCCATTGAGGaa                                          >  1:2283553/1‑62 (MQ=255)
                                  aTGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGaa                                          >  1:3715896/1‑62 (MQ=255)
                                  aTGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcattcat                                 <  1:2212346/71‑1 (MQ=255)
                                    gCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                                    <  1:1862397/50‑1 (MQ=255)
                                    gCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTc                                                    <  1:1494944/50‑1 (MQ=255)
                                     cccTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataa                               >  1:1394440/1‑70 (MQ=255)
                                     cccTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataa                               >  1:1081670/1‑70 (MQ=255)
                                      ccTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcattca                                  <  1:2670769/66‑1 (MQ=255)
                                      ccTGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataa                               >  1:3281668/1‑69 (MQ=255)
                                      ccTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCataa                               >  1:533356/1‑69 (MQ=255)
                                        tGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACgg                                                         <  1:92054/41‑1 (MQ=255)
                                        tGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTg                                               >  1:393815/1‑51 (MQ=255)
                                        tGGCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTg                                               >  1:1204238/1‑51 (MQ=255)
                                         ggCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCa                                                  >  1:2431735/1‑47 (MQ=255)
                                         ggCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa                           >  1:2359772/1‑70 (MQ=255)
                                          gCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcatt                                    >  1:964902/1‑60 (MQ=255)
                                          gCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtca                                      <  1:2897445/58‑1 (MQ=255)
                                          gCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGaa                          <  1:2209736/70‑1 (MQ=255)
                                          gCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa                           <  1:736684/69‑1 (MQ=255)
                                          gCAGTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGa                           <  1:1827862/69‑1 (MQ=255)
                                          gCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGt                                        <  1:2594546/56‑1 (MQ=255)
                                             gTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtc                                       <  1:864108/54‑1 (MQ=255)
                                             gTATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtc                                       <  1:830625/54‑1 (MQ=255)
                                              tATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCa                     <  1:291537/71‑1 (MQ=255)
                                              tATGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCa                     <  1:2854270/71‑1 (MQ=255)
                                               aTGTCACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCAt                    >  1:2561140/1‑71 (MQ=255)
                                                 gTCACGTTCTCGCGTTTCTGCACGTGGAACGGCGCTCCATTGAGGAAGtc                                       <  1:2616992/50‑1 (MQ=255)
                                                   cACGTTCTCGCGTTTGTGCACGGGGAACGGCGCTCCATTGAGGAAGtcat                                     >  1:2176642/1‑50 (MQ=255)
                                                   cACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcat                                     >  1:2271518/1‑50 (MQ=255)
                                                   cACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcat                                     >  1:1073937/1‑50 (MQ=255)
                                                   cACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtcat                                     >  1:1189566/1‑50 (MQ=255)
                                                   cACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGtca                                      >  1:416504/1‑49 (MQ=255)
                                                   cACGTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAg                            <  1:2308150/59‑1 (MQ=255)
                                                    aCGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGaa               <  1:69990/71‑1 (MQ=255)
                                                     cGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGaaa              <  1:1646152/71‑1 (MQ=255)
                                                     cGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGaaa              <  1:1738819/71‑1 (MQ=255)
                                                      gTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGaaaa             <  1:151660/71‑1 (MQ=255)
                                                      gTTCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGaaaa             <  1:2682894/71‑1 (MQ=255)
                                                       ttCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAgg                                            <  1:196758/39‑1 (MQ=255)
                                                       ttCTCGCGTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAg            <  1:2463737/71‑1 (MQ=255)
                                                         ctcGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAgg           >  1:1985728/1‑70 (MQ=255)
                                                         ctcGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAgg           >  1:2573960/1‑70 (MQ=255)
                                                            gcgTTTCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAggg          <  1:2072724/68‑1 (MQ=255)
                                                              gTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAatga     >  1:1243953/1‑71 (MQ=255)
                                                              gTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAatga     >  1:1265365/1‑71 (MQ=255)
                                                              gTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAatga     >  1:2173912/1‑71 (MQ=255)
                                                                 tCTGCACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAatgatg   <  1:3364975/70‑1 (MQ=255)
                                                                  cTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAg                            >  1:842225/1‑44 (MQ=255)
                                                                  cTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGc  <  1:2663418/70‑1 (MQ=255)
                                                                  cTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGc  <  1:1964246/70‑1 (MQ=255)
                                                                   tGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCata                                <  1:2612864/39‑1 (MQ=255)
                                                                   tGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCata                                <  1:2495626/39‑1 (MQ=255)
                                                                     |                                                                  
TTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGCGTTTCTGAACGGGGAACGGCGCTCCATTGAGGAAGTCATTCATAATAGAAACTCATTCGAAAAGGGAATGATGC  >  Exported/2823492‑2823627

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: