Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,695,039	C→T	16.4%	F31F (TTC→TTT)	nudE →	ADP‑ribose diphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,695,039	0	C	T	16.4%	144.4 / 18.8	61	F31F (TTC→TTT)	nudE	ADP‑ribose diphosphatase
Reads supporting (aligned to +/- strand):  ref base C (22/29);  new base T (7/3);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

CCACCATTCTGAATGTTGAAACTGTAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGTGATGATTGTGC  >  minE/2694969‑2695106                                                                       |                                                                    ccaccaTTCTGAATGTTGAAACTGTAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTc                                                                     <  1:486025/71‑1 (MQ=255)   accaTTCTGAATGTTGAAACTGTAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAg                                                                   <  1:2564587/71‑1 (MQ=255)   accaTTCTGAATGTTGAAACTGTAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAg                                                                   <  1:3808633/71‑1 (MQ=255)     caTTCTGAATGTTGAAACTGTAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTTAGc                                                                  >  1:2717012/1‑70 (MQ=255)       ttCTGAATGTTGAAACTGTAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCaa                                                                <  1:2413815/70‑1 (MQ=255)            aaTGTTGAAACTGTAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGt                                                          >  1:1304983/1‑71 (MQ=255)                  gAAACTGTAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCg                                                           <  1:3882433/64‑1 (MQ=255)                         tAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGttt                                               <  1:1213489/69‑1 (MQ=255)                         tAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTa                                              <  1:2919667/70‑1 (MQ=255)                         tAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTAt                                             <  1:1337178/71‑1 (MQ=255)                           gCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTTAGCAATGGCGTGCGGCGTGTTTATg                                            <  1:1304003/70‑1 (MQ=255)                            cccGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGaa                                          >  1:3008310/1‑71 (MQ=255)                              cGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAAc                                         >  1:3734453/1‑70 (MQ=255)                               gTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGa                                       <  1:1640498/71‑1 (MQ=255)                               gTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGa                                       <  1:2647934/71‑1 (MQ=255)                                ttCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGc                                                                  <  1:631672/43‑1 (MQ=255)                                ttCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATgg                                                             <  1:3078953/48‑1 (MQ=255)                                ttCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTgcggc                                                     <  1:3383539/56‑1 (MQ=255)                                 tCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGa                                       <  1:2253713/69‑1 (MQ=255)                                  cccGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTTAGCAATgg                                                             >  1:1103337/1‑46 (MQ=255)                                  cccGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCaa                                                                >  1:143888/1‑43 (MQ=255)                                   ccGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGa                                       >  1:3568539/1‑67 (MQ=255)                                     gACTGTTTACCGTCGAGAGCTTGGATCTGGAGTTCAGCAATGGc                                                            <  1:3192365/44‑1 (MQ=255)                                     gACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGc                                                            <  1:2088042/44‑1 (MQ=255)                                     gACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGc                                                            <  1:2149486/44‑1 (MQ=255)                                     gACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGt                                 <  1:1927230/71‑1 (MQ=255)                                     gACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGt                                 <  1:2183043/71‑1 (MQ=255)                                      aCTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGc                                                                  <  1:1928385/37‑1 (MQ=255)                                       cTGTTTACCGTCGAGAGCGTGGATCTGGAGTTTAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTcc                               <  1:3361534/71‑1 (MQ=255)                                       cTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTcc                               <  1:1116367/71‑1 (MQ=255)                                       cTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTcc                               <  1:1117170/71‑1 (MQ=255)                                       cTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTc                                >  1:2423721/1‑70 (MQ=255)                                        tGTTTACCGTCGAGAGCGTGGATCTGGAGTTTAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTcc                               >  1:1358887/1‑70 (MQ=255)                                          tttACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTccaac                            >  1:3356969/1‑71 (MQ=255)                                          tttACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTccaac                            >  1:3035107/1‑71 (MQ=255)                                          tttACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTcca                              <  1:3597188/69‑1 (MQ=255)                                           ttACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTccca                             >  1:251955/1‑67 (MQ=255)                                            tACCGTCGAGAGCGTGGATCTGGAGTTTAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTccaacc                           >  1:3593708/1‑70 (MQ=255)                                            tACAGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGtt                                                >  1:2424792/1‑49 (MQ=255)                                             aCCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTccaacca                          <  1:3250929/70‑1 (MQ=255)                                              ccGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCcaaccaac                        >  1:986626/1‑71 (MQ=255)                                                gTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCg                      >  1:945349/1‑71 (MQ=255)                                                 tCGAGAGCTTGGATCTGGAGTTTAGCAATGGCGTGCGGCGTGTTTATGAACGa                                       <  1:1966791/53‑1 (MQ=255)                                                 tCGAGAGCGTGGATCTGGAGTTTAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTccaa                             >  1:2822862/1‑63 (MQ=255)                                                 tCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCgg                     >  1:378981/1‑71 (MQ=255)                                                   gagaGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCggg                    <  1:815317/70‑1 (MQ=255)                                                         gTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGt             <  1:2934964/71‑1 (MQ=255)                                                          tGGATCTGGAGTTTAGCAATGGCGTGCGGCGTGTTTATGAACGaa                                      >  1:1285016/1‑45 (MQ=255)                                                           ggATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGa                                           <  1:3265835/39‑1 (MQ=255)                                                           ggATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGtga           >  1:3317439/1‑71 (MQ=255)                                                           ggATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGt             >  1:1973401/1‑69 (MQ=255)                                                            gATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGtgat          >  1:1619248/1‑71 (MQ=255)                                                             aTCTGGAGTTTAGCAATGGCGTGCGGTGTGTTTATGAACGAATGCGTccaa                             >  1:1166000/1‑51 (MQ=255)                                                                 ggAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACg                                        <  1:300321/36‑1 (MQ=255)                                                                 ggAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCggg                    <  1:1809149/56‑1 (MQ=255)                                                                 ggAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGTGATGATtg     >  1:1545151/1‑71 (MQ=255)                                                                  gAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGtg            >  1:2444347/1‑63 (MQ=255)                                                                  gAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGTGATGATtgt    <  1:2007907/71‑1 (MQ=255)                                                                  gAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGTGATGATtgt    >  1:1000786/1‑71 (MQ=255)                                                                   aGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGTGATGATtgtg   >  1:3517702/1‑71 (MQ=255)                                                                    gTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGTGATGATTGTGc  >  1:637304/1‑71 (MQ=255)                                                                     ttCAGCAATGGCGTGCGGCGTGTTTATGATCGAATGCGTccacc                            >  1:444630/1‑42 (MQ=39)                                                                       cAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCaaccaacc                       <  1:3660663/47‑1 (MQ=255)                                                                       |                                                                    CCACCATTCTGAATGTTGAAACTGTAGCCCGTTCCCGACTGTTTACCGTCGAGAGCGTGGATCTGGAGTTCAGCAATGGCGTGCGGCGTGTTTATGAACGAATGCGTCCAACCAACCGGGAAGCAGTGATGATTGTGC  >  minE/2694969‑2695106

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