Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,484,894	T→C	81.6%	M316V (ATG→GTG)	menF ←	isochorismate synthase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,484,894	0	T	C	81.6%	137.4 / 22.3	65	M316V (ATG→GTG)	menF	isochorismate synthase 2
Reads supporting (aligned to +/- strand):  ref base T (9/3);  new base C (27/26);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAG  >  Exported/1484826‑1484964                                                                     |                                                                       cTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCActt                                                                      <  1:2646895/71‑1 (MQ=255) cTGATGTAAACAGATCACATCATCCGCTTTGTTGAGGGAAGTCCAGATACAGCGGCGAAGATGCTGCAttt                                                                      <  1:1241072/71‑1 (MQ=255)  tGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCActt                                                                      >  1:710570/1‑70 (MQ=255)   gATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTa                                                                    >  1:1561087/1‑71 (MQ=255)    aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCAtttt                                                                     <  1:192956/69‑1 (MQ=255)    aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTac                                                                   >  1:5081432/1‑71 (MQ=255)    aTGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTac                                                                   >  1:3223716/1‑71 (MQ=255)    aTGTAAACAGATCACATCATCCGCTTTGTT‑AGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacg                                                                  >  1:5399660/1‑71 (MQ=255)     tGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacg                                                                  <  1:2054417/71‑1 (MQ=255)      gTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacg                                                                  >  1:2817261/1‑70 (MQ=255)         aaCAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcaga                                                              >  1:2267276/1‑71 (MQ=255)           cAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTacgcagacg                                                            <  1:1598365/71‑1 (MQ=255)           cAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacg                                                            <  1:3474297/71‑1 (MQ=255)           cAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacg                                                            <  1:4211750/71‑1 (MQ=255)              aTCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacgcag                                                         >  1:5185619/1‑71 (MQ=255)              aTCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacgcag                                                         >  1:5566577/1‑71 (MQ=255)               tCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacg                                                                  <  1:156262/61‑1 (MQ=255)                cacaTCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTa                                                       >  1:5080572/1‑71 (MQ=255)                cacaTCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTa                                                       >  1:1321123/1‑71 (MQ=255)                 acaTCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacgc                                                           <  1:4359252/66‑1 (MQ=255)                 acaTCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTa                                                       <  1:2628214/70‑1 (MQ=255)                 acaTCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTAc                                                      >  1:5387568/1‑71 (MQ=255)                  catcatCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTAc                                                      >  1:3711945/1‑70 (MQ=255)                  catcatCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTAc                                                      >  1:2807809/1‑70 (MQ=255)                   atcatcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacgcag                                                         >  1:391355/1‑66 (MQ=255)                   atcatcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcagacgcag                                                         >  1:5450171/1‑66 (MQ=255)                    tcatcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcaga                                                              <  1:124239/60‑1 (MQ=255)                      atcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTgcg                                                 <  1:2549878/71‑1 (MQ=255)                      atcCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTgc                                                  <  1:4701096/70‑1 (MQ=255)                             ttgttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTa                                           >  1:3497263/1‑70 (MQ=255)                              tgttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTaa                                          <  1:2309862/70‑1 (MQ=255)                              tgttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTaa                                          <  1:1958275/70‑1 (MQ=255)                              tgttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTaa                                          <  1:4781499/70‑1 (MQ=255)                               gttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTacgcaga                                                              >  1:5676481/1‑49 (MQ=255)                                ttgAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTaaaa                                        <  1:2926830/70‑1 (MQ=255)                                  gAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTa                                                                    >  1:2324694/1‑40 (MQ=255)                                  gAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACa                                      <  1:2252377/70‑1 (MQ=255)                                   aGTGAAGTCCAGATACAGCGGCGAAGATGATGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATc                                    <  1:3580520/71‑1 (MQ=255)                                    gTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATcc                                   <  1:1939265/71‑1 (MQ=255)                                    gTGAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATcc                                   <  1:2346247/71‑1 (MQ=255)                                      gAAGTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGt                                                        >  1:3319777/1‑48 (MQ=255)                                       aaGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                   >  1:243469/1‑36 (MQ=255)                                       aaGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                   >  1:139551/1‑36 (MQ=255)                                       aaGTCCAGATACAGCGGCGAAGATGCTGCATTTTac                                                                   >  1:1968206/1‑36 (MQ=255)                                       aaGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGc                                >  1:3474869/1‑71 (MQ=255)                                         gTCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCt                                                    <  1:2628766/49‑1 (MQ=255)                                          tCCAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTaaa                                         <  1:1836438/59‑1 (MQ=255)                                           ccAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTaaa                                         <  1:1595759/58‑1 (MQ=255)                                            cAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCt                            <  1:2344709/70‑1 (MQ=255)                                            cAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTg                           >  1:3892664/1‑71 (MQ=255)                                            cAGATACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTg                           >  1:5377672/1‑71 (MQ=255)                                               aTACAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGAGGTAAAACATCCAGCGTCTGGGt                        <  1:4797586/71‑1 (MQ=255)                                                tACAGCGGCGAAGATGCTGCATTTTacgcagacgcag                                                         >  1:1810862/1‑37 (MQ=255)                                                 aCAGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTAt                      >  1:5222770/1‑71 (MQ=255)                                                 aCAGCGGCGAAGATACTGCACTTTacgcagacgca                                                          >  1:3383896/1‑35 (MQ=38)                                                   aGCGGCGAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCg                    <  1:4345916/71‑1 (MQ=255)                                                    gcggcgAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCgg                   >  1:901315/1‑71 (MQ=255)                                                    gcggcgAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGt                        >  1:1445508/1‑66 (MQ=255)                                                    gcggcgAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGt                        >  1:677519/1‑66 (MQ=255)                                                     cggcgAAGATGCTGCACTTTACGCAGACGCAGTACCTgc                                                  <  1:3944595/39‑1 (MQ=255)                                                     cggcgAAGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGc                  <  1:1764479/71‑1 (MQ=255)                                                          aaGATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTa             <  1:1902223/71‑1 (MQ=255)                                                           aGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAg                                 >  1:3961135/1‑50 (MQ=255)                                                            gATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTaaaa                                        <  1:5290930/42‑1 (MQ=255)                                                            gATGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCt                >  1:673664/1‑66 (MQ=255)                                                             aTGCTGCACTTTACGCAGACGCAGTACCTGCGGCGGTaa                                          <  1:2159032/39‑1 (MQ=255)                                                              tgctgcACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTa                       >  1:5471341/1‑57 (MQ=255)                                                               gctgcACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACa                                      >  1:1758372/1‑41 (MQ=255)                                                                  gcACTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGa     >  1:419023/1‑71 (MQ=255)                                                                     ttttACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAg  <  1:2848449/71‑1 (MQ=255)                                                                     |                                                                       CTGATGTAAACAGATCACATCATCCGCTTTGTTGAGTGAAGTCCAGATACAGCGGCGAAGATGCTGCATTTTACGCAGACGCAGTACCTGCGGCGGTAAAACATCCAGCGTCTGGGTATCGGCCTGTAATCGTTGACAG  >  Exported/1484826‑1484964

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: ‑