Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,758,108	T→C	23.0%	P76P (CCA→CCG)	def ←	peptide deformylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,758,108	0	T	C	23.0%	104.4 / 25.4	82	P76P (CCA→CCG)	def	peptide deformylase
Reads supporting (aligned to +/- strand):  ref base T (19/44);  new base C (8/11);  total (27/55)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CACGTTGTTCAGGGATCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTT  >  minE/2758041‑2758169                                                                    |                                                              cACGTTGTTCAGGGATCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTg                                                              <  1:759532/69‑1 (MQ=255)          cAGGGATCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTg                                                              <  1:154929/60‑1 (MQ=255)          cAGGGATCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTg                                                              <  1:787003/60‑1 (MQ=255)          cAGGGATCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTg                                                              <  1:1625676/60‑1 (MQ=255)          cAGGGATCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTg                                                              <  1:475647/60‑1 (MQ=255)          cAGGGATCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTg                                                              <  1:3368054/60‑1 (MQ=255)                tCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCCGga                                                            >  1:341250/1‑56 (MQ=255)                tCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCCGga                                                            >  1:2001396/1‑56 (MQ=255)                tCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCCGga                                                            >  1:2002746/1‑56 (MQ=255)                tCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCCGga                                                            >  1:1389616/1‑56 (MQ=255)                tCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCCGga                                                            >  1:3430314/1‑56 (MQ=255)                tCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCCGga                                                            >  1:1311607/1‑56 (MQ=255)                tCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCCGga                                                            >  1:778537/1‑56 (MQ=255)                tCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCCGga                                                            >  1:797926/1‑56 (MQ=255)                        cAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTaa                                                    <  1:1385408/56‑1 (MQ=255)                        cAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTaa                                                    <  1:226478/56‑1 (MQ=255)                        cAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTaa                                                    <  1:1127111/56‑1 (MQ=255)                        cAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTc                                        >  1:3213786/1‑68 (MQ=255)                                      gCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGt                                      >  1:989928/1‑56 (MQ=255)                                      gCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCg                          >  1:2158494/1‑68 (MQ=255)                                      gCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         <  1:1456760/69‑1 (MQ=255)                                      gCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:1411746/1‑69 (MQ=255)                                      gCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         <  1:2344123/69‑1 (MQ=255)                                      gCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:578721/1‑69 (MQ=255)                                      gCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         <  1:1667389/69‑1 (MQ=255)                                      gCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         <  1:284736/69‑1 (MQ=255)                                      gCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:3125501/1‑69 (MQ=255)                                       gcTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:1679949/2‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:3424160/1‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:202930/1‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:2413017/1‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:1825955/1‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:1782225/1‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:3021341/1‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:1302319/1‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:1614874/1‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:980557/1‑68 (MQ=255)                                       ccTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGa                         >  1:1501820/1‑68 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACa                     <  1:1527709/66‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACa                     <  1:489026/66‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACa                     <  1:2514574/66‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACa                     <  1:128217/66‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACa                     <  1:2359567/66‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACa                     <  1:1098847/66‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACa                     <  1:1668105/66‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATc                   >  1:1918171/1‑68 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATc                   <  1:1883663/68‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATc                   <  1:144059/68‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATc                   <  1:652111/68‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATc                   <  1:1708214/68‑1 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATc                   >  1:340221/1‑68 (MQ=255)                                             tCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATc                   <  1:1352751/68‑1 (MQ=255)                                                     tttttCTAAAAGCTCTTGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:1827712/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:1437778/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:1003216/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:907126/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:1075604/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:1636865/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:528528/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:388000/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:3390668/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:1429452/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:332521/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:1848124/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:2876720/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATAAATAACAAt          <  1:899171/69‑1 (MQ=255)                                                     tttttCTAAAAGCTCTGGATTGATTAACACTAGACGTTCGTCACGGTTTTCCGAAACATCAATAACAAt          <  1:1351348/69‑1 (MQ=255)                                                        ttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGa        <  1:3416729/68‑1 (MQ=255)                                                        ttCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGa        <  1:1232390/68‑1 (MQ=255)                                                               aGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGtt  <  1:3355842/67‑1 (MQ=255)                                                               aGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGtt  <  1:1363977/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:1748522/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:1719936/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:68191/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:2637298/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:418886/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:1849190/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:1324229/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:2306324/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:1533249/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGTCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:1719093/67‑1 (MQ=255)                                                               aGCTCCGGATTGATTAACACCAGCCGTTCGGCACGGTTTTCCGAAACATCAATTACAATGATACGtt  <  1:1667239/67‑1 (MQ=255)                                                                    |                                                              CACGTTGTTCAGGGATCGACAGGCAACCTTCTTCAATGCCTGTTTCGCCGCTTTTTTCTAAAAGCTCTGGATTGATTAACACTAGCCGTTCGTCACGGTTTTCCGAAACATCAATAACAATGATACGTT  >  minE/2758041‑2758169

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