Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,044,240	A→G	9.1%	S350S (AGT→AGC)	PP_5297 ←	amino acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,044,240	0	A	G	9.1%	150.2 / 5.2	55	S350S (AGT→AGC)	PP_5297	amino acid transporter
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base G (3/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

CGGCCAGGCCGGCGCAGATGGTGGCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAGGGTGGCGAT  >  NC_002947/6044092‑6044383                                                                                                                                                     |                                                                                                                                                cGGCCAGGCCGGCGCAGATGGTGGCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTa                                                                                                                                                 <  4:265366/149‑1 (MQ=255)  ggCCAGGCCGGCGCAGATGGTGGCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAATAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTAc                                                                                                                                                <  1:236825/149‑1 (MQ=255)   gCCAGGCCGGCGCAGATGGTGGCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACt                                                                                                                                               >  4:196576/1‑149 (MQ=255)                  aTGGTGGCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGAt                                                                                                                                <  2:55688/149‑1 (MQ=255)                  aTGGTGGCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGAt                                                                                                                                <  6:194156/149‑1 (MQ=255)                    ggtggCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTGCTGGCCAGTACGGCgg                                                                                                                                 <  2:72407/146‑2 (MQ=255)                    ggtggCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTGCTGGCCAGTACGGCgg                                                                                                                                 >  1:72407/1‑145 (MQ=255)                    ggtggCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATg                                                                                                                               <  2:208496/148‑1 (MQ=255)                                     tttGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGTGATGGACATCAg                                                                                                                       <  1:277518/139‑1 (MQ=255)                                     tttGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGAGATGGACATCAg                                                                                                                       >  2:277518/1‑139 (MQ=255)                                       tGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCg                                                                                                            <  3:196576/148‑1 (MQ=255)                                                gggCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTGCTGGCCAGTACGGCg                                                                                                                                  >  4:57448/1‑117 (MQ=255)                                                gggCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTGCTGGCCAGTACGGCg                                                                                                                                  <  3:57448/117‑1 (MQ=255)                                                gggCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAAACCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGAGACGCCGCCGAACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGAc                                                                                                  >  6:155411/1‑149 (MQ=255)                                                    aGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCggg                                                                                              <  2:150143/149‑1 (MQ=255)                                                    aGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCgg                                                                                               <  8:135014/148‑1 (MQ=255)                                                    aGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCgg                                                                                               >  7:135014/1‑148 (MQ=255)                                                      acacTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGggcc                                                                                            <  8:110878/149‑1 (MQ=255)                                                         cTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGggccggc                                                                                         <  1:225868/149‑1 (MQ=255)                                                          tGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATc                                                                                                                   <  6:64183/122‑1 (MQ=255)                                                          tGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATc                                                                                                                   >  5:64183/1‑122 (MQ=255)                                                            gCGAACTTCTGCGGTACTGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAg                                                                                      >  3:57356/1‑149 (MQ=255)                                                              gAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGt                                                                                    >  2:1384/1‑149 (MQ=255)                                                                  ttCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTggg                                                                                 >  3:177331/1‑148 (MQ=255)                                                                       cGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGa                                                                           >  7:14399/1‑149 (MQ=255)                                                                               cACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTcc                                                                   >  4:18236/1‑149 (MQ=255)                                                                               cACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTcc                                                                   <  8:14399/149‑1 (MQ=255)                                                                                  cGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTg                                                                 >  3:262376/1‑148 (MQ=255)                                                                                    taataaCCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGAt                                                              <  3:267109/149‑1 (MQ=255)                                                                                           ccATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCggg                                                       >  8:11364/1‑149 (MQ=255)                                                                                           ccATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCggg                                                       >  1:148423/1‑149 (MQ=255)                                                                                           ccATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATAGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCggg                                                       >  7:192420/1‑149 (MQ=255)                                                                                             aTGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCg                                                                                                >  8:144610/1‑106 (MQ=255)                                                                                             aTGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCg                                                                                                <  7:144610/106‑1 (MQ=255)                                                                                             aTGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCaaa                                                                             <  5:331526/125‑1 (MQ=255)                                                                                             aTGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCaaa                                                                             >  6:331526/1‑125 (MQ=255)                                                                                                gCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGgttgtt                                                  <  2:148423/149‑1 (MQ=255)                                                                                                gCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGgttgtt                                                  <  1:1384/149‑1 (MQ=255)                                                                                                 ccAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGgttgtt                                                  >  5:178479/1‑148 (MQ=255)                                                                                                                ccGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCggt                                  >  8:144031/1‑147 (MQ=255)                                                                                                                    cgcCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCAc                              <  4:177331/149‑1 (MQ=255)                                                                                                                     gcCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCCGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACt                             >  6:89364/1‑149 (MQ=255)                                                                                                                           aaGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATAcc                       <  3:270297/149‑1 (MQ=255)                                                                                                                                      cAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGcc                                            <  2:24859/117‑1 (MQ=255)                                                                                                                                      cAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGcc                                            >  1:24859/1‑117 (MQ=255)                                                                                                                                      cAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCa            <  3:18236/149‑1 (MQ=255)                                                                                                                                        gCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAgg          >  4:24323/1‑149 (MQ=255)                                                                                                                                         cGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAg                                               >  8:228360/1‑111 (MQ=255)                                                                                                                                         cGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAg                                               <  7:228360/111‑1 (MQ=255)                                                                                                                                         cGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAggg         >  5:163076/1‑149 (MQ=255)                                                                                                                                           aCGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAGGGTg       <  6:318627/149‑1 (MQ=255)                                                                                                                                           aCGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAGGGTg       <  4:262376/149‑1 (MQ=255)                                                                                                                                            cgccgcCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAGGGTgg      <  4:57356/149‑1 (MQ=255)                                                                                                                                             gccgccGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAGGGTGGc     >  1:123180/1‑149 (MQ=255)                                                                                                                                              tcgccgTGCTGGCCAGTACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAGGGTGGc     >  1:136393/2‑148 (MQ=255)                                                                                                                                                gacgTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGTCAAAGACGTTCTCCTGTATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAGGGTGGCGAt  <  5:155411/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CGGCCAGGCCGGCGCAGATGGTGGCATAGCGCGGCGTTTGCGAGCGCGGGCAGACACTGGCGAACTTCTGCGGTACCGCACCGTAATAACCCATGGCCAACAGGGTGCGCGCCGGCGCCACGAAGGTGGACTGCAGCGACGCCGCCGTACTGGCCAGCACGGCGATGGACATCAGGATCGCCAGCGGCCCCATGACCGGGCCGGCCAGGTGGGCAAAGACGTTCTCCTGGATGCGCGGGTTGTTCAGGCCCAGACCCGTTTCACTGATACCGGCAAACTGCAGGGTGGCGAT  >  NC_002947/6044092‑6044383

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