Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,190,621	G→A	7.1%	Q265* (CAG→TAG)	cysK ←	cysteine synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,190,621	0	G	A	7.1%	162.1 / 2.5	56	Q265* (CAG→TAG)	cysK	cysteine synthase A
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base A (2/2);  total (27/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 = 5.39e-01

CAGCATGCTCGACAAGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCTGGATTTTGTGCGGGCTTGGCTTGAGTTCTTCGCCG  >  NC_002947/5190478‑5190763                                                                                                                                                |                                                                                                                                                cAGCATGCTCGACAAGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCa                                                                                                                                            <  2:340633/149‑1 (MQ=255)  aGCATGCTCGACAAGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAg                                                                                                                                           >  4:30942/1‑149 (MQ=255)      tGCTCGACAAGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGa                                                                                                                                       <  3:30942/149‑1 (MQ=255)         tCGACAAGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGACCTCTTCCTGCATCAGGCGaa                                                                                                                                      >  7:273576/1‑147 (MQ=255)          cGACAAGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCg                                                                                                                                   >  6:286510/1‑149 (MQ=255)           gACAAGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCg                                                                                                                                   >  7:189127/1‑148 (MQ=255)           gACAAGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGc                                                                                                                                  >  7:183898/1‑149 (MQ=255)               aGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATg                                                                                                                              >  5:263117/1‑149 (MQ=255)                    cgctcgcCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCtt                                                                                                                         <  4:188200/149‑1 (MQ=255)                     gctcgcCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTg                                                                                                                        <  8:189127/149‑1 (MQ=255)                     gctcgcCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTg                                                                                                                        <  8:183898/149‑1 (MQ=255)                     gctcgcCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCGCTTCCTGCATCAGGCGAATCGCCATGGCCtt                                                                                                                         <  3:305140/148‑1 (MQ=255)                                  gggCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTtcgtcg                                                                                                            <  2:351215/148‑1 (MQ=255)                                  gggCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCgg                                                                                                           <  7:290968/149‑1 (MQ=255)                                        gATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCg                                                                                                     >  2:92848/1‑149 (MQ=255)                                                 gaTGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACtt                                                                                           <  8:7045/149‑1 (MQ=255)                                                 gaTGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTACATCAGGCGATTCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACtt                                                                                           >  2:124790/1‑149 (MQ=255)                                                          aCCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACa                                                                                   <  8:133580/148‑1 (MQ=255)                                                             cTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATc                                                                                >  4:92984/1‑148 (MQ=255)                                                              tGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGa                                                                              >  4:228829/1‑149 (MQ=255)                                                                cATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGAc                                                                             <  1:92848/148‑1 (MQ=255)                                                                 aTTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAg                                                                           <  8:63068/149‑1 (MQ=255)                                                                    tCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTc                                                                        >  1:222285/1‑149 (MQ=255)                                                                        gTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAgg                                                                    >  3:333585/1‑149 (MQ=255)                                                                           tttCGGCCAGACGCACCGCTGCCGCCATTGCAGCACAGCAGGAAATACCAAAAAGGATGCACTCATCCTGCATCAGGCGAATCGCCATGGCCTTCGACTCTTAGACGGTCACCGTCTC‑CACTTGGTAGACAATCGACAGGTCGAGGttt                                                                 >  6:177553/1‑149 (MQ=255)                                                                            ttCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGtttt                                                                >  6:229426/1‑149 (MQ=255)                                                                             tCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGtttt                                                                <  5:286510/148‑1 (MQ=255)                                                                              cGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGttttt                                                               <  3:228829/148‑1 (MQ=255)                                                                               ggCCAGACGCACCGCTGCCGCCATTGCCGGACCGCAGGAAATACCACAAAGGATGCCCTCTTCGTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGGCGAGGTTTTTCg                                                             <  5:177958/149‑1 (MQ=255)                                                                               ggCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTc                                                              <  3:92984/148‑1 (MQ=255)                                                                               ggCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTACATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTc                                                              <  1:124790/148‑1 (MQ=255)                                                                                       gCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTACATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGaaa                                                     >  8:301664/1‑149 (MQ=255)                                                                                         aCCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCg                                                             <  8:155301/139‑1 (MQ=255)                                                                                         aCCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCg                                                             >  7:155301/1‑139 (MQ=255)                                                                                            gCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTACATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGgc                                                <  7:301664/149‑1 (MQ=255)                                                                                                ccgccATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGcc                                             >  6:7604/1‑148 (MQ=255)                                                                                                      ttGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGccc                                      <  2:222285/149‑1 (MQ=255)                                                                                                      ttGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGCCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGccc                                      <  4:342348/149‑1 (MQ=255)                                                                                                        gCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCt                                     <  5:229426/148‑1 (MQ=255)                                                                                                             accgcaGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATc                                                                                <  6:253222/100‑1 (MQ=255)                                                                                                             accgcaGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATc                                                                                >  5:253222/1‑100 (MQ=255)                                                                                                             accgaaGCAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCGTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCTGGAttt                               >  6:284686/1‑149 (MQ=255)                                                                                                               cgcaGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCTGGATTTtg                             >  2:44607/1‑149 (MQ=255)                                                                                                                 caGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCTGGATTTtgtg                           >  4:216292/1‑149 (MQ=255)                                                                                                                          ccACAACCGATGCCCTCACCCTGCATCAGGCGAATCGCCAGGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCTGGATTTTGTGCGggcttgg                  <  6:263117/149‑1 (MQ=255)                                                                                                                           cacaAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCTGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATTCCCTGGATTTTGTGCGggcttggc                 <  5:177553/149‑1 (MQ=255)                                                                                                                           cacaAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGgcgc                                              >  6:199042/1‑120 (MQ=255)                                                                                                                           cacaAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGgcgc                                              <  5:199042/120‑1 (MQ=255)                                                                                                                           cacaAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCTGGATTTTGTGCGggcttggc                 >  7:310564/1‑149 (MQ=255)                                                                                                                                aGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCTGGATTTTGTGCGGGCTTGGCTTGAg            <  8:46824/149‑1 (MQ=255)                                                                                                                                  gATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGttttt                                                               <  4:77342/96‑1 (MQ=255)                                                                                                                                  gATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGttttt                                                               >  3:77342/1‑96 (MQ=255)                                                                                                                                       ccTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATg                                         <  3:131407/113‑1 (MQ=255)                                                                                                                                       ccTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATg                                         >  4:131407/1‑113 (MQ=255)                                                                                                                                        ctctTCCTGCATCAGGCGAATAGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTGGACAATCGAAAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGCTGCCCTGGATTGTGTGCGGGCTTGTCTTGAGTTCTTCGc    >  6:243452/1‑149 (MQ=255)                                                                                                                                           ttCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTCTCACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCTGGATTTTGTGCGGGCTTGGCTTGAGTTCTTCGCCg  >  1:94541/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                CAGCATGCTCGACAAGTATCGCTCGCCCGAATCGGGCAGGATTACGACGATGGTCTTACCCTGCATTTCCGGTTTTTCGGCCAGACGCACCGCTGCCGCCATTGCCGCACCGCAGGAAATACCACAAAGGATGCCCTCTTCCTGCATCAGGCGAATCGCCATGGCCTTGGACTCTTCGTCGGTCACCGTCTC‑CACTTGGTCGACAATCGACAGGTCGAGGTTTTTCGGCACGAAACCGGCGCCGATGCCCTGGATTTTGTGCGGGCTTGGCTTGAGTTCTTCGCCG  >  NC_002947/5190478‑5190763

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