Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,176,738	C→T	8.2%	T680M (ACG→ATG)	PP_2787 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,176,738	0	C	T	8.2%	173.1 / 3.1	61	T680M (ACG→ATG)	PP_2787	transporter
Reads supporting (aligned to +/- strand):  ref base C (32/23);  new base T (2/3);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTACGCCTGGGCATCCTCGGCATGCTGA  >  NC_002947/3176601‑3176875                                                                                                                                          |                                                                                                                                          tGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACatat                                                                                                                                <  3:199014/149‑1 (MQ=255)       ttACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCg                                                                                                                          >  7:175677/1‑149 (MQ=255)           tCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGgt                                                                                                                      >  8:158654/1‑149 (MQ=255)             gAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCggg                                                                                                                        <  4:345940/145‑1 (MQ=255)             gAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCggg                                                                                                                        >  3:345940/1‑145 (MQ=255)             gAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGgtgt                                                                                                                    >  1:278569/1‑149 (MQ=255)             gAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGgtgt                                                                                                                    >  4:125648/1‑149 (MQ=255)               aGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTg                                                                                                                  >  5:175984/1‑149 (MQ=255)                     tGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGaa                                                                                                            <  6:175984/149‑1 (MQ=255)                       aCGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAAt                                                                                                          <  7:10214/149‑1 (MQ=255)                         gCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCa                                                                                                        >  5:346705/1‑149 (MQ=255)                          cTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCAc                                                                                                       <  2:278569/149‑1 (MQ=255)                              tCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCa                                                                                                    <  7:158654/148‑1 (MQ=255)                                aCCATGACCTCCTCTGTGGCGGGAGATTACCAGGCGATGCACGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAggg                                                                                                 >  2:154717/1‑149 (MQ=255)                                  cATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACata                                                                                                                                 >  4:183703/1‑115 (MQ=255)                                  cATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACata                                                                                                                                 <  3:183703/115‑1 (MQ=255)                                  cATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGt                                                                                                >  6:203015/1‑148 (MQ=255)                                  cATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTc                                                                                               <  3:266694/149‑1 (MQ=255)                                                   cgggcgATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCa                                                                                                                >  6:18933/1‑115 (MQ=255)                                                   cgggcgATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCa                                                                                                                <  5:18933/115‑1 (MQ=255)                                                      gcgATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTg                                                                           <  8:175677/149‑1 (MQ=255)                                                              cAGGCGATGCTCGATAAGCTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCa                                                                                                        >  8:321283/1‑112 (MQ=255)                                                              cAGGCGATGCTCGATAAGCTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCa                                                                                                        <  7:321283/112‑1 (MQ=255)                                                              cAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATAAACAAACCGGCCAACCTGCAGGTCACACATGCCGGTACGCCGTCCATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGATCTTCTCGCTGGCGTTGCt                                                                   >  7:267338/1‑149 (MQ=255)                                                               aGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCt                                                                            >  6:35764/1‑139 (MQ=255)                                                               aGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCt                                                                            <  5:35764/139‑1 (MQ=255)                                                                       cTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGAc                                                          >  6:43481/1‑149 (MQ=255)                                                                        tCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGAACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGAcc                                                         >  5:20049/1‑149 (MQ=255)                                                                          taTAAACTACACGCCTGGCTGGCAATAAAAAAACCGGCCAAACTGCAGGTAACCCATGCCCGTACGCCGTACATATGGACCGGGGTGTTCCAGGAAATCACACAGTATCAGCTGCTGAGCTTATCTCTGGCGTTGATGGACg                                                              >  8:287674/2‑142 (MQ=255)                                                                          gATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCg                                                              <  7:287674/142‑1 (MQ=255)                                                                               aCTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGa                                                           >  7:310486/1‑140 (MQ=255)                                                                               aCTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGa                                                           <  8:310486/140‑1 (MQ=255)                                                                                cTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGc                                                                         >  4:256505/1‑125 (MQ=255)                                                                                cTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGc                                                                         <  3:256505/125‑1 (MQ=255)                                                                                 tAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGCCGTTGCTGGCCGTGACCTtgatgat                                                 >  5:246720/1‑148 (MQ=255)                                                                                      gCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGtt                                           >  4:307029/1‑149 (MQ=255)                                                                                        ctggctggCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCa                                                                                                    <  6:223099/90‑1 (MQ=255)                                                                                        ctggctggCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCa                                                                                                    >  5:223099/1‑90 (MQ=255)                                                                                           gctggcCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTATGCCGTACATATGGACCGGGGTGTTGCAGGAAATCCCCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTTCCCTTGATGATGATGGTCTGGc                                      >  4:348604/1‑149 (MQ=255)                                                                                           gctggcCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTgg                                       <  5:203015/148‑1 (MQ=255)                                                                                           gctggcCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTAAGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGc                                      >  4:342982/1‑149 (MQ=255)                                                                                              ggcCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTc                                    >  1:178107/1‑148 (MQ=255)                                                                                               gcCATCAACAAACCGGGCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTc                                    <  6:246720/147‑1 (MQ=255)                                                                                               gcCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCa                                   <  5:43481/148‑1 (MQ=255)                                                                                                 cATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGt                                >  3:195767/1‑149 (MQ=255)                                                                                                    caacaaACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAg                             >  4:166444/1‑149 (MQ=255)                                                                                                     aacaaACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAg                             >  3:247146/1‑148 (MQ=255)                                                                                                       caaACCGGCCAACTTGCAGGTCACCCATGCCGGTATGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTAc                          >  5:307838/1‑149 (MQ=255)                                                                                                        aaaCCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTACg                         >  6:160181/1‑149 (MQ=255)                                                                                                          aCCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGt                                >  4:245043/1‑140 (MQ=255)                                                                                                          aCCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGt                                <  3:245043/140‑1 (MQ=255)                                                                                                             ggCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTACGCCTgg                    <  1:154717/149‑1 (MQ=255)                                                                                                                 aaCCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGGGACCTTGATGATGATGTGCTGGCTCAAGTCCGTACGCCTGGGCAt                >  2:233243/1‑149 (MQ=255)                                                                                                                   ccTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTtgatgatg                                                <  6:387357/115‑1 (MQ=255)                                                                                                                   ccTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCGGGCCGTGACCTtgatgatg                                                >  5:387357/1‑115 (MQ=255)                                                                                                                    cTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCAtgatgat                                                 >  6:286150/1‑113 (MQ=255)                                                                                                                        aGGTCACCCATGCCGGTATGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTACGCCTGGGCATCCTCGGc         <  3:348604/149‑1 (MQ=255)                                                                                                                        aGGTCACCCATGCCGGTATGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTACGCCTGGGCATCCTCGGc         <  3:342982/149‑1 (MQ=255)                                                                                                                         ggTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTACGCCTGGGCATCCTCGGCa        <  3:125648/149‑1 (MQ=255)                                                                                                                           tCACCCATGCCGGTATGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTACGCCTGGGCATCCTCGGCATg      <  6:307838/149‑1 (MQ=255)                                                                                                                               ccATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTACGCCTGGGCATCCTCGGCATGCTGa  >  1:355350/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                          TGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCTGGCTCAAGTCAGTACGCCTGGGCATCCTCGGCATGCTGA  >  NC_002947/3176601‑3176875

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