Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,286,328	(T)7→6	intergenic (+41/+62)	proL → / ← yejO	tRNA‑Pro/adhesin‑like autotransporter YejO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,286,322	0	T	.	78.3%	32.2 / ‑0.4	23	intergenic (+35/+68)	proL/yejO	tRNA‑Pro/adhesin‑like autotransporter YejO
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base . (6/12);  total (8/15)
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.66e-01

AGACGAGAAGCG‑TTTTATCGCTAACTGATTAATTAT‑AAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACC‑AAAAATCCCAAG‑AAAAAACCAACCCTTACGGTTGG‑TTTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑AAA‑CGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCAGCGTCGTACCGTTTTGCAGGTC  >  NC_000913/2286109‑2286647                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                                                               aGACGAGAAGCGTTTTTATCGCTAACTGATTAATTAT‑AAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACC‑AAAAATCCCAAGAAAAAAACCAACCCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACC‑TGACACATCACAGCC‑TG‑TTTATTTTCTGTTATc                                                                                                                                                                                                                                                                       >  1:1260254/1‑283 (MQ=255)        aaGCG‑TTTTATCGCTAACTGATTAATTAT‑AAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACC‑AAAAATCCCAAGAAAAAAACCAACCCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAg                                                                                                                                                                                                                                                                     >  1:1260255/1‑278 (MQ=255)                            gATTAATTATAAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAAATCCCAAG‑‑‑AAAACCAACCCTTACGGTTGG‑TTTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGa                                                                                                                                                                                                                               >  1:1260256/1‑297 (MQ=255)                                                                                                      aGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACC‑AAAAATCCCAAG‑‑‑AAAACCAACCCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCATA‑AAACAGACCGTG‑CACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑AAA‑c                                                                                                                                   >  1:1260264/1‑307 (MQ=255)                                                                                                                            gTAGCGCACCGTCATGGGGTGTC‑GGGGTCGGAGGTTCAAATCCTCTCGTGCCGACC‑‑AAAATCCCAAG‑AAAAAACCAACCCTTACGGTTG‑‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑a                                                                                                                                       <  1:1260263/282‑1 (MQ=255)                                                                                                                                   aCCGTCATGGGGTGTC‑GGGGTCGGAGGTTCAAATCCTCTCGTGCCGACC‑‑AAAATCCCAAG‑‑AAAAACCAACCCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TGTTTTATTTTCTGtt                                                                                                                                                                                                                                                                          <  1:1260265/149‑1 (MQ=255)                                                                                                                                                       gTCGGAGGTTCAAATCCTCTCGTGCC‑ACC‑AAAAATCCCAAG‑AAAAAACCAACCCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTTAACGTCGGGGGTAAACGATGACCTTATACCAGCCTGAT‑AAA‑CGCCACt                                                                                                                             <  1:1260266/267‑1 (MQ=255)                                                                                                                                                                tCAAATCCTCTCGTGCCGACCAAAAAATCCCAAG‑AAAAAACCAACCCTTACGGTTGG‑‑‑TTTTTATATCTGC‑AATTAATTCGATAAAACAGACCGTGACACATCACAGCC‑TG‑TTTA‑TTTCTG‑TATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGAt                                                                                                                                                                                                                              <  1:413201/165‑1 (MQ=255)                                                                                                                                                                tCAAATCCTCTCGTGCCGACC‑AAAAATCCCAAG‑AAAAAACCAACCCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑AAA‑CGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAAC‑TTCACCt                                                                           >  1:1260271/1‑307 (MQ=255)                                                                                                                                                                  aaaTCCTCTCGTGCCGACC‑AAAAATCCCAAG‑‑AAAAACCAACCCTTACGGTTGG‑‑‑TTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAG‑ACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑‑AATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑a                                                                                                                                       <  1:1260267/243‑1 (MQ=255)                                                                                                                                                                        tctcGTGCCGACC‑‑AAAATCCCAAG‑‑AAAAACCAACCCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TGTTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑cccc                                                                                                                                                                                                                     <  1:1260268/165‑1 (MQ=255)                                                                                                                                                                           cGTGCCGACC‑‑AAAATCCCAAG‑AAAAAACCAACCCTTACGGTTG‑‑‑‑TTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCCTGCGCCATTGCCATAGCTGACTGAC‑aatg                                                                                                                                                                                       <  1:1260269/190‑3 (MQ=255)                                                                                                                                                                              gCCGACC‑AAAAATCCCAAG‑‑AAAAACCAACCCTTACGGTTG‑‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TGTTTTATTTTCTG‑TATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTTAACGTC‑GGGGTAAACGATGACCTtat                                                                                                                                                    <  1:413202/221‑1 (MQ=255)                                                                                                                                                                               ccGACC‑AAAAATCCCAAG‑AAAAAACCAACCCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTa                                                                                                                                                                                                                                          >  1:1260274/1‑139 (MQ=255)                                                                                                                                                                               ccGACC‑AAAAATCCCAAG‑‑‑AAAACCAACCCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑AAA‑CGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCAt                                                                                       >  1:1260273/1‑279 (MQ=255)                                                                                                                                                                                  aCC‑AAAAATCCCAAG‑AAAAAACCAACCCTTACGGTTG‑‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑‑AATGACCGCTT‑ACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑AAA‑CGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGc                                                 <  1:1260270/313‑1 (MQ=255)                                                                                                                                                                                                         ccaaccCTTACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGct                                                                                                                                                                                                   <  1:1260272/153‑1 (MQ=255)                                                                                                                                                                                                            accCTTACGGTT‑G‑TTTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCCTTGTTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTT‑ACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑AAACCGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCAGCGTCGCA‑CGTTTTGCAGGTc  <  1:413203/340‑1 (MQ=255)                                                                                                                                                                                                             cccTTACGGTTG‑‑‑‑TTTTTATATCTGCAAATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TGTTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑‑AATGACCGCTT‑ACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGATAAAA‑CGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGcc                                 <  1:1260275/306‑1 (MQ=255)                                                                                                                                                                                                              ccTTACGGTTGG‑TTTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCCTTGTTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑AAA‑CGCCACTGGTTCCAGcc                                                                                                                   <  1:413208/227‑1 (MQ=255)                                                                                                                                                                                                              ccTTACGGTTGG‑TTTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TGTTTTATTTTCTG‑TATCAGAACGTCCAGACCACACCCGCCTg                                                                                                                                                                                                                                              <  1:413205/107‑1 (MQ=255)                                                                                                                                                                                                                ttACGGTTGG‑‑TTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TGTTTTATTTTCTG‑TATCAG‑ACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGACAAAATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGATAAAA‑CGCCACTGGTTCCAGCCACATCATTATTtaa                                                                                                     <  1:413204/237‑1 (MQ=255)                                                                                                                                                                                                                 tACGGTTGG‑TTTTTTTATATCTGC‑AATT‑ATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTA‑TTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑‑CCCTGCGCC‑TTGCCATAGCTGACTGAC‑‑AATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑AAA‑CGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATc                                                                 >  1:413213/1‑267 (MQ=255)                                                                                                                                                                                                                    ggttggttTTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TGTTTTATTTTCTGTTATCAGAACGTCCAGACcaca                                                                                                                                                                                                                                                      <  1:413206/95‑1 (MQ=255)                                                                                                                                                                                                                        tg‑ttTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCCTTGTTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGa                                                                                                                                                                                                                               <  1:413209/113‑1 (MQ=255)                                                                                                                                                                                                                          ttttttttATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACA‑CACAGCC‑TG‑TTTATTTTCTGTTATCAGa                                                                                                                                                                                                                                                                    <  1:413211/72‑1 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                                                               AGACGAGAAGCG‑TTTTATCGCTAACTGATTAATTAT‑AAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACC‑AAAAATCCCAAG‑AAAAAACCAACCCTTACGGTTGG‑TTTTTTTATATCTGC‑AATTAATTCGAT‑AAACAGACCGTGACACATCACAGCC‑TG‑TTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTA‑CCCCTGCGCCATTGCCATAGCTGACTGAC‑AAATGACCGCTTAACGTC‑GGGGTAAACGATGACCTTATACCAGCCTGAT‑AAA‑CGCCACTGGTTCCAGCCACATCATTATTAAATTTGCCATCGTCATTAACTTTCACCTGGTTAGAATCGGCGTATTCCTGACGCACGGCCGCTTTCAGCCAGGGTTCCAGCGTCGTACCGTTTTGCAGGTC  >  NC_000913/2286109‑2286647

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