Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,456,137	T→G	13.8%	K109T (AAG→ACG)	eco ←	ecotin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,456,137	0	T	G	13.8%	127.1 / 15.0	51	K109T (AAG→ACG)	eco	ecotin
Reads supporting (aligned to +/- strand):  ref base T (26/18);  new base G (5/2);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

CTTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGC  >  NC_002947/3455991‑3456285                                                                                                                                                   |                                                                                                                                                     ctTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTg                                                                                                                                                    <  2:169911/149‑1 (MQ=255)  tttAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCAAGAAGACCGAGGCCTTCGTGc                                                                                                                                                   >  2:84322/1‑149 (MQ=255)       ccTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCAAGAAGACCGAGGCCTTCGTGCCggtg                                                                                                                                              <  3:337876/149‑4 (MQ=255)         ttCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCgg                                                                                                                                            >  3:148628/1‑149 (MQ=255)               acaTCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGc                                                                                                                                      >  6:12375/1‑149 (MQ=255)                  tCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCAt                                                                                                                                   >  7:279247/1‑149 (MQ=255)                  tCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCAt                                                                                                                                   >  7:225494/1‑149 (MQ=255)                                 cacaCGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTggccgg                                                                                                                    >  5:8584/1‑149 (MQ=255)                                     cgcgGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTggccgggccg                                                                                                                <  3:79090/149‑1 (MQ=255)                                          tAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGAc                                                                                                           <  6:8584/149‑1 (MQ=255)                                                aCTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCAAGAAGACCGAGGCCTTCGTGCCgg                                                                                                                                                >  4:276780/1‑105 (MQ=38)                                                aCTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCAAGAAGACCGAGGCCTTCGTGCCgg                                                                                                                                                <  3:276780/106‑2 (MQ=38)                                                      aCATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCg                                                                                               <  8:279247/149‑1 (MQ=255)                                                      aCATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCg                                                                                               <  8:225494/149‑1 (MQ=255)                                                            ttGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCAATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGgtagt                                                                                          >  7:3732/1‑148 (MQ=255)                                                                    aTACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGc                                                                                 >  8:109796/1‑149 (MQ=255)                                                                      acacCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGccc                                                                               >  5:32208/1‑149 (MQ=255)                                                                      acacAACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGccc                                                                               >  8:188860/1‑149 (MQ=255)                                                                        accaccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGccc                                                                               >  6:219666/1‑147 (MQ=255)                                                                        accaccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGccc                                                                               >  6:219663/1‑147 (MQ=255)                                                                        accaccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCa                                                                             <  5:219663/149‑1 (MQ=255)                                                                        accaccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCa                                                                             <  5:219666/149‑1 (MQ=255)                                                                            ccaccGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCt                                                                         >  8:91002/1‑149 (MQ=255)                                                                                  gCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAgg                                                                   >  7:71724/1‑149 (MQ=255)                                                                                         gCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTg                                                                                                                            <  4:40380/85‑1 (MQ=255)                                                                                         gCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTg                                                                                                                            >  3:40380/1‑85 (MQ=255)                                                                                         gCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGgtgtgt                                                             <  7:188860/148‑1 (MQ=255)                                                                                           tgttgtAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGc                                                           >  3:44659/1‑148 (MQ=255)                                                                                                 aGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAgg                                                    >  8:186329/1‑149 (MQ=255)                                                                                                   cgcAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGtt                                                  <  5:12375/149‑1 (MQ=255)                                                                                                    gcagcaGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTg                                                 <  8:3732/149‑1 (MQ=255)                                                                                                    gcagcaGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTg                                                 <  7:109796/149‑1 (MQ=255)                                                                                                    gcagcaGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCAGACGATGTACACGCGGTAGTAATTGAAGCACCAGCCTTCCAGGGTGTGTGCTTCCAGGTTg                                                 >  2:272318/1‑149 (MQ=255)                                                                                                            aaGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCCGCGTGCTGGCCGGGCCGCTGACGTTGTCCCGGCGGTAGGAAGTGTAGCCCCAGCCTTCCAGGGTGTGTGGTTCCAGGTTGCCGccccg                                         >  8:216848/1‑147 (MQ=255)                                                                                                              gCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTggccggg                                                                                                                   >  8:116428/1‑73 (MQ=255)                                                                                                                cATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGgcgc                                     >  5:307295/1‑149 (MQ=255)                                                                                                                    gCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGAc                                 <  6:307295/149‑1 (MQ=255)                                                                                                                     ccGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGAc                                 >  7:210342/1‑148 (MQ=255)                                                                                                                     ccGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACg                                <  6:32208/149‑1 (MQ=255)                                                                                                                      cGACCACCGGCAGGAAGGCCTCGGTCTTCTTGCCGTCCGGTCAGGCCATCAGCGTGCTGGCCGGGGCGCTGAGGTTGTCCAGGCGGTAGTTATTGTAGCCCCAGCCTTGCTGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGa                               <  1:272318/149‑1 (MQ=255)                                                                                                                      cGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGAAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGa                               <  4:148628/149‑1 (MQ=255)                                                                                                                       gACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGAt                              >  8:259249/1‑149 (MQ=255)                                                                                                                         cggacGGCAAGAAGACCGAGGCCTTCGTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTg                            >  3:42377/5‑149 (MQ=255)                                                                                                                         cggacGGCAAGAAGACCGAGGCCTTCGTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTg                            >  3:218950/5‑149 (MQ=255)                                                                                                                            ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGgc                                       <  7:178248/135‑1 (MQ=255)                                                                                                                            ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGgc                                       >  8:178248/1‑135 (MQ=255)                                                                                                                            ccGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAg                         >  5:317587/1‑149 (MQ=255)                                                                                                                             cGGCAAGAAGACCGAGGCCTTCGTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGt                        >  2:221966/1‑149 (MQ=255)                                                                                                                                   gAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCt                  >  1:328314/1‑149 (MQ=255)                                                                                                                                        ccTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGc             <  8:210342/149‑1 (MQ=255)                                                                                                                                         cTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCg            >  3:78875/1‑149 (MQ=255)                                                                                                                                                ttcttGCCGTCCGGGCAGGCCATCAGCGTGATGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGc      >  1:205104/1‑148 (MQ=255)                                                                                                                                                   ttGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGc  >  6:71351/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CTTTAGCCTTCTGCACATCCTGCGATGCCGACCACACGCGGTAGCGCACTTCGACATCCTTGGGCACATACACCACCACCGGCAGCTTGCTGTTGTAGCGCAGCAGGAAGCCATCGCCGACCACCGGCACGAAGGCCTCGGTCTTCTTGCCGTCCGGGCAGGCCATCAGCGTGCTGGCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCAGCGTCTTGCCAGC  >  NC_002947/3455991‑3456285

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