Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,991,839	G→A	W229* (TGG→TGA)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	3,991,839	0	G	A	97.8%	128.7 / NA	46	W229* (TGG→TGA)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/23);  minor base T (0/1);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACC  >  NC_000913_3_bme_pgi/3991705‑3991978                                                                                                                                       |                                                                                                                                            gAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGAGTGCTGACGCCAAATGAATGACTGGATCTCGGTgg                                                                                                                               >  1:127184/1‑149 (MQ=255)    ttttATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCtt                                                                                                                            >  1:12410/1‑149 (MQ=255)    ttttATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCtt                                                                                                                            >  2:677118/1‑149 (MQ=255)           gTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCGAATGAATGACTGGATCTCGGTGGCTTAAGCTCg                                                                                                                     >  1:232419/1‑149 (MQ=255)           gTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCg                                                                                                                     <  1:464240/149‑1 (MQ=255)               cgccgTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAGGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCttt                                                                                                                 >  1:613891/1‑149 (MQ=255)                 ccgTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTct                                                                                                               <  2:579112/149‑1 (MQ=255)                 ccgTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTct                                                                                                               >  1:519337/1‑149 (MQ=255)                         ctcGCCTGGAAGCGTATTGTTGGGAATGTGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCCTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAaga                                                                                                       <  2:632023/149‑1 (MQ=255)                         ctcGCCGGTAAGCGTATTCTGGGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAaga                                                                                                       <  1:77040/149‑1 (MQ=255)                            gCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTa                                                                                                    <  2:426161/149‑1 (MQ=255)                                  aaGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTgg                                                                                              <  1:111780/149‑1 (MQ=255)                                    gCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTg                                                                                            <  1:673540/149‑1 (MQ=255)                                       tatgttgtgGAATATGGTGCCGTTCGTCGAAGAAGAGCGTTACGACGACTATGTGATGTCGCTTTACGCGGAGGGCGTGCTGACGACAAATGAATGTCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCa                                                                                         <  2:794435/144‑1 (MQ=255)                                            tgtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCtt                                                                                    >  2:220556/1‑149 (MQ=255)                                            tgtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCtt                                                                                    >  1:270438/1‑149 (MQ=255)                                            tgtgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGATTTCTGCGGAAGGGTCATTTGGTGCCAGTCtt                                                                                    >  1:513310/1‑149 (MQ=255)                                              tgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTACGCTCGCTTTCTGCTGAAGAGTACTGTGGTGCCAGCCTTTg                                                                                  >  2:264583/1‑149 (MQ=255)                                              tgGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTACGCCCGCTTTCTGCCGAAGAGTACTTTCGTGCCCGCCttca                                                                                  >  2:607367/1‑147 (MQ=255)                                                gAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTACGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGc                                                                                >  2:72452/1‑149 (MQ=255)                                                 aaTATGGTGCCGTGCGACGAAGAAGAGCATTACGAAGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTAAATAGTACTTTGGTGCCAGCCTTTGGCa                                                                               >  2:543689/1‑149 (MQ=255)                                                   tatGGTTCCGTACGACGAAGCAGAGCATCACGACGACTCTGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCCAATGAACGACTGGATCTCTGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGc                                                                             <  2:513310/149‑1 (MQ=255)                                                      ggTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGctct                                                                          >  2:484101/1‑149 (MQ=255)                                                       gTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGctctc                                                                         >  1:777874/1‑148 (MQ=255)                                                          cccagCGACTAAGAAGAGCATTACGACGACTATGTCATGTCGCTTTACGCGCAGGGCGTGCCGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACaa                                                                      <  1:570703/145‑1 (MQ=255)                                                          ccGTGCGACGAAGAAGAGCTGTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACaa                                                                      <  1:265230/149‑1 (MQ=255)                                                                 acgaGGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGATCGCTTTATTCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATc                                                               <  2:613891/149‑1 (MQ=255)                                                                   gaagaaGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTACCAGCCTTTGGCAGCTCTACAAGAGTATCGa                                                             >  1:314025/1‑149 (MQ=255)                                                                              ttACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATCCa                                                  >  2:560349/1‑149 (MQ=255)                                                                                acgacgACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACaaa                                                <  1:458468/149‑1 (MQ=255)                                                                                     gacTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGt                                           >  2:189413/1‑149 (MQ=255)                                                                                        tATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTAGGCAGCTCTACAAGAGTATCGATTCTCCATACAAAGCGGTACt                                        >  2:715281/1‑149 (MQ=255)                                                                                                gatgCTTTACGCGCTGGGCGTGTTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCGGAAGAGTACTTTGGGGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAacac                                <  2:222380/146‑1 (MQ=255)                                                                                                gaCGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAacac                                >  2:763036/1‑149 (MQ=255)                                                                                                  ctcTTTACGCGCAGGGCGTGCTGATGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActg                              <  1:524710/147‑1 (MQ=255)                                                                                                  cGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActg                              >  1:346582/1‑149 (MQ=255)                                                                                                  cGCATTACGCGCAGGGCGGGCTTTCGCCAAATGAATGACGGGATCTCGGGGGATTAAGCTCGCTTTATGCTGAAGAGTACCTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAACCGGTACTGAAAACActg                              <  2:261632/149‑1 (MQ=255)                                                                                                           cgcAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGaa                     <  2:12410/149‑1 (MQ=255)                                                                                                            gcAGGGCGGGGTGGCGCCAAATGAATGACTGGATCTCCGTGGCTGGACCGCGCTTTCTGTTGAAGAGTACTTTGGGGCCAGCCTGTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAg                    <  2:95644/149‑1 (MQ=255)                                                                                                             cAGGGCGTGCTGACGCCAAATGAATGACTGGATCTCGGTGGGTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATCCAAAGCGGTACTGAAAACACTGCTGCTGGAAGc                   <  2:646871/149‑1 (MQ=255)                                                                                                                     gCTGACGCCAAATGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCt           <  1:189413/149‑1 (MQ=255)                                                                                                                       tGACGCCAAATGAATGACTGGATCTCGTTGGCGTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTgg         <  1:608879/149‑1 (MQ=255)                                                                                                                       tGACGCCAAATGAATGACTGGATCTCGGTGGCTTACGCGCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTgg         <  2:520955/149‑1 (MQ=255)                                                                                                                              aaaTGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATa                                                    >  2:702585/1‑99 (MQ=255)                                                                                                                              aaaTGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATa                                                    <  1:702585/99‑1 (MQ=255)                                                                                                                              aaaTGAATGACTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAcc  <  2:758368/149‑1 (MQ=255)                                                                                                                                       |                                                                                                                                            GAATTTTATCGTACCGCCGTGCGTCTCGCCGGTAAGCGTATTCTGTGGAATATGGTGCCGTGCGACGAAGAAGAGCATTACGACGACTATGTGATGACGCTTTACGCGCAGGGCGTGCTGACGCCAAATGAATGGCTGGATCTCGGTGGCTTAAGCTCGCTTTCTGCTGAAGAGTACTTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACC  >  NC_000913_3_bme_pgi/3991705‑3991978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑