Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,867,733:1	+C	8.6%	intergenic (‑28/+24)	thyA ← / ← lgt	thymidylate synthase/prolipoprotein diacylglyceryl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,867,733	1	.	C	8.6%	216.9 / 14.0	58	intergenic (‑28/+24)	thyA/lgt	thymidylate synthase/prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base . (31/22);  new base C (2/3);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

AGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAA  >  NC_002947/5867593‑5867874                                                                                                                                              |                                                                                                                                              aGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGccc                                                                                                                                       >  6:198636/1‑149 (MQ=255)      gAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGcc                                                                                                                                  >  5:91640/1‑149 (MQ=255)          cGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGATC‑CCCCTGCCGttt                                                                                                                              <  1:34738/149‑3 (MQ=255)            cagcaTCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGGCGAT‑‑CCCCTGCCGtttttg                                                                                                                           <  6:27453/149‑6 (MQ=255)                    cgcCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTa                                                                                                                         <  3:279404/144‑1 (MQ=255)                     gcCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAgg                                                                                                                   >  3:23128/1‑149 (MQ=255)                        cGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCgg                                                                                                                <  6:91640/149‑1 (MQ=255)                        cGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCgg                                                                                                                >  4:88452/1‑149 (MQ=255)                          gCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGt                                                                                                              >  6:137639/1‑149 (MQ=255)                                   tGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTa                                                                                                     <  5:189878/149‑1 (MQ=255)                                     gTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAg                                                                                                    <  3:126355/148‑1 (MQ=255)                                             gCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGa                                                                                            <  5:198636/148‑1 (MQ=255)                                                     ggTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGcc                                                                                   >  6:242171/1‑149 (MQ=255)                                                       tttCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCcca                                                                                 >  3:131756/1‑149 (MQ=255)                                                           ccTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCccac                                                                                >  2:97926/1‑146 (MQ=255)                                                           ccTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCccac                                                                                >  1:31991/1‑146 (MQ=255)                                                           ccTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAg                                                                             <  3:88452/149‑1 (MQ=255)                                                            cTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCccac                                                                                >  3:168673/1‑145 (MQ=255)                                                             tGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCccac                                                                                >  7:240915/1‑144 (MQ=255)                                                               cAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGAt                                                                           <  1:97926/147‑1 (MQ=255)                                                                             cAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCAcc                                                            >  8:13659/1‑148 (MQ=255)                                                                              aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGt                                                                                                      >  2:166905/1‑105 (MQ=255)                                                                              aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGt                                                                                                      <  1:166905/105‑1 (MQ=255)                                                                                      cgcgGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcatc                                                  >  6:126978/1‑149 (MQ=255)                                                                                        cgGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCgg                                                <  2:31991/149‑1 (MQ=255)                                                                                             cAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGc                                           <  4:23128/149‑1 (MQ=255)                                                                                                gtacAGATACTGTTTAATGGATTACCTGTAGCAAAAACGGCAGGTG‑GATAGCACCTTCCGTGGGATTAATAATGAGGTTTTCGCCGTAGGCTTGCGATTATAATCCCAAAAGATCAAGACTAGGACACAGAATATCATCGGAATGca                                          >  4:120433/4‑147 (MQ=255)                                                                                                gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGcac                                         >  6:218455/1‑148 (MQ=255)                                                                                                gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGcac                                         <  5:218455/148‑1 (MQ=255)                                                                                                gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGca                                          >  5:159109/1‑147 (MQ=255)                                                                                                gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGca                                          >  4:115829/1‑147 (MQ=255)                                                                                                gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGACCACAAGATCAATACGAGGCCACCGACTATCATCGGAATTcaca                                        >  3:13028/1‑149 (MQ=255)                                                                                                gTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCATGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTTCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATgaa                                          >  1:202716/1‑145 (MQ=255)                                                                                                  ccAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTgg                                                                                                                              <  2:93643/61‑1 (MQ=255)                                                                                                  ccAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTgg                                                                                                                              >  1:93643/1‑61 (MQ=255)                                                                                                  ccAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAg                                       >  8:89502/1‑148 (MQ=255)                                                                                                   cAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGc                                      >  4:261389/1‑148 (MQ=255)                                                                                                       tACTGTTTCATGGCTTACCTGTTGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGAGGGTTTCGACGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTg                                 >  7:45865/1‑149 (MQ=255)                                                                                                         cTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGc                               >  4:292898/1‑149 (MQ=255)                                                                                                          tGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAAc                                   <  1:82950/144‑1 (MQ=255)                                                                                                          tGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAAc                                   >  2:82950/1‑144 (MQ=255)                                                                                                            tttCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCtt                            <  8:240915/149‑1 (MQ=255)                                                                                                            tttCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCt                             <  5:137639/148‑1 (MQ=255)                                                                                                               cATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcatc                                                  <  6:154351/124‑1 (MQ=255)                                                                                                               cATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcatc                                                  >  5:154351/1‑124 (MQ=255)                                                                                                               cATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtga                          <  3:261389/148‑1 (MQ=255)                                                                                                                       aCCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGc                                           <  7:22559/123‑1 (MQ=255)                                                                                                                       aCCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGc                                           >  8:22559/1‑123 (MQ=255)                                                                                                                             aGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAc           >  5:6043/1‑149 (MQ=255)                                                                                                                              gCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCACGCGGGATTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAcg          <  3:207244/149‑1 (MQ=255)                                                                                                                                aaaaaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAcgc         >  3:170651/1‑148 (MQ=255)                                                                                                                                 aaaaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGCCGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa       >  3:80384/1‑149 (MQ=255)                                                                                                                                  ccaCGGCAGGGGCGATC‑CCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagt                        >  3:234944/3‑131 (MQ=255)                                                                                                                                  ccaCGGCAGGGGCGATC‑CCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagt                        <  4:234944/129‑1 (MQ=255)                                                                                                                                   aaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATg     <  3:120433/149‑1 (MQ=255)                                                                                                                                   aaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATg     <  3:115829/149‑1 (MQ=255)                                                                                                                                     cGGCAGGGGCGATC‑CCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGt    >  6:210250/1‑148 (MQ=255)                                                                                                                                      ggCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTaa  <  6:159109/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                              AGGTCGAAACGCAGCATCGCGCCCGGCAGGCTGATGGTGCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAA  >  NC_002947/5867593‑5867874

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