Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,256,416	A→T	87.8%	intergenic (‑464/+305)	ycgV ← / ← ychF	putative adhesin/catalase inhibitor protein; ATPase, K+‑dependent, ribosome‑associated

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,256,416	0	A	T	87.8%	122.5 / 8.5	41	intergenic (‑464/+305)	ycgV/ychF	putative adhesin/catalase inhibitor protein; ATPase, K+‑dependent, ribosome‑associated
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base T (15/21);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CTAATTAGTTATGGGAATTTCGCTTGATGTATCTGCATAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTACTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGATATCCGGCATCTGAACACTTATGTCTTTAATTAT  >  NC_000913/1256271‑1256564                                                                                                                                                  |                                                                                                                                                     cTAATTAGTTATGGGAATTTCGCTTGATGTATCTGCATAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTg                                                                                                                                                   <  2:473172/149‑1 (MQ=255)   aaTTAGTTATGGGAATTTCGCTTGATGTATCTGCATAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGt                                                                                                                                                 >  1:624570/1‑149 (MQ=255)                        ttGATGTATCTGCATAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTACTGGTTATCATTTTCATCTTGCCAGc                                                                                                                            >  1:554123/1‑149 (MQ=255)                             gTATCTGCATAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATAtt                                                                                                                       <  3:44387/149‑1 (MQ=255)                                  tGCATAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGc                                                                                                                  >  3:615186/1‑149 (MQ=255)                                     aTAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTg                                                                                                               >  2:383926/1‑149 (MQ=255)                                      tAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGa                                                                                                              <  3:295027/149‑1 (MQ=255)                                      tAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTACTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGa                                                                                                              <  2:554123/149‑1 (MQ=255)                                        aCAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATc                                                                                                            <  2:180465/149‑1 (MQ=255)                                        aCAAAAACAATATGCGCCACTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATc                                                                                                            <  2:10282/149‑1 (MQ=255)                                         cAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCa                                                                                                           >  4:162374/1‑149 (MQ=255)                                           aaaaCAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAAt                                                                                                         <  4:111008/149‑1 (MQ=255)                                            aaaCAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAAtt                                                                                                        >  2:167415/1‑149 (MQ=255)                                            aaaCAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAAtt                                                                                                        <  1:531748/149‑1 (MQ=255)                                                 atatGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCCATTTTGAt                                                                                                   >  2:536761/1‑149 (MQ=255)                                                       gcCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCt                                                                                              <  4:561705/148‑1 (MQ=255)                                                       gcCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCt                                                                                              >  3:561705/1‑148 (MQ=255)                                                       gcCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATGGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCtg                                                                                             >  3:516741/1‑149 (MQ=255)                                                            ttttGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAg                                                                                        <  3:162374/149‑1 (MQ=255)                                                             tttGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGGAgg                                                                                       <  1:536761/149‑1 (MQ=255)                                                                gCATTACTAATAAGAAAAAGCAAAAAACACCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTTTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAggggg                                                                                    <  2:346418/149‑1 (MQ=255)                                                                         ataaGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATgtgt                                                                           >  2:595927/1‑149 (MQ=255)                                                                           aaGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAg                                                                         <  4:325238/149‑1 (MQ=255)                                                                             gAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTg                                                                       >  3:42173/1‑149 (MQ=255)                                                                             gAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTg                                                                       <  2:624570/149‑1 (MQ=255)                                                                             gAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTg                                                                       >  1:160597/1‑149 (MQ=255)                                                                                aaaGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATa                                                                    >  2:425510/1‑149 (MQ=255)                                                                                                 tatTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTtata                                                   <  1:425510/149‑1 (MQ=255)                                                                                                   ttCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTACTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATc                                                 >  4:343404/1‑149 (MQ=255)                                                                                                       tttttGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTg                                             >  4:90664/1‑149 (MQ=255)                                                                                                             gTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATg                                       <  1:383926/149‑1 (MQ=255)                                                                                                             gTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATg                                       <  4:51637/149‑1 (MQ=255)                                                                                                                 aTGGGTTTTCTGGCCTATCGTTGTTATTGTGCTACTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGat                                   >  4:10578/1‑149 (MQ=255)                                                                                                                          cTGGCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGATATCCGGCAt                          <  1:595927/149‑1 (MQ=255)                                                                                                                            ggCCTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGATATCCGGCATCt                        <  2:400022/149‑1 (MQ=255)                                                                                                                               cTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGATATCCGGCATCTGaa                     <  3:377599/149‑1 (MQ=255)                                                                                                                               cTATCGTTGTTATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGATATCCGGCATCTGaa                     <  1:268384/149‑1 (MQ=255)                                                                                                                                    gttgttATTGTGCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGATATCCGGCATCTGAACACtt                >  1:499168/1‑149 (MQ=255)                                                                                                                                             gtgCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTa                                                          <  2:301864/98‑1 (MQ=255)                                                                                                                                             gtgCTTCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGAGCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTa                                                          >  1:301864/1‑98 (MQ=255)                                                                                                                                              tgCTACTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGATATCCGGCATCTGAACACTTATGTCTTTaa      <  3:10578/149‑1 (MQ=255)                                                                                                                                                  aCTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGATATCCGGCATCTGAACACTTATGTCTTTAattat  <  3:343404/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CTAATTAGTTATGGGAATTTCGCTTGATGTATCTGCATAACAAAAACAATATGCGCCATTTTTGCATTACTAATAAGAAAAAGCAAAAAACATCGATATTCATTTTTGGTCAATGGGTTTTCTGGCCTATCGTTGTTATTGTGCTACTGGTTATCATTTTCATCTTGCCAGCATATTGGAGCGTGATCAATTTTGATCAGCTGTGAAGGGGGCTGATGTGTAGTGATACAAAAACTTACGTTATATCACTGCTAATGAGATATCCGGCATCTGAACACTTATGTCTTTAATTAT  >  NC_000913/1256271‑1256564

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