Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	980,236	C→A	10.6%	P23P (CCG→CCT)	cspI ←	cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	980,236	0	C	A	10.6%	236.9 / 9.7	75	P23P (CCG→CCT)	cspI	cold shock protein
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base A (2/6);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAA  >  minE/980169‑980306                                                                    |                                                                       ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                     >  1:2636296/1‑71 (MQ=255) ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                     <  1:1732959/71‑1 (MQ=255) ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                     <  1:3820552/71‑1 (MQ=255)    aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                                                  <  1:3160479/71‑1 (MQ=255)    aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                                                  <  1:1122383/71‑1 (MQ=255)         ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                              <  1:2133592/70‑1 (MQ=255)         ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                              <  1:621809/70‑1 (MQ=255)         ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                                              <  1:266283/70‑1 (MQ=255)            aaaTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                          <  1:3570037/71‑1 (MQ=255)            aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                          <  1:3829845/71‑1 (MQ=255)             aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                     >  1:3214165/1‑59 (MQ=255)             aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                                                     >  1:107367/1‑59 (MQ=255)             aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                          <  1:79450/70‑1 (MQ=255)              aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                        <  1:1164842/71‑1 (MQ=255)              aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                        <  1:74153/71‑1 (MQ=255)               tCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                        >  1:2807331/1‑70 (MQ=255)                 gTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                       <  1:1274742/69‑1 (MQ=255)                    gCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                                                  >  1:3253011/1‑71 (MQ=255)                       ctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGaa                                                             >  1:2839880/1‑57 (MQ=255)                        tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                                                 >  1:2051136/1‑52 (MQ=255)                        tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                          >  1:1881350/1‑59 (MQ=255)                          aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                            >  1:2917685/1‑71 (MQ=255)                          aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                            >  1:1779617/1‑71 (MQ=255)                           aTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                           >  1:1634379/1‑71 (MQ=255)                           aTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                           >  1:1237321/1‑71 (MQ=255)                            ttGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                                                   <  1:1176248/46‑1 (MQ=255)                               cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                       <  1:2521743/71‑1 (MQ=255)                               cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                       <  1:3571031/71‑1 (MQ=255)                                 gagaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                        <  1:2008069/52‑1 (MQ=255)                                  agaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGaaa                                                            >  1:1054258/1‑47 (MQ=255)                                   gaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                                                  >  1:3158922/1‑40 (MQ=255)                                   gaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaaa                                                      >  1:2162911/1‑52 (MQ=255)                                   gaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                    >  1:3664975/1‑70 (MQ=255)                                    aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa                                  >  1:1007396/1‑71 (MQ=255)                                     aaTGGACAAACACATCTTTGCTGCCATCTTTCGGAGTGATGAAACCAAAACCTTTTTCAg                                            <  1:494509/60‑1 (MQ=25)                                      atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                <  1:987032/68‑1 (MQ=255)                                       tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                          >  1:232045/3‑60 (MQ=255)                                       tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                          >  1:1273661/3‑60 (MQ=255)                                       tAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                               <  1:2211814/71‑1 (MQ=255)                                        gcaCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAttt                               >  1:2548954/3‑70 (MQ=255)                                         gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                        >  1:2141863/1‑44 (MQ=255)                                         gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                        >  1:1701962/1‑44 (MQ=255)                                         caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAtttta                             <  1:658298/70‑3 (MQ=255)                                           cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                          <  1:2518770/40‑1 (MQ=255)                                           cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                            >  1:1534400/1‑70 (MQ=255)                                           cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACt                           <  1:1450160/71‑1 (MQ=255)                                           cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACt                           <  1:2728011/71‑1 (MQ=255)                                            aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACTa                          <  1:3525091/71‑1 (MQ=255)                                             aaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaa                         <  1:1430417/71‑1 (MQ=255)                                             aaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaa                         <  1:3226106/71‑1 (MQ=255)                                              acacATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                          >  1:2235155/1‑53 (MQ=255)                                              acacATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaaa                        >  1:3760251/1‑71 (MQ=255)                                                acgtCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAcc                      <  1:2092435/68‑1 (MQ=9)                                                acgtCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAc                       >  1:3505348/4‑70 (MQ=9)                                                acaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAcc                      <  1:3379269/71‑1 (MQ=255)                                                        gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt              >  1:1066236/1‑71 (MQ=255)                                                        gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt               <  1:2909648/70‑1 (MQ=255)                                                         gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtttt             <  1:2536224/69‑1 (MQ=255)                                                         gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtttt             <  1:3804793/69‑1 (MQ=255)                                                          tgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                               <  1:2182225/52‑1 (MQ=255)                                                          tgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAt                                 <  1:531124/50‑1 (MQ=255)                                                          agcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtttt             <  1:2864534/69‑1 (MQ=255)                                                            ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                            <  1:833754/53‑1 (MQ=255)                                                            ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt           <  1:2761567/70‑1 (MQ=255)                                                             cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                             <  1:3158695/35‑1 (MQ=255)                                                             cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaa                                      >  1:2866869/1‑42 (MQ=255)                                                             cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCa                             >  1:400185/1‑51 (MQ=255)                                                             cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa         <  1:1984849/71‑1 (MQ=255)                                                              aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                            <  1:2663073/35‑1 (MQ=255)                                                              aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                       <  1:2536476/56‑1 (MQ=255)                                                               tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                           <  1:1487346/35‑1 (MQ=255)                                                               tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                           <  1:1740353/35‑1 (MQ=255)                                                                cTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg            >  1:1881533/1‑65 (MQ=255)                                                                 tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt               >  1:3310296/1‑61 (MQ=255)                                                                 tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa     >  1:1521667/1‑71 (MQ=255)                                                                    cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa  >  1:1203435/1‑71 (MQ=255)                                                                    cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa  >  1:3118615/1‑71 (MQ=255)                                                                    |                                                                       TTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAA  >  minE/980169‑980306

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