Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,113,805	A→C	7.0%	*94Y (TAA→TAC)	hupN →	DNA‑binding protein HU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,113,805	0	A	C	7.0%	172.0 / 2.9	57	*94Y (TAA→TAC)	hupN	DNA‑binding protein HU
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base C (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

GCATCGGCAAGCTGAAAGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAGGCGACCATGCGGCTCTGTTTCTGCCCCTGGCCCA  >  NC_002947/1113666‑1113950                                                                                                                                            |                                                                                                                                                  gCATCGGCAAGCTGAAAGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCtt                                                                                                                                          <  7:87228/149‑1 (MQ=255) gCATCGGCAAGCTGAAAGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCt                                                                                                                                           >  3:19505/1‑148 (MQ=255) gCATCGGCAAGCTGAAAGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCt                                                                                                                                           <  4:19505/148‑1 (MQ=255)  cATCGGCAAGCTGAAAGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCttt                                                                                                                                         >  2:70198/1‑149 (MQ=255)        cAAGCTGAAAGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGaa                                                                                                                                   <  1:35228/149‑1 (MQ=255)           gCTGAAAGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAAttgtt                                                                                                                                              >  7:164199/1‑135 (MQ=255)           gCTGAAAGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAAttgtt                                                                                                                                              <  8:164199/135‑1 (MQ=255)                aaGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGcc                                                                                                                           <  4:18529/149‑1 (MQ=255)                    ctcCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGc                                                                                                                       >  1:101148/1‑149 (MQ=255)                      ccGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCtt                                                                                                                     >  2:106374/1‑149 (MQ=255)                       cGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCtt                                                                                                                     <  8:339724/148‑1 (MQ=255)                                tGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGgcgc                                                                                                           >  8:84795/1‑149 (MQ=255)                                  cccGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCgg                                                                                                         <  2:101148/149‑1 (MQ=255)                                     gTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGttt                                                                                                      >  1:216904/1‑149 (MQ=255)                                        ccGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGtttttt                                                                                                   <  1:106374/149‑1 (MQ=255)                                         cGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAATTACttgttgatg                                                                                                                                          <  5:164918/109‑4 (MQ=255)                                         cGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAATTACttgttgatg                                                                                                                                          >  6:164918/1‑106 (MQ=255)                                             cGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGc                                                                                              <  7:128884/149‑1 (MQ=255)                                                   ccTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCtg                                                                                           <  8:249804/146‑1 (MQ=255)                                                   ccTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCtg                                                                                           >  7:249804/1‑146 (MQ=255)                                                    cTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGAt                                                                                       >  7:113334/1‑149 (MQ=255)                                                     ttCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGAtt                                                                                      >  4:138926/1‑149 (MQ=255)                                                             gcgcTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCgg                                                                                                               >  3:31627/1‑116 (MQ=255)                                                             gcgcTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCgg                                                                                                               <  4:31627/116‑1 (MQ=255)                                                             gcgcTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTc                                                                                                  >  2:309929/1‑129 (MQ=255)                                                             gcgcTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTc                                                                                                  <  1:309929/129‑1 (MQ=255)                                                             gcgcTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTTTGATTTACGCTgg                                                                              <  3:138926/149‑1 (MQ=255)                                                             gcgcTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTg                                                                               <  5:236637/148‑1 (MQ=255)                                                              cgcTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCAGTGATTTACGCTggg                                                                             <  7:84795/149‑1 (MQ=255)                                                                   ccATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGc                                                                        >  5:330553/1‑149 (MQ=255)                                                                   ccATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACACCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGc                                                                        >  1:102212/1‑149 (MQ=255)                                                                                  ccAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAATTACttgttgatg                                                                                                                                          >  2:60743/1‑65 (MQ=255)                                                                                  ccAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAATTACttgttgatg                                                                                                                                          <  1:60743/68‑4 (MQ=39)                                                                                  ccAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTg                                                                               >  5:263991/1‑127 (MQ=255)                                                                                  ccAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTg                                                                               <  6:263991/127‑1 (MQ=255)                                                                                  ccAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTgc                                                          <  2:102212/148‑1 (MQ=255)                                                                                  ccAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTgc                                                          <  6:330553/148‑1 (MQ=255)                                                                                   cAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTgcgg                                                        >  2:104017/1‑149 (MQ=255)                                                                                    aaGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTgcggc                                                       >  8:220643/1‑149 (MQ=255)                                                                                      gAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGcgt                                                     >  1:104774/1‑149 (MQ=255)                                                                                             gtcAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGcgcg                                                                                                          >  7:123947/1‑89 (MQ=255)                                                                                             gtcAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGcgcg                                                                                                          <  8:123947/89‑1 (MQ=255)                                                                                             gtcAAGTTCGTTCCACCAAAAGTGCTCACCGACGCCATCAACAAGTAATTGTTGCTTTGACTAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCaa                                               <  1:70198/148‑1 (MQ=255)                                                                                                       ttcCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAg                                    >  7:128351/1‑149 (MQ=255)                                                                                                            gCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAGGCGAc                               <  7:220643/149‑1 (MQ=255)                                                                                                              cAAAGTGCTGACAGACGACAACAACAACTAATTGTTGCTTTGACGAAAACGCGCCCGGCTTTACCTTGCCCGTTTATTTCATGCCAGTTAATTACTCTTGGCGTGCCACGCTTGGCTTGCCTCCTCGACAAGGAAGCTCCAggacaaca                             >  4:253664/1‑143 (MQ=255)                                                                                                                          cGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAGGCGACCATGCGGCTCTGtt                 <  8:113334/149‑1 (MQ=255)                                                                                                                           gACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAGGCGACCATGCGGCTCTGttt                <  2:233981/149‑1 (MQ=255)                                                                                                                              gCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAGGCGACCATGCGGCTCTGTTTCTg             >  4:275495/1‑149 (MQ=255)                                                                                                                                 aTCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTgcggcg                                                      <  2:114337/105‑1 (MQ=255)                                                                                                                                  tCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTTGGCGTGCCAGGCTTGGCGTTCGGCGTCGTCAAGGAAGCTCCAGGCGACCATGCGGCTCTGTTTCTGcccc         >  5:44373/1‑149 (MQ=255)                                                                                                                                   caacaaGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCAAGCCTGTCATTTACGATGGGCGTGCCAGGCTTGGCGTGCGGCGTCGACAAGGAAGCTCCAGGCGACCATGCGGCTCTGTATCTGCCACt        >  8:28998/1‑149 (MQ=255)                                                                                                                                    aactaGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAGGCGACCATGCGGCTCTGTTTCTGCCCCTg       <  7:28998/149‑1 (MQ=255)                                                                                                                                    aacaaGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAGGCGACCATGCGGCTCTGTTTCTGCCCCTg       <  6:44373/149‑1 (MQ=255)                                                                                                                                         gTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAg                                          <  5:299798/109‑1 (MQ=255)                                                                                                                                         gTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAg                                          >  6:299798/1‑109 (MQ=255)                                                                                                                                         gTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAGGCGACCATGCGGCTCTGTTTCTTGCCCTGGGCCa  >  4:336240/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                  GCATCGGCAAGCTGAAAGTCTCCGAGCGTCCTGCCCGTACCGGCCGCAACCCTTCGACTGGCGCTGCCATCGAAATCGCTGCCAAGAAAGTCGTCAAGTTCGTTCCAGCCAAAGTGCTGACCGACGCCATCAACAAGTAATTGTTGCTTTGACGAAACCGCGCCCGGCTTGCCCGGGCGCGGTTTTTTCATGCCTGTGATTTACGCTGGGCGTGCCAGGCTTGGCGTGCGGCGTCGTCAAGGAAGCTCCAGGCGACCATGCGGCTCTGTTTCTGCCCCTGGCCCA  >  NC_002947/1113666‑1113950

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