Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,411,718 T→C 10.7% L330L (TTG→CTG)  hemY → predicted protoheme IX synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,411,7180TC10.7% 156.2 / 9.4 56L330L (TTG→CTG) hemYpredicted protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base T (34/16);  new base C (3/3);  total (37/19)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

CAAACAATCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCCGCGC  >  Exported/2411650‑2411788
                                                                    |                                                                      
caaacaaTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgttg                                                                      <  1:3513659/71‑1 (MQ=255)
  aacaaTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg                                                                    <  1:581218/71‑1 (MQ=255)
   acaaTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg                                                                    >  1:3361774/1‑70 (MQ=255)
     aaTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAg                                                                 >  1:3133828/1‑71 (MQ=255)
      aTCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAg                                                                 >  1:1880583/1‑70 (MQ=255)
        ccGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgtt                                                                       >  1:1867890/1‑62 (MQ=255)
             aCAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTgg                                                         <  1:1904804/71‑1 (MQ=255)
             aCAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTgg                                                         <  1:3073392/71‑1 (MQ=255)
              cAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTggg                                                        >  1:3605242/1‑71 (MQ=255)
                 cTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGc                                                       >  1:336028/1‑69 (MQ=255)
                 cTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCa                                                     <  1:3518169/71‑1 (MQ=255)
                  tGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgt                                                                        >  1:2646774/1‑51 (MQ=255)
                  tGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgt                                                                        >  1:1918312/1‑51 (MQ=255)
                   ggAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg                                                                    <  1:3671641/54‑1 (MQ=255)
                   ggAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg                                                                    <  1:1679766/54‑1 (MQ=255)
                    gAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTc                                                  <  1:1074281/71‑1 (MQ=255)
                     aaaaaGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCa                                                 <  1:2686118/71‑1 (MQ=255)
                       aaaGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTc                                                  >  1:912222/1‑68 (MQ=255)
                       aaaGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACt                                               >  1:474833/1‑71 (MQ=255)
                        aaGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACt                                               >  1:973108/1‑70 (MQ=255)
                         aGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACtga                                             >  1:1401185/1‑71 (MQ=255)
                            gctgcGCCAGCAAATCAAAACCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACt                                               >  1:1323870/1‑66 (MQ=255)
                            gctgcGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTtgtg                                                                    <  1:1894841/45‑1 (MQ=255)
                                cgcCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAg                                        <  1:2142666/69‑1 (MQ=255)
                                cgcCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCa                                      >  1:2050828/1‑71 (MQ=255)
                                 gcCAGCAAATCAAAAACGTCGGCGATCGCCCGCtgt                                                                        >  1:2909586/1‑36 (MQ=255)
                                 gcCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCAc                                     >  1:2623543/1‑71 (MQ=255)
                                  ccAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGcaca                                                             >  1:1156764/1‑46 (MQ=255)
                                    aGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCAc                                     >  1:2447990/1‑68 (MQ=255)
                                     gCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAg                                        >  1:3117118/1‑64 (MQ=255)
                                     gCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACg                                    >  1:348709/1‑68 (MQ=255)
                                         aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACtga                                             >  1:714841/1‑55 (MQ=255)
                                         aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACtga                                             >  1:1230856/1‑55 (MQ=255)
                                         aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATg                             >  1:2645248/1‑71 (MQ=255)
                                         aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATg                             >  1:3410132/1‑71 (MQ=255)
                                         aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATg                             >  1:659351/1‑71 (MQ=255)
                                         aTCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATg                             <  1:1714804/71‑1 (MQ=255)
                                          tCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTgg                                                         <  1:1427471/42‑1 (MQ=255)
                                          tCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTgg                                                         <  1:959998/42‑1 (MQ=255)
                                           cAAAAACGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGc                                                       <  1:3022912/43‑1 (MQ=255)
                                           cAAAAACGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGc                                                       <  1:3489391/43‑1 (MQ=255)
                                           cAAAAACGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGc                                                       <  1:2632948/43‑1 (MQ=255)
                                            aaaaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATgg                            >  1:3172759/1‑69 (MQ=255)
                                            aaaaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGc                           <  1:2556074/70‑1 (MQ=255)
                                             aaaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAACCACGGAGAATGGCAg                         >  1:2355897/1‑71 (MQ=255)
                                              aaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGgag                                 >  1:625178/1‑62 (MQ=255)
                                              aaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGgag                                 >  1:3487461/1‑62 (MQ=255)
                                              aaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAgg                        <  1:1081144/71‑1 (MQ=255)
                                              aaaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAg                         >  1:3564143/1‑70 (MQ=255)
                                               aaCGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAgg                        >  1:573777/1‑70 (MQ=255)
                                                aCGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAgg                        >  1:1395580/1‑69 (MQ=255)
                                                aCGTCGGCGATCGCCCGCTGCTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAgg                        >  1:3684795/1‑69 (MQ=255)
                                                      gCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACt                                               >  1:439275/1‑40 (MQ=255)
                                                       cGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCAtcgc               <  1:2266518/71‑1 (MQ=255)
                                                          tCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCa                   >  1:1707663/1‑64 (MQ=255)
                                                           cGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCAtcgctcgc           >  1:2727824/1‑71 (MQ=255)
                                                           cGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCAtcgctcgc           >  1:2070958/1‑71 (MQ=255)
                                                           cGCCCGCTGCTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCAtcgctcgc           >  1:1464952/1‑71 (MQ=255)
                                                             cccGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGc                                       >  1:1027414/1‑41 (MQ=255)
                                                             cccGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCt         >  1:874216/1‑71 (MQ=255)
                                                             cccGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCt         >  1:1876233/1‑71 (MQ=255)
                                                                  tgttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCcgc    >  1:1061061/1‑71 (MQ=255)
                                                                  tgttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCcgc    >  1:2945427/1‑71 (MQ=255)
                                                                    ttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGcc          <  1:1295227/63‑1 (MQ=255)
                                                                    ttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGcc          <  1:1304644/63‑1 (MQ=255)
                                                                    ttgtGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCCgcgc  <  1:1211084/71‑1 (MQ=255)
                                                                    |                                                                      
CAAACAATCCGGAACAGCTGGAAAAAGTGCTGCGCCAGCAAATCAAAAACGTCGGCGATCGCCCGCTGTTGTGGAGCACACTGGGCCAGTCACTGATGAAGCACGGAGAATGGCAGGAAGCATCGCTCGCCTTCCGCGC  >  Exported/2411650‑2411788

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: