Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,072,378	G→C	E39Q (GAA→CAA)	unknown →	Fumarate and nitrate reduction regulatory protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,072,378	0	G	C	96.8%	80.7 / ‑3.7	31	E39Q (GAA→CAA)	unknown	Fumarate and nitrate reduction regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (17/13);  minor base A (0/1);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GTTTTTCAGCAGCCATATATCCTCACTTATGTTTATACTTCTTAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATGAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAGTACTTTTTAAAGCGGGCGATGAGCTGACCTCACTCTACGCTATCCGTTCCGGCACAATCAAATCTTATACCATTAGCGAAACAGGCGAAGAACAAATTACCTCTTTTCATTTACCGGGCGATTTGGTGGGGTTCGATGCCATTATGAATATGCAACATCCGAGTT  >  1/1072140‑1072606                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     gTTTTTCAGCAGCCATATATCCTCACTTATGTTTATACTTCTTAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAg                                                                                                                                                                                                                          >  1:253194/1‑251 (MQ=255)  tttttCAGCAGCCATATATCCTCACTTATGTTTATACTTCTTAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAAtt                                                                                                                                                                                                                            <  1:860914/248‑1 (MQ=255)            gcCATATATCCTCACTTATGTTTATACTTCTTAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGAt                                                                                                                                                                                                                        >  1:130730/1‑242 (MQ=255)                        cACTTATGTTTATACTTCTTAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGataat                                                                                                                                                                                                            <  1:511618/242‑1 (MQ=255)                             aTGTTTATACTTCTTAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATc                                                                                                                                                                                                                       >  1:165178/1‑226 (MQ=255)                              tGTTTATACTTCTTAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTaaa                                                                                                                                                                                             >  1:18233/1‑251 (MQ=255)                                        tcttAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATc                                                                                                                                                                                                                       >  1:356633/1‑215 (MQ=255)                                        tcttAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATAc                                                                                                                                                                                   <  1:146046/251‑1 (MQ=255)                                                    gAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAg                                                                                                                                                                       >  1:251297/1‑251 (MQ=255)                                                        tAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATc                                                                                                                                                                                                                       <  1:616773/199‑1 (MQ=255)                                                        tAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTa                                                                                                                                                                                                                  <  1:496473/204‑1 (MQ=255)                                                         aGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATATATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTg                                                                                                                                                                                                     >  1:199040/1‑216 (MQ=255)                                                          gTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTg                                                                                                                                                                                                     >  1:719143/1‑215 (MQ=255)                                                          gTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAGTACttt                                                                                                                                                                 >  1:374334/1‑251 (MQ=255)                                                             gTACTATTGAGCTAGTTCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGTGGCATAAATGAAAATTTTAGCGACAGAAATAAAATTAGGACGTCGCGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTTTAGTATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAATACGAATTTGATCAGTTAGATAATATTATTGAACGTCAAAAACCTATATAGTAATCCCAAGTACTTTTTa                                                                                                                                                              >  1:293302/1‑251 (MQ=255)                                                             gTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATc                                                                                                                                                                                                                       >  1:688226/1‑194 (MQ=255)                                                             gTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGt                                                                                                                                                                                                >  1:182106/1‑217 (MQ=255)                                                                cTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAAtattat                                                                                                                                                                                                       <  1:277919/207‑1 (MQ=255)                                                                cTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATAAACACGAATTAGa                                                                                                                                                                                                                         <  1:879155/189‑1 (MQ=255)                                                                           gATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAg                                                                                                                                                                                 >  1:901159/1‑218 (MQ=255)                                                                                                           cTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGataa                                                                                                                                                                                                             <  1:666493/158‑1 (MQ=255)                                                                                                           cTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAatat                                                                                                                                                                                                          <  1:772809/161‑1 (MQ=255)                                                                                                                                   aTTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAGTACTTTTTAAAGCGGGCGATGAGCTGACCTCACTCTACg                                                                                                                                <  1:20651/211‑1 (MQ=255)                                                                                                                                                                    gTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAGTACTTTTTAAAGCGGGCGATGAGCTGACCTCACTCTACGCTATCCGTTCCGGCACAATCAAATCTTATACCATTAGCGAAACAGGCGAAGAACAAATTACCTCTTTTCAttt                                                       >  1:52833/1‑251 (MQ=255)                                                                                                                                                                      cggtaccgggcatttgcgtaataattGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAGTACTTTTTAAAGCGGGCGATGAGCTGACCTCACTCTACGCTATCCGTTCCGGCACAATCAAATCTTATACCATTAGCGAAACAGGCGAAGAACAAATTACCTCTTTTCATTTAc                                                     <  1:217632/228‑1 (MQ=255)                                                                                                                                                                         aTCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAAcc                                                                                                                                                                                        <  1:854007/117‑1 (MQ=255)                                                                                                                                                                          tCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGataa                                                                                                                                                                                                             >  1:599661/1‑95 (MQ=255)                                                                                                                                                                                                 ccAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAGTACTTTTTAAAGCGGGCGATGAGCTGACCTCACTCTACGCTATCCGTTCCGGCACAATCAAATCTTATACCATTAGCGAAACAGGCGAAGAACAAATTACCTCTTTTCATTTACCGGGCGATTTGGTGGGGTTCGATGCCa                          <  1:32738/251‑1 (MQ=255)                                                                                                                                                                                                            gCATTAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAg                                                                                                                                                                                 >  1:791157/1‑89 (MQ=255)                                                                                                                                                                                                                tAGCCAGCTATGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAGTACTTTTTAAAGCGGGCGATGAGCTGACCTCACTCTACGCTATCCGTTCCGGCACAATCAAATCTTATACCATTAGCGAAACAGGCGAAGAACAAATTACCTCTTTTCATTTACCGGGCGATTTGGTGGGGTTCGATGCCATTATGAATATGCAAc           <  1:113197/251‑1 (MQ=255)                                                                                                                                                                                                                         aTGTATTCCTTTTACTTTGAATCAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAGTACTTTTTAAAGCGGGCGATGAGCTGACCTCACTCTACGCTATCCGTTCCGGCACAATCAAATCTTATACCATTAGCGAAACAGGCGAAGAACAAATTACCTCTTTTCATTTACCGGGCGATTTGGTGGGGTTCGATGCCATTATGAATATGCAACATCCGAGtt  >  1:4953/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                     GTTTTTCAGCAGCCATATATCCTCACTTATGTTTATACTTCTTAAATTAATGAATTAGTGGTACTATTGAGCTAGATCAATTTTACACGATTTTTTTTCAATTTTTCTGTAATTTGGGGCATAAATGAAAATTTTAGCGACAGAAACAAAATTAGGACGTCGAGTTCAATCCGGGGGGTGTGCCATTCATTGCCAAAATTGTAGCATTAGCCAGCTATGTATTCCTTTTACTTTGAATGAACACGAATTAGATCAGTTAGATAATATTATTGAACGTAAAAAACCTATACAGAAATCCCAAGTACTTTTTAAAGCGGGCGATGAGCTGACCTCACTCTACGCTATCCGTTCCGGCACAATCAAATCTTATACCATTAGCGAAACAGGCGAAGAACAAATTACCTCTTTTCATTTACCGGGCGATTTGGTGGGGTTCGATGCCATTATGAATATGCAACATCCGAGTT  >  1/1072140‑1072606

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑