Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	980,296	G→A	5.9%	N3N (AAC→AAT)	cspI ←	cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	980,296	0	G	A	5.9%	280.5 / 4.2	84	N3N (AAC→AAT)	cspI	cold shock protein
Reads supporting (aligned to +/- strand):  ref base G (39/40);  new base A (3/2);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAG  >  minE/980228‑980363                                                                     |                                                                    ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                                                    <  1:1421616/70‑1 (MQ=255)      tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa                                                                  >  1:2373808/1‑67 (MQ=255)        tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGa                                                                <  1:282312/67‑1 (MQ=255)        tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                                                             <  1:1702742/70‑1 (MQ=255)        tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                                                             <  1:1899011/70‑1 (MQ=255)        tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                                                             <  1:1131223/70‑1 (MQ=255)         cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGa                                                                <  1:2707690/66‑1 (MQ=255)          ggCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat                                                          <  1:3668892/71‑1 (MQ=255)           gagTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCATTTTATTAGACATAttt                                                         <  1:3732265/69‑1 (MQ=255)           gCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat                                                          >  1:1403659/1‑70 (MQ=255)           gCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat                                                          <  1:2375025/70‑1 (MQ=255)              tgatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                                                    <  1:1443865/57‑1 (MQ=255)               gatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAg                                                                 <  1:141036/59‑1 (MQ=255)                atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAACCCAGTCATTTTGTTAGACATAATTATTAc                                                    >  1:1118278/1‑71 (MQ=255)                atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTAc                                                    >  1:2435792/1‑71 (MQ=255)                atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTAc                                                    >  1:3689580/1‑71 (MQ=255)                atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTAc                                                    >  1:1791914/1‑71 (MQ=255)                atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTAc                                                    >  1:681991/1‑71 (MQ=255)                   aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                                                    >  1:3279076/1‑52 (MQ=255)                   aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtt                                                 <  1:3966583/71‑1 (MQ=255)                    aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAc                                                               >  1:3125964/1‑56 (MQ=255)                    aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCttt                                                <  1:2234344/71‑1 (MQ=255)                     aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAc                                                               >  1:3450140/1‑55 (MQ=255)                     aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAc                                                               >  1:1521040/1‑55 (MQ=255)                     aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAGATAATTATTACCtttt                                               >  1:3562023/1‑71 (MQ=255)                     aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt                                                         <  1:275604/61‑1 (MQ=255)                     aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtttt                                               <  1:1097479/71‑1 (MQ=255)                     aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtttt                                               <  1:2667077/71‑1 (MQ=255)                     aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtttt                                               <  1:2858166/71‑1 (MQ=255)                         aaaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                                                                     <  1:1747438/45‑1 (MQ=255)                          aaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaaga                                          <  1:1520341/71‑1 (MQ=255)                          aaCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaaga                                          <  1:2374128/71‑1 (MQ=255)                           aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                                                        >  1:3751185/1‑56 (MQ=255)                           aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaagaa                                         <  1:224202/71‑1 (MQ=255)                           aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaaga                                          >  1:2403280/1‑70 (MQ=255)                           aCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaag                                           >  1:2613212/1‑69 (MQ=255)                            ccTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                                                        >  1:1797803/1‑55 (MQ=255)                             cTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAt                                       <  1:3350055/71‑1 (MQ=255)                              tttttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                                                                    <  1:3269716/41‑1 (MQ=255)                              tttttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaag                                           <  1:681908/66‑1 (MQ=255)                              tttatCAGCGTTAAACCATTTTACTAAACCAGTCATTTTATTTGa                                                                <  1:4067809/45‑1 (MQ=255)                               ttttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAt                                       >  1:1281556/1‑69 (MQ=255)                                tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAAtt                                      <  1:3362705/69‑1 (MQ=255)                                tttCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAg                                    >  1:909849/1‑71 (MQ=255)                                tatCAGCGTTAAACCATTTTACTAAACCAGTCATTTTATTTGACa                                                              >  1:2393450/3‑45 (MQ=255)                                tatCAGCGTTAAACCATTTTACTAAACCAGTCATTTTATTTGACa                                                              >  1:2297928/3‑45 (MQ=255)                                tatCAGCGTTAAACCATTTTACTAAACCAGTCATTTTATTTGACa                                                              >  1:3999664/3‑45 (MQ=255)                                    aGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                                                             >  1:3624282/1‑42 (MQ=255)                                    aGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                                                             >  1:3894000/1‑42 (MQ=255)                                    aGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                                                             >  1:4118427/1‑42 (MQ=255)                                      ggTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat                                                          <  1:1486013/43‑1 (MQ=255)                                      ggTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa                                        >  1:1631578/1‑61 (MQ=255)                                        ttAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTga                                             <  1:2276828/54‑1 (MQ=255)                                        ttAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                              <  1:1074534/53‑1 (MQ=255)                                         tAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGaaa                                        <  1:2924969/58‑1 (MQ=255)                                          aaaCCATTTCACTAAACCAGTCATTTTGTTAGACATaa                                                           >  1:1884872/1‑38 (MQ=255)                                          aaaCCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                                                        >  1:313392/1‑41 (MQ=255)                                          aaaCCATTTCACTAAACCAGTCATTTTGTTAGACATAatta                                                        >  1:1082944/1‑41 (MQ=255)                                          aaaCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGcc                                  <  1:3582828/63‑1 (MQ=255)                                          aaaCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTggg                            >  1:1012658/1‑69 (MQ=255)                                          aaaCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTggg                            >  1:3031080/1‑69 (MQ=255)                                            aCCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCttt                                                >  1:1177447/1‑47 (MQ=255)                                             ccATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGaa                       <  1:313743/71‑1 (MQ=255)                                                  tCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTc                  >  1:1406659/1‑71 (MQ=255)                                                  tCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATCGTc                  >  1:3376049/1‑71 (MQ=255)                                                   cACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                              <  1:3570927/42‑1 (MQ=255)                                                    acaaaaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGa                        <  1:2762640/60‑1 (MQ=255)                                                    aCTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTcc                 <  1:2045876/70‑1 (MQ=255)                                                     cTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaa                                            >  1:548789/1‑42 (MQ=255)                                                     cTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTcc                 <  1:4045292/69‑1 (MQ=255)                                                        aaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATg                     >  1:1671788/1‑62 (MQ=255)                                                        aaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGaaa             >  1:950381/1‑70 (MQ=255)                                                        aaCCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGaaa             >  1:1709010/1‑70 (MQ=255)                                                          ccAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGaaaaaa          <  1:4025640/71‑1 (MQ=255)                                                          ccAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGaaaaaa          <  1:2628239/71‑1 (MQ=255)                                                           cAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGc                                   <  1:991627/45‑1 (MQ=255)                                                                aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCt                                >  1:2656199/1‑43 (MQ=255)                                                                aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGaaaaaaaa        >  1:3972381/1‑67 (MQ=255)                                                                aTTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAatat     >  1:4098067/1‑70 (MQ=255)                                                                 ttttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCa   <  1:903030/71‑1 (MQ=255)                                                                 ttttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCa   <  1:307229/71‑1 (MQ=255)                                                                  tttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGa                        >  1:470547/1‑49 (MQ=255)                                                                   ttGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGaa                       >  1:4102079/1‑49 (MQ=255)                                                                    tGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAata      <  1:9289/65‑1 (MQ=255)                                                                    tGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCag  <  1:3072719/69‑1 (MQ=255)                                                                     |                                                                    CCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAG  >  minE/980228‑980363

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