Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,962,637 T→A 9.1% intergenic (‑22/‑218) cyaY ← / → PP_5226 frataxin/Lppl family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,962,6370TA 9.1% 122.8 / 2.1 44intergenic (‑22/‑218)cyaY/PP_5226frataxin/Lppl family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (17/23);  new base A (2/2);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CACCGGCGGAGTTCTCCATGTCCAGGTCCAGGCCGCTCTCGTCGAACAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGGGCTGCGCAGAGGCTAAGCATAG  >  NC_002947/5962504‑5962764
                                                                                                                                     |                                                                                                                               
cacCGGCGGAGTTCTCCATGTCCAGGTCCAGGCCGCTCTCGTCGAACAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCAGc                                                                                                                               >  1:254889/1‑136 (MQ=255)
cacCGGCGGAGTTCTCCATGTCCAGGTCCAGGCCGCTCTCGTCGAACAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCAGc                                                                                                                               >  1:257108/1‑136 (MQ=255)
 acCGGCGGAGTTCTCCATGTCCAGGTCCAGGCCGCTCTCGTCGAACAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCg                                                                                                                              <  2:391543/136‑1 (MQ=255)
      cggAGTTCTCCATGTCCAGGTCCAGGCCGCTCTCGTCGAACATATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGATTCACTCAAACTAATTTCAGGAACCTCGAAATTGGCTGCTTTACa                                                                                                                         >  1:481964/1‑136 (MQ=255)
                      caggtccaggCCGCTCTCGTCGAACAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGa                                                                                                         >  2:351354/1‑136 (MQ=255)
                        ggtccaggCCGCTCTCGTCGAACAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATa                                                                                                       >  2:291894/1‑136 (MQ=255)
                           ccaggCCGCTCTCGTCGAACAGATCTTCCAGAGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACgg                                                                                                    >  2:343683/1‑136 (MQ=255)
                                 cGCTCTCGTCGAACAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCtt                                                                                              >  1:144181/1‑136 (MQ=255)
                                 cGCTCTCGTCGAACAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCtt                                                                                              >  1:410666/1‑136 (MQ=255)
                                     ctcGTCGAACAGATCTTCCAGAGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTgccg                                                                                          <  1:343683/136‑1 (MQ=255)
                                            aaCAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATg                                                                                   <  1:332245/136‑1 (MQ=255)
                                                 aTCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAg                                                                              >  2:475108/1‑136 (MQ=255)
                                                     tCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGttt                                                                          >  2:544937/1‑136 (MQ=255)
                                                     tCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGttt                                                                          >  2:260504/1‑136 (MQ=255)
                                                              tGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGATTAGCCGTTGc                                                                 >  1:282868/1‑136 (MQ=255)
                                                              tGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATAAGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGc                                                                 >  1:325817/1‑136 (MQ=255)
                                                                ttGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCt                                                               <  2:41758/136‑1 (MQ=255)
                                                                  gggTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTc                                                             <  2:487380/136‑1 (MQ=255)
                                                                  gggTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTc                                                             <  2:186172/136‑1 (MQ=255)
                                                                  gggTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCAATGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTc                                                             >  1:456030/1‑136 (MQ=255)
                                                                        cgTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGc                                                       >  2:361642/1‑136 (MQ=255)
                                                                            gACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCctct                                                   <  2:467216/136‑1 (MQ=255)
                                                                                aGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTcgc                                               <  1:179347/136‑1 (MQ=255)
                                                                                      tGGAAACGCGCTTGACTCAATCTCTTTGCAGGAACCTCGAATTTGGCTGCGTTACAGTCAAGCCGGGCAAGATAGGGGCGCTTGCCGCCGATTGCATTGGTTTAGCCGTTGCCTTCACCGGGCTCTTCGCGGGCaa                                         <  1:66259/136‑1 (MQ=255)
                                                                                      tGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCaa                                         <  2:5686/136‑1 (MQ=255)
                                                                                      tGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCaa                                         <  2:22703/136‑1 (MQ=255)
                                                                                      tGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCaa                                         <  2:195092/136‑1 (MQ=255)
                                                                                      tGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCaa                                         <  1:475108/136‑1 (MQ=255)
                                                                                        gAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGc                                       <  1:467772/136‑1 (MQ=255)
                                                                                        gAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGc                                       <  2:40752/136‑1 (MQ=255)
                                                                                         aaaCGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGcc                                      <  1:257810/136‑1 (MQ=255)
                                                                                               gcTTCCCTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGAGGGCAAGCCCGCAcc                                >  2:51871/1‑136 (MQ=255)
                                                                                               gcTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTcc                                >  1:198593/1‑136 (MQ=255)
                                                                                                cTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCt                               <  1:351354/136‑1 (MQ=255)
                                                                                                   cACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGGTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACa                            <  2:381513/136‑1 (MQ=255)
                                                                                                     cTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAgg                          <  2:410666/136‑1 (MQ=255)
                                                                                                     cTCAAACTCATTGCAGGAACCTCGAAATTGGCAGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAgg                          <  2:254889/136‑1 (MQ=255)
                                                                                                     cTCAAACTCATTGCAGGAACCTCGAAATTGGCAGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCATCGCGGGCAAGCCCGCTCCTACAgg                          <  2:257108/136‑1 (MQ=255)
                                                                                                            tCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGGGCTgc                   <  1:478519/136‑1 (MQ=255)
                                                                                                            tCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCGTAACGGGCAACCCCCCACCAAAAGGGGGccgg                   >  2:449948/1‑133 (MQ=255)
                                                                                                            tCATTGCAGGAACCTCGAAATTGGCTGAGTTACAGGCAAGCCGGGCAAGATACGGGCGCTAGCCGCCGAAAGCAAAGGTATAGCAGTTGCCTTCACCGGCCTCATAGAGGGCAAACACCCTACAACAAGGGGCTgc                   >  2:199263/1‑136 (MQ=255)
                                                                                                             cATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGGGCTgcg                  <  1:107458/136‑1 (MQ=255)
                                                                                                                 gCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGGGCTGCGCaga              <  1:146834/136‑1 (MQ=255)
                                                                                                                             aaaTTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGGGCTGCGCAGAGGCTAAGCATAg  <  1:233664/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
CACCGGCGGAGTTCTCCATGTCCAGGTCCAGGCCGCTCTCGTCGAACAGATCTTCCAGGGCCTGTTGGGTCGCGTCGACCAGATCATGGAAACGCGCTTCACTCAAACTCATTGCAGGAACCTCGAAATTGGCTGCGTTACAGGCAAGCCGGGCAAGATACGGGCGCTTGCCGCCGAATGCAAAGGTTTAGCCGTTGCCTTCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGGGGCTGCGCAGAGGCTAAGCATAG  >  NC_002947/5962504‑5962764

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