Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,153,755	(G)5→6	coding (193/660 nt)	yghB →	DedA family protein YghB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,153,750	1	.	G	97.4%	130.6 / ‑3.2	38	coding (188/660 nt)	yghB	DedA family protein YghB
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (19/18);  total (19/19)
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.33e-01

CGCGCTGGCGGATCCTCATATTGTTAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑A‑GGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCA  >  NC_000913/3153613‑3153874                                                                                                                                            |                                                                                                                             cgcgCTGGCGGATCCTCATATTGTTAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑A‑gg                                                                                                                            <  1:24654/140‑1 (MQ=255) cgcgCTGGCGGATCCTCATATTGTTAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑A‑gg                                                                                                                            <  1:801041/140‑1 (MQ=255)    gcTGGCGGATCCTCATATTGTTAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑A‑ggggg                                                                                                                         <  2:68917/140‑1 (MQ=255)        gCGGATCCTCATATTGTTAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑A‑ggggggtta                                                                                                                     <  1:97421/140‑5 (MQ=255)                tCATATTGTTAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGAtttt                                                                                                              >  1:345388/1‑140 (MQ=255)                    atTGTTAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGc                                                                                                          <  1:91015/140‑1 (MQ=255)                        ttAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTAc                                                                                                      >  1:122240/1‑140 (MQ=255)                        ttAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTAc                                                                                                      >  2:358802/1‑140 (MQ=255)                        ttAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTAc                                                                                                      >  1:22728/1‑140 (MQ=255)                        ttAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTAc                                                                                                      <  1:605158/140‑1 (MQ=255)                        ttAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTAc                                                                                                      >  1:631640/1‑140 (MQ=255)                                        tGTCATGTTTGCCACGCTGTTGTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGAcc                                                                                      >  2:819062/1‑140 (MQ=255)                                          tCATTTTTGCCCCGTTGTTTTTAGAAACCGGCCTGCGGCCCGCCTCATTTTTGCCAGGCGACAGTTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGAccgc                                                                                    <  2:102291/140‑1 (MQ=255)                                                ttGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCgcagc                                                                              >  2:128434/1‑140 (MQ=255)                                                ttGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCgcagc                                                                              >  1:8957/1‑140 (MQ=255)                                                           tttttAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGgct                                                                   <  2:345388/140‑1 (MQ=255)                                                           tttttAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGgct                                                                   >  1:748098/1‑140 (MQ=255)                                                                aGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTgg                                                              <  1:358802/140‑1 (MQ=255)                                                                  aaaaCGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCt                                                            <  1:855920/140‑1 (MQ=255)                                                                  aaaaCGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCt                                                            <  2:631640/140‑1 (MQ=255)                                                                  aaaaCGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCt                                                            <  2:724306/140‑1 (MQ=255)                                                                  aaaaCGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCt                                                            <  2:748098/140‑1 (MQ=255)                                                                  aaaaCGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCt                                                            <  1:161328/140‑1 (MQ=255)                                                                     aCGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGt                                                                          >  1:750802/1‑123 (MQ=255)                                                                     aCGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGt                                                                          <  2:750802/123‑1 (MQ=255)                                                                     aCGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAg                                                         >  1:447123/1‑140 (MQ=255)                                                                     aCGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAg                                                         >  2:457462/1‑140 (MQ=255)                                                                            gctgcCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATAtt                                                  <  1:719577/140‑1 (MQ=255)                                                                                   cccTCTTTTTTCCCAGGCGACAGTTTGTTGATACTGGCAGGCGCATTGATTGCCCGG‑GGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGgc                                           <  2:865490/139‑1 (MQ=255)                                                                                    ccTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGgcg                                          <  1:353896/140‑1 (MQ=255)                                                                                       cATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAg                                               <  2:554986/132‑1 (MQ=255)                                                                                       cATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAg                                               >  1:554986/1‑132 (MQ=255)                                                                                       cATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTg                                       >  1:865490/1‑140 (MQ=255)                                                                                              gCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAggg                                <  2:22728/140‑1 (MQ=255)                                                                                              gCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAggg                                <  1:310828/140‑1 (MQ=255)                                                                                                        aGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCaaaa                      >  2:164465/1‑140 (MQ=255)                                                                                                             gttgATACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTg                 >  1:194709/1‑140 (MQ=255)                                                                                                                  tACTGGCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAgg            <  2:194709/140‑1 (MQ=255)                                                                                                                      ggCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAggctgg        <  1:128434/140‑1 (MQ=255)                                                                                                                       gCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAggctggc       >  2:614537/1‑140 (MQ=255)                                                                                                                       gCAGGCGCATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAggctggc       >  1:757340/1‑140 (MQ=255)                                                                                                                            cgcATTGATTGCCC‑AGGGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGca  >  2:328552/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                             CGCGCTGGCGGATCCTCATATTGTTAGCGTTGTTTACTTTGTCATGTTTGCCACGCTGTTTTTAGAAAACGGCCTGCTGCCCGCCTCATTTTTGCCAGGCGACAGCTTGTTGATACTGGCAGGCGCATTGATTGCCC‑A‑GGGGGTTATGGATTTTCTGCCTACGATTGCGATTCTGACCGCCGCAGCAAGTCTGGGCTGCTGGCTAAGTTATATTCAGGGGCGCTGGTTAGGGAATACCAAAACGGTGAAAGGCTGGCTGGCA  >  NC_000913/3153613‑3153874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑