Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,772,254	C→T	15.5%	P527S (CCA→TCA)	aceB →	malate synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,772,254	0	C	T	15.5%	139.1 / 16.1	58	P527S (CCA→TCA)	aceB	malate synthase A
Reads supporting (aligned to +/- strand):  ref base C (27/22);  new base T (4/5);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TTTTGACGATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACACAA  >  minE/2772184‑2772324                                                                       |                                                                       ttttGACGATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGc                                                                        <  1:2971736/71‑1 (MQ=255)   ttGACGATGCCGCACGCTTGATGGAACAGATCCCCACTTCCGATGAGTTAATTGATTTCCTGACCCTGc                                                                        >  1:2489890/1‑69 (MQ=255)   ttGACGATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGcc                                                                       <  1:2514777/70‑1 (MQ=255)   ttGACGATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGcc                                                                       <  1:155478/70‑1 (MQ=255)   ttGACGATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGcc                                                                       <  1:3495110/70‑1 (MQ=255)    tGACGATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAg                                                                     >  1:3138654/1‑71 (MQ=255)     gACGATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAg                                                                     >  1:2070513/1‑70 (MQ=255)        gATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGc                                                                        <  1:1642565/64‑1 (MQ=255)        gATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTa                                                                 <  1:3257187/71‑1 (MQ=255)          tGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTAcc                                                               >  1:1296895/1‑71 (MQ=255)          tGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTAcc                                                               >  1:223612/1‑71 (MQ=255)          tGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTAcc                                                               >  1:1976570/1‑71 (MQ=255)               cACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTAc                                                                >  1:1098469/1‑65 (MQ=255)               cACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTAc                                                                >  1:3144562/1‑65 (MQ=255)               cACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTg                                                          >  1:2167480/1‑71 (MQ=255)                 cGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGt                                                         >  1:493617/1‑70 (MQ=255)                  gCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGt                                                         <  1:4412/69‑1 (MQ=255)                  gCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTa                                                       <  1:167551/71‑1 (MQ=255)                  gCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTa                                                       <  1:2838010/71‑1 (MQ=255)                    ttGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGTCAGGCTACCGCCTGTTAGc                                                     >  1:2980448/1‑71 (MQ=255)                    ttGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGc                                                     >  1:2141006/1‑71 (MQ=255)                      gATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCg                                                    >  1:1569852/1‑70 (MQ=255)                       aTGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTa                                                  <  1:3774298/71‑1 (MQ=255)                        tGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTaa                                                 <  1:954901/71‑1 (MQ=255)                        tGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTa                                                  >  1:2606510/1‑70 (MQ=255)                          gAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGTCAGGCTACCGCCTGTTAGCGTaaa                                                <  1:304819/70‑1 (MQ=255)                          gAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGTCAGGCTACCGCCTGTTAGCGTaaa                                                <  1:349542/70‑1 (MQ=255)                          gAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGTCAGGCTACCGCCTGTTAGCGTaaa                                                <  1:3727103/70‑1 (MQ=255)                             cAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAaccac                                            <  1:3149993/71‑1 (MQ=255)                             cAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAaccac                                            <  1:1799047/71‑1 (MQ=255)                                aTCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGc                                                                        <  1:2471277/40‑1 (MQ=255)                                aTCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGc                                                                        <  1:3271661/40‑1 (MQ=255)                                aTCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAAccacca                                          >  1:1382373/1‑70 (MQ=255)                                 tCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGcc                                                                       <  1:476548/40‑1 (MQ=255)                                 tCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGcc                                                                       <  1:610201/40‑1 (MQ=255)                                 tCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGcc                                                                       <  1:987326/40‑1 (MQ=255)                                   accacTTCCGATGAGTTAATTGATTTCCTGACCCTGTCAGGCTACCGCCTGTTAGCGTAAACCACCACAt                                       >  1:1187931/1‑70 (MQ=255)                                      acTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTaa                                                 <  1:443497/57‑1 (MQ=255)                                      acTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACAt                                       >  1:2720335/1‑67 (MQ=255)                                       cTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGc                                                     <  1:2651644/52‑1 (MQ=255)                                       cTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTa                                  >  1:1796408/1‑71 (MQ=255)                                         tCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGt                                                   >  1:1893376/1‑52 (MQ=255)                                         tCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATg                                >  1:2733981/1‑71 (MQ=255)                                          ccGATGAGTTAATTGATTTCCTGACCCTGTCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATg                                <  1:2931453/70‑1 (MQ=255)                                          ccGATGAGTTAATTGATTTCCTGACCCTGTCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATg                                <  1:3630956/70‑1 (MQ=255)                                          ccGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGt                                                         <  1:1559635/45‑1 (MQ=255)                                          ccGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATg                                <  1:2202147/70‑1 (MQ=255)                                          ccGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCATAAACCACCACATAACTATg                                >  1:920411/1‑70 (MQ=255)                                             aTGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAg                             >  1:2360167/1‑70 (MQ=255)                                              tGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCa                           <  1:2631652/71‑1 (MQ=255)                                                aGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCAt                          <  1:3452307/70‑1 (MQ=255)                                                 gTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACAt                                       <  1:3931070/56‑1 (MQ=255)                                                    aaTTGATTTCCTGACCCTGTCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGca                     >  1:1990920/1‑71 (MQ=255)                                                    aaTTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGa                              >  1:1685906/1‑62 (MQ=255)                                                    aaTTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGc                      <  1:3439954/70‑1 (MQ=255)                                                    aaTTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGc                      >  1:3912083/1‑70 (MQ=255)                                                         aTTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATg                 <  1:2241421/70‑1 (MQ=255)                                                         aTTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATg                 >  1:916887/1‑70 (MQ=255)                                                          tttCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTAt                                 >  1:1644546/1‑53 (MQ=255)                                                          tttCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGaa               <  1:3329012/71‑1 (MQ=255)                                                          tttCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGaa               <  1:2397038/71‑1 (MQ=255)                                                           ttCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGaa               >  1:2156141/1‑70 (MQ=255)                                                           ttCCTGACCCTGCCAGGCTAACGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGaa               >  1:1378647/1‑70 (MQ=255)                                                            tCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGaaa              <  1:988585/70‑1 (MQ=255)                                                            tCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGaaa              <  1:214592/70‑1 (MQ=255)                                                              cTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAAcc           >  1:1102560/1‑71 (MQ=255)                                                                  cccTGTCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGaa               >  1:2576272/1‑63 (MQ=255)                                                                  cccTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTa       <  1:3824359/71‑1 (MQ=255)                                                                     tGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTaca     >  1:3447614/1‑70 (MQ=255)                                                                       ccAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGaa               <  1:2387876/58‑1 (MQ=255)                                                                       ccAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTacac    <  1:1715038/69‑1 (MQ=255)                                                                       ccAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTAcaca   >  1:2514473/1‑70 (MQ=255)                                                                       ccAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACacaa  >  1:2478457/1‑71 (MQ=255)                                                                       ccAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACacaa  >  1:1251528/1‑71 (MQ=255)                                                                       |                                                                       TTTTGACGATGCCGCACGCTTGATGGAACAGATCACCACTTCCGATGAGTTAATTGATTTCCTGACCCTGCCAGGCTACCGCCTGTTAGCGTAAACCACCACATAACTATGGAGCATCTGCACATGAAAACCCGTACACAA  >  minE/2772184‑2772324

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