Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	805,466	T→C	9.1%	N901D (AAC→GAC)	ycgV ←	predicted adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	805,466	0	T	C	9.1%	150.0 / 5.8	55	N901D (AAC→GAC)	ycgV	predicted adhesin
Reads supporting (aligned to +/- strand):  ref base T (33/17);  new base C (2/3);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGC  >  minE/805398‑805534                                                                     |                                                                     cTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTAGCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCa                                                                    <  1:2086231/71‑1 (MQ=255)         tGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAAc                                                                  >  1:3574225/1‑65 (MQ=255)           aCTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTc                                                          >  1:1216701/1‑71 (MQ=255)           aCTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTc                                                          >  1:2998589/1‑71 (MQ=255)           aCTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTCCAACAGATAtt                                                           <  1:3040438/70‑1 (MQ=255)             tCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCgg                                                        >  1:3704716/1‑71 (MQ=255)               ccACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCgg                                                        <  1:3816014/69‑1 (MQ=255)               ccACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTa                                                      >  1:2402713/1‑71 (MQ=255)               ccACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGCTGTTCAACAGATATTCGGTa                                                      >  1:1051396/1‑71 (MQ=255)                cACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTAt                                                     >  1:2576850/1‑71 (MQ=255)                 aCGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATc                                                    <  1:1740483/71‑1 (MQ=255)                  cGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATcc                                                   >  1:3886591/1‑71 (MQ=255)                  cGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATcc                                                   >  1:3540855/1‑71 (MQ=255)                  cGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATcc                                                   >  1:2139618/1‑71 (MQ=255)                    ccGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATcccc                                                 >  1:962568/1‑71 (MQ=255)                     cGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCt                                                >  1:2742163/1‑71 (MQ=255)                         aTTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGaa                                             >  1:3014946/1‑70 (MQ=255)                          ttCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGaaaa                                           >  1:2007482/1‑71 (MQ=255)                          ttCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGaaaa                                           <  1:3814278/71‑1 (MQ=255)                          ttCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGaaaa                                           <  1:1227889/71‑1 (MQ=255)                          ttCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGaaa                                            <  1:1169127/70‑1 (MQ=255)                          ttCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGa                                              <  1:594201/68‑1 (MQ=255)                          ttCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTCCAACAGATATTCGGTATCCCCTGaaaa                                           >  1:2728345/1‑71 (MQ=255)                            ccAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATccc                                                  <  1:275832/62‑1 (MQ=255)                             cAGCCATTACCTTTGAAACTGTACTTCTCCCTGGAGTTGTTCAACAGATATTCGGTATcccc                                                 <  1:2318297/62‑1 (MQ=255)                             cAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCaa                                                                   >  1:1636261/1‑44 (MQ=255)                              aGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTa                                                      >  1:2627743/1‑56 (MQ=255)                               gCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTg                                               >  1:1009674/1‑62 (MQ=255)                               gCCATTACCTTTGAAACAGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTcgc                                      >  1:1824644/1‑71 (MQ=255)                                     aCCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATc                                                    >  1:853860/1‑51 (MQ=255)                                       cTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGata                                                             <  1:1796442/40‑1 (MQ=255)                                          tGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGaa                                             >  1:991915/1‑53 (MQ=255)                                          tGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTc                                        <  1:1723464/58‑1 (MQ=255)                                           gAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTcg                                       >  1:75314/1‑58 (MQ=255)                                           gAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCa                           <  1:2570219/70‑1 (MQ=255)                                            aaaCTGTACTTCTCCCGGGAGTTGTCCAACAGATAATCGGTATCCCCTGAAAACTCGCGGATAGCGCCAg                          <  1:1127972/70‑1 (MQ=255)                                              aCTGCACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCt                       <  1:3596245/71‑1 (MQ=255)                                               cTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTcgcg                                     >  1:1940557/1‑56 (MQ=255)                                               cTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCtt                      <  1:3401744/71‑1 (MQ=255)                                               cTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCtt                      <  1:2723822/71‑1 (MQ=255)                                                tGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTc                                        >  1:665151/1‑52 (MQ=255)                                                tGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTc                                        >  1:660438/1‑52 (MQ=255)                                                tGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGAATCCCCTGAAAACt                                         >  1:1110779/1‑51 (MQ=255)                                                  tACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCt                                                >  1:549796/1‑42 (MQ=255)                                                     ttCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTc                                        <  1:245706/47‑1 (MQ=255)                                                      tctcCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGaaaa                                           >  1:2828880/1‑43 (MQ=255)                                                      tctcCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCgg                                    <  1:786667/50‑1 (MQ=255)                                                      tctcCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTcaca                  >  1:3030491/1‑68 (MQ=255)                                                          ccGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACAt           >  1:1562534/1‑71 (MQ=255)                                                           cGGGAGTTGTCCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAg               <  1:1634499/66‑1 (MQ=255)                                                            gggAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACAtt          >  1:1719027/1‑70 (MQ=255)                                                            gggAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACAtt          >  1:3015450/1‑70 (MQ=255)                                                            gggAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACAtt          >  1:3527777/1‑70 (MQ=255)                                                            gggAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACAtt          >  1:335727/1‑70 (MQ=255)                                                               aGTTGTCCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgc      >  1:1422837/1‑71 (MQ=255)                                                                gttgtcCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgct     <  1:2770073/71‑1 (MQ=255)                                                                 ttgttCAACAGATATTCGGTATCCCCTGAAAACTcg                                       <  1:511390/36‑1 (MQ=255)                                                                 ttgttCAACAGATATTCGGTATCCCCTGAAAACTcg                                       <  1:2886056/36‑1 (MQ=255)                                                                   gttCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACAt           >  1:1441277/1‑62 (MQ=255)                                                                   gttCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAgctggc  >  1:170204/1‑71 (MQ=255)                                                                     |                                                                     CTGTGCACTGACTCCCACGCCGTTATTCCAGCCATTACCTTTGAAACTGTACTTCTCCCGGGAGTTGTTCAACAGATATTCGGTATCCCCTGAAAACTCGCGGATAGCGCCAGTCTTCACATAGACATTCAGCTGGC  >  minE/805398‑805534

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