Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,434,228	G→A	7.1%	D72N (GAC→AAC)	PP_3929 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,434,228	0	G	A	7.1%	161.7 / 2.1	56	D72N (GAC→AAC)	PP_3929	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base A (2/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TATCGTGCTCACGCGTCTTGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATCATGCGTCTGACTCTGCCTTCCC  >  NC_002947/4434092‑4434376                                                                                                                                           |                                                                                                                                                     tatCGTGCTCACGCGTCTTGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGcgacgac                                                                                                                                           <  6:353303/148‑1 (MQ=255)     gTGCTCACGCGTCTTGCAAGCGTTTCTCTGGTTGCTGTTGTCCTGCTGAGCGCAGGTTGCCATAGCCATCGCTACCATGATGAGGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACgacgatga                                                                                                                                      <  6:98579/149‑1 (MQ=255)      tGCTCACGCGTCTTGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACgacgatgac                                                                                                                                     >  2:297688/1‑149 (MQ=255)           aCGCGTCTTGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGCCCGCGACGACGATGACGATCg                                                                                                                                >  7:78967/1‑149 (MQ=255)           aCGCGTCTTGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATc                                                                                                                                 <  5:282982/148‑1 (MQ=255)               gTCTTGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATc                                                                                                                                 <  8:286298/144‑1 (MQ=255)               gTCTTGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATc                                                                                                                                 >  7:286298/1‑144 (MQ=255)                 cTTGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATAACCGCGACGACGATGACGATCGGCAATa                                                                                                                          >  5:18190/1‑149 (MQ=255)                   tGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATAcc                                                                                                                        >  3:87693/1‑149 (MQ=255)                      aaGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACgacgatg                                                                                                                                       >  5:361148/1‑131 (MQ=255)                      aaGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACgacgatg                                                                                                                                       <  6:361148/131‑1 (MQ=255)                      aaGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACcgcg                                                                                                                     >  8:32658/1‑149 (MQ=255)                               tctGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGc                                                                                                             >  3:89967/1‑148 (MQ=255)                                        tgttgtCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGc                                                                                                    >  2:361407/1‑148 (MQ=255)                                          ttgtCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGc                                                                                                    <  1:361407/146‑1 (MQ=255)                                          ttgtCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGa                                                                                                 <  4:87693/149‑1 (MQ=255)                                            gtCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGaa                                                                                                <  4:281184/148‑1 (MQ=255)                                            gtCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATAACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGaa                                                                                                <  6:18190/148‑1 (MQ=255)                                             tCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAg                                                                                               <  2:288949/148‑1 (MQ=255)                                               ctgctgACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGt                                                                                            >  6:226699/1‑149 (MQ=255)                                                  atgACCGCAGGTTTCCATAGCCATCGCTACCATGATGACGACGGCGGGCGTGTCAAGGCGACTAAAATCTTCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTc                                                                                           <  1:297688/148‑1 (MQ=255)                                                      ccGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCa                                                                                     <  2:252549/149‑1 (MQ=255)                                                         cAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGcc                                                                                  <  2:344478/149‑1 (MQ=255)                                                            gTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCt                                                                                >  5:107783/1‑148 (MQ=255)                                                                ccAAAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATAACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCt                                                                           >  1:315883/1‑149 (MQ=255)                                                                  aTAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTgc                                                                         >  5:116838/1‑149 (MQ=255)                                                                   tAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTgcc                                                                        <  4:89967/149‑1 (MQ=255)                                                                         tCGCTACCATGATGACGACGGCTGGTGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGt                                                                  >  5:8628/1‑149 (MQ=255)                                                                           gCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTg                                                                <  7:32658/149‑1 (MQ=255)                                                                             tACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACAATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTgcg                                                              >  2:31024/1‑149 (MQ=255)                                                                                      gacgacGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGAc                                                       <  6:160375/147‑1 (MQ=255)                                                                                      gacgacGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGAc                                                       >  5:160375/1‑147 (MQ=255)                                                                                      gacgacGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACt                                                      >  6:339654/1‑148 (MQ=255)                                                                                       acgacgGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTa                                                    >  2:238493/1‑149 (MQ=255)                                                                                         gacgGCTGGCGTG‑‑AACGCGACCGTGACCATCGCCACCACCACCGCCATGACCTCGACGACGATGACGATCGGCAATACCGCGTTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATc                                                  <  5:150645/149‑1 (MQ=255)                                                                                             gCTGGCGTG‑‑AACGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCtga                                              <  8:78967/149‑1 (MQ=255)                                                                                             gCTGGCGGG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCtga                                              <  6:228032/149‑1 (MQ=255)                                                                                               tGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCtgat                                             >  3:7610/1‑148 (MQ=255)                                                                                                 gCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGtt                                                                 >  3:168486/1‑126 (MQ=255)                                                                                                 gCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGtt                                                                 <  4:168486/126‑1 (MQ=255)                                                                                                  cGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGAtt                                         >  3:335749/1‑149 (MQ=255)                                                                                                         ccGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATTATTCCttt                                    <  4:7610/149‑1 (MQ=255)                                                                                                           gcgACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGa                                  <  1:31024/149‑1 (MQ=255)                                                                                                           gcgACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGa                                  <  6:107783/149‑1 (MQ=255)                                                                                                           gcgACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGa                                  <  6:116838/149‑1 (MQ=255)                                                                                                                 gTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTc                             >  8:334283/1‑148 (MQ=255)                                                                                                                     ccATCGCCACCACCACCGCCATGAGCGCGACGACGATGACGATCGGAAATACCGCGTTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCTCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATc                        <  5:226699/149‑1 (MQ=255)                                                                                                                     ccATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATc                        <  1:238493/149‑1 (MQ=255)                                                                                                                       aTCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATCAt                      <  5:339654/149‑1 (MQ=255)                                                                                                                          gCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATCATGc                    <  7:334283/148‑1 (MQ=255)                                                                                                                           ccaccaccaCCGCCATAACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATCATGCGt                  <  2:315883/149‑1 (MQ=255)                                                                                                                            caccaccacCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCttt                                    >  3:289638/1‑130 (MQ=255)                                                                                                                            caccaccacCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCttt                                    <  4:289638/130‑1 (MQ=255)                                                                                                                              ccaccaccGCCATGACCGCGACGACGATCACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATCATGCGTCt                >  1:204572/1‑148 (MQ=255)                                                                                                                               caccaccGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATCATGCGTCTGa              >  5:130980/1‑149 (MQ=255)                                                                                                                                  caccGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTG‑AGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATCATGCGTCTGActct          >  3:240139/1‑149 (MQ=255)                                                                                                                                           gACCGCGACGACGATGACGATCGCCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATCATGCGTCTGACTCTGCCTTccc  >  3:104155/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                     TATCGTGCTCACGCGTCTTGCAAGCCTTTCTCTGGTTGCTGTTGTCCTGCTGACCGCAGGTTGCCATAGCCATCGCTACCATGATGACGACGGCTGGCGTG‑‑ACCGCGACCGTGACCATCGCCACCACCACCGCCATGACCGCGACGACGATGACGATCGGCAATACCGCGGTGACCGCTACTACCGCTGAAGCGTCCTTTCAGCCCTGACCTGCCGCCCGTTGCGCGGCGACTTATCCTGATGATTCCTTTGAGGTTCGTATCATGCGTCTGACTCTGCCTTCCC  >  NC_002947/4434092‑4434376

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