Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,692,667 T→C 12.5% K103K (AAA→AAG iscU ← scaffold protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,692,6670TC12.5% 132.2 / 10.0 48K103K (AAA→AAGiscUscaffold protein
Reads supporting (aligned to +/- strand):  ref base T (15/27);  new base C (2/4);  total (17/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCGAGAG  >  Exported/1692598‑1692734
                                                                     |                                                                   
ttCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                                    >  1:199710/1‑71 (MQ=255)
  cACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                                    <  1:955476/69‑1 (MQ=255)
  cACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                                    <  1:6758/69‑1 (MQ=255)
  cACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                                    <  1:2710043/69‑1 (MQ=255)
   aCGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCa                                                                 >  1:1347495/1‑71 (MQ=255)
      tttGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACcg                                                              >  1:1702770/1‑71 (MQ=255)
      tttGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACcg                                                              <  1:2410453/71‑1 (MQ=255)
        tGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACcggc                                                            >  1:3355487/1‑71 (MQ=255)
                     cGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGtt                                               >  1:3424777/1‑71 (MQ=255)
                     cGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGtt                                               >  1:1497342/1‑71 (MQ=255)
                     cGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGt                                                >  1:2907000/1‑70 (MQ=255)
                        aaTGGCGGCTTTGATCGCGTCTTCTGCCAGATTAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGttctt                                            <  1:3380813/71‑1 (MQ=255)
                        aaTGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGttctt                                            >  1:1908823/1‑71 (MQ=255)
                         aTGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAtttt                                                                   >  1:1544570/1‑47 (MQ=255)
                         aTGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCa                                                   <  1:276956/63‑1 (MQ=255)
                           ggcggcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAg                                         >  1:3100382/1‑71 (MQ=255)
                              ggcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                                      <  1:2096005/71‑1 (MQ=255)
                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAAttt                                                                    <  1:1166503/40‑1 (MQ=255)
                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTGAAGTTCTTCAGCaa                                      <  1:2495944/70‑1 (MQ=255)
                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                                      <  1:1854766/70‑1 (MQ=255)
                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAat                                     <  1:2015791/71‑1 (MQ=255)
                               gcTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAat                                     <  1:2523917/71‑1 (MQ=255)
                                cTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAat                                     <  1:3324123/70‑1 (MQ=255)
                                  ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                            <  1:1873922/45‑1 (MQ=255)
                                  ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                            <  1:2998770/45‑1 (MQ=255)
                                  ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                            <  1:1181764/45‑1 (MQ=255)
                                  ttGATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggc                                                            <  1:1043055/45‑1 (MQ=255)
                                    gATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAata                                    <  1:1494418/67‑1 (MQ=255)
                                    gATCGCGTCTTCGGCCAGGATCGAGCAGTGAATCTTGACcggcgg                                                          >  1:1606743/1‑45 (MQ=255)
                                        gcgTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGtt                            <  1:846428/71‑1 (MQ=255)
                                         cgTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGttt                           <  1:3058147/71‑1 (MQ=255)
                                             ttctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGc                                        <  1:366366/54‑1 (MQ=255)
                                             ttctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTg                        <  1:2511449/70‑1 (MQ=255)
                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAg                                                 >  1:2668090/1‑43 (MQ=255)
                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAg                                         <  1:836346/51‑1 (MQ=255)
                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                                      >  1:3256441/1‑54 (MQ=255)
                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCaa                                      >  1:1956211/1‑54 (MQ=255)
                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTg                        <  1:2089608/68‑1 (MQ=255)
                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTg                        <  1:1436081/68‑1 (MQ=255)
                                               ctGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGa                       <  1:1455435/69‑1 (MQ=255)
                                               cggCCAGGATCGAGCAGTGAATCTTGACCGGCGGCAGGGCCAGttc                                              >  1:2866152/3‑46 (MQ=255)
                                                  ccAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc                     <  1:2870129/68‑1 (MQ=255)
                                                     gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc                     <  1:1315728/65‑1 (MQ=255)
                                                     gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc                     <  1:2593588/65‑1 (MQ=255)
                                                     gAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATc                     <  1:1075108/65‑1 (MQ=255)
                                                       aTAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTgcgc             >  1:1221539/1‑71 (MQ=255)
                                                           aaCAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTgc               <  1:1600955/65‑1 (MQ=255)
                                                             cAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTtcgtc       <  1:3302325/71‑1 (MQ=255)
                                                             cAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTtcgtc       >  1:1116636/1‑71 (MQ=255)
                                                               gTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTtcgtc       >  1:2885693/1‑69 (MQ=255)
                                                                  aaTTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGcc                  >  1:280225/1‑55 (MQ=255)
                                                                  aaTTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCgaga   <  1:2815253/70‑1 (MQ=255)
                                                                   aTTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAat                                     <  1:1034173/35‑1 (MQ=255)
                                                                   aTTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCgagag  <  1:2215765/70‑1 (MQ=255)
                                                                     |                                                                   
TTCACGTTTGCTTTTATAGTCCGCAATGGCGGCTTTGATCGCGTCTTCTGCCAGAATAGAACAGTGAATTTTCACCGGCGGCAGTTCAAGTTCTTCAGCAATATCGGTGTTTTTGATCGCCTGCGCTTCGTCGAGAG  >  Exported/1692598‑1692734

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