Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	980,245	G→A	5.5%	F20F (TTC→TTT)	cspI ←	cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	980,245	0	G	A	5.5%	236.0 / 3.0	72	F20F (TTC→TTT)	cspI	cold shock protein
Reads supporting (aligned to +/- strand):  ref base G (30/38);  new base A (2/2);  total (32/40)
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.43e-01

TTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTAT  >  Exported/980176‑980310                                                                      |                                                                  tttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatg                                                                   >  1:2412359/1‑70 (MQ=255)  ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                  <  1:1671398/70‑1 (MQ=255)  ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                  <  1:1110907/70‑1 (MQ=255)  ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                  <  1:1787966/70‑1 (MQ=255)  ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                  <  1:539108/70‑1 (MQ=255)  ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                  <  1:1574133/70‑1 (MQ=255)  ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCCTTCGGCGtgatga                                                                  <  1:2636730/70‑1 (MQ=255)  ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCCTTCGGCGtgatga                                                                  <  1:2632424/70‑1 (MQ=255)     aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                                                               >  1:82546/1‑70 (MQ=255)     aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                                                               >  1:2938397/1‑70 (MQ=255)     aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                                                               >  1:3358330/1‑70 (MQ=255)      aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                              >  1:588817/1‑70 (MQ=255)      aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                             <  1:3354565/71‑1 (MQ=255)      aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                             <  1:3296875/71‑1 (MQ=255)      aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                             <  1:3600398/71‑1 (MQ=255)      aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                             <  1:39410/71‑1 (MQ=255)       aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                                                               >  1:35912/1‑68 (MQ=255)       aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                                                               >  1:2756414/1‑68 (MQ=255)       aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                            <  1:3532434/71‑1 (MQ=255)       aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGAGTGATGAAAc                                                               >  1:2949511/1‑68 (MQ=255)           ttGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                        <  1:1195259/71‑1 (MQ=255)           ttGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                        <  1:2986436/71‑1 (MQ=255)           ttGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                        <  1:2912701/71‑1 (MQ=255)             gCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                                                      >  1:2256026/1‑71 (MQ=255)               tctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                     <  1:1807807/70‑1 (MQ=255)                 tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                  <  1:2438203/55‑1 (MQ=255)                 tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                                  <  1:324109/55‑1 (MQ=255)                  gAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATGTTTCGGCGTGATGAAACCaaa                                                           >  1:1341911/1‑61 (MQ=38)                  gAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                     >  1:1440872/1‑67 (MQ=255)                  gAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATATTTCGGCGTGATGAAACCAAAACCttt                                                     >  1:1281004/1‑67 (MQ=255)                   aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCGAAACCTTTTTc                                                  <  1:2127203/69‑1 (MQ=255)                   aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                  <  1:1026343/69‑1 (MQ=255)                   aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                 <  1:642689/70‑1 (MQ=255)                        cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                           <  1:547626/71‑1 (MQ=255)                          gagaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                         >  1:256963/1‑71 (MQ=255)                            gaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                        >  1:1501188/1‑70 (MQ=255)                             aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                               >  1:1246287/1‑62 (MQ=255)                               aTAGACAAACACATCTTTTCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                    <  1:2231721/71‑1 (MQ=255)                               aTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                               <  1:601740/60‑1 (MQ=255)                               aTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                    >  1:363315/1‑71 (MQ=255)                               aTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                    >  1:2115194/1‑71 (MQ=255)                                  gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCa                                 >  1:479987/1‑71 (MQ=255)                                   aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                   <  1:2791493/52‑1 (MQ=255)                                   aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                   <  1:2482070/52‑1 (MQ=255)                                     aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                              <  1:2277477/71‑1 (MQ=255)                                     aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                              <  1:2727818/71‑1 (MQ=255)                                       acacATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaaa                            <  1:2227212/71‑1 (MQ=255)                                        cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                           <  1:3266168/71‑1 (MQ=255)                                        cacaTCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACTAAAc                           >  1:1670177/1‑71 (MQ=255)                                        cacaTCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACTAAAc                           >  1:2214762/1‑71 (MQ=255)                                         acaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                           >  1:1952714/1‑70 (MQ=255)                                         acaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                           >  1:68072/1‑70 (MQ=255)                                         acaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                           >  1:801449/1‑70 (MQ=255)                                           aTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                               <  1:1034334/48‑1 (MQ=255)                                             cTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                              <  1:3571749/63‑1 (MQ=255)                                              tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                           >  1:1066563/1‑65 (MQ=255)                                               ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                   >  1:1842473/1‑40 (MQ=255)                                               ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                   >  1:2823302/1‑40 (MQ=255)                                               ttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCa                                 >  1:1605349/1‑58 (MQ=255)                                                tgctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGt                       <  1:2879782/67‑1 (MQ=255)                                                  ctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtttt                 <  1:841682/71‑1 (MQ=255)                                                     ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                        <  1:1463075/45‑1 (MQ=255)                                                     ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt               >  1:1378120/1‑70 (MQ=255)                                                     ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt               <  1:3508796/70‑1 (MQ=255)                                                     ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt               <  1:1004229/70‑1 (MQ=255)                                                     ccATCTTTCGGCGTCATGAAACCAAAACCTTTTTCAGGGTTAAAc                                        <  1:3410776/45‑1 (MQ=255)                                                       aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAg            <  1:2090371/71‑1 (MQ=255)                                                        tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                        >  1:2238685/1‑42 (MQ=255)                                                        tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGa           >  1:1199119/1‑71 (MQ=255)                                                         cTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                      <  1:3114156/43‑1 (MQ=255)                                                         cTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                      <  1:2309491/43‑1 (MQ=255)                                                         cTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGa           >  1:1120168/1‑70 (MQ=255)                                                             cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                                >  1:22512/1‑45 (MQ=255)                                                              ggCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                                  <  1:2020739/42‑1 (MQ=255)                                                              ggCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                                  <  1:1159330/42‑1 (MQ=255)                                                              ggCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                                  <  1:1717763/42‑1 (MQ=255)                                                               gagTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCATTTTATTAGACATAttt    <  1:3226215/69‑1 (MQ=255)                                                               gagTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCATTTTATTAGACATAttt    <  1:1639919/69‑1 (MQ=255)                                                               gCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt    <  1:2358719/71‑1 (MQ=255)                                                                cGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                                >  1:628724/1‑42 (MQ=255)                                                                cGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta   <  1:1436152/71‑1 (MQ=255)                                                                 gTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat  <  1:2450815/71‑1 (MQ=255)                                                                      gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtttt                 <  1:1065795/51‑1 (MQ=255)                                                                      gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtttt                 <  1:759609/51‑1 (MQ=255)                                                                      gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtttt                 <  1:3716857/51‑1 (MQ=255)                                                                      |                                                                  TTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTAT  >  Exported/980176‑980310

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