Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	985,846	(A)8→7	intergenic (‑137/+48)	aspC ← / ← ompF	aspartate aminotransferase/outer membrane porin F

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	985,839	0	A	.	77.8%	27.4 / ‑2.3	18	intergenic (‑130/+55)	aspC/ompF	aspartate aminotransferase/outer membrane porin F
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base . (6/8);  total (6/12)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TTTATAGACACCAATCCCGAGGTTAATTTTGCCGGGACG‑TTCATCGGCACGAAACAGATCGGCCAGGCCCAGAATC‑GGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑AAAAAAAACAGGACC‑A‑AAGTCCTG‑TTTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCCACAGCAACGGT‑GTCGTCTGAACCTACGCCCAGTTTGT‑TGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCG  >  NC_000913/985599‑986122                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                      tttATAGACACCAATCCCGAGGTTAATTTTGCCGGGACG‑TTCATCGGCACGAAACAGATCGGCCAGGCCCAGAATC‑GGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAACATATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑‑‑AAAAAACAGGACC‑A‑AAGTCCTG‑TTTTTTCGGCA‑TTT‑AAC‑AAAGAGGtgtg                                                                                                                                                                                                                                              >  1:869803/1‑291 (MQ=255)            cAATCCCGAGGTTAATTTTGCCGGGACG‑TTCATCGGCACGAAACAGATCGGCCAGGCCCAGAATC‑GGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑‑‑‑AAAAACAGGTCC‑A‑AAGTCCTG‑TTTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGCTTATTAG‑AACTGGTAAACGATACCCACAGCAACGGT‑GTCGTCTGAACCTACg                                                                                                                                                                                       >  1:869807/1‑331 (MQ=255)                     ggTTAATTTTGCCGGGACG‑TTCATCGGCACGAAACAGATCGGCCAGGCCCAGAATCGGGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑aaaaaaaa                                                                                                                                                                                                                                                                                                <  1:869801/229‑1 (MQ=255)                        tAATTTTGCCGGGACGTTTCATCGGCACGAAACAGATCGGCCAGGCCCAGAATC‑GGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTGTC‑TCTGG‑‑AAGGTAT‑aaaaa                                                                                                                                                                                                                                                                                                   >  1:869809/1‑224 (MQ=255)                          aTTTTGCCGGGACG‑TTCATCGGCACGAAACAGATCGGCCAGGCCCAGAATC‑GGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGGA‑AAGGTAT‑‑AAAAAAACAGGACC‑A‑AAGTCCTG‑‑TTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCCACAGCAAc                                                                                                                                                                                                           <  1:869804/299‑1 (MQ=255)                                               ggCACGAAACAGATCGGCCAGGCCCAGAATC‑GGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑‑‑‑AAAAACAGGACC‑T‑AAGTCCTG‑TTTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCCACAGCAACGGTCGTCGTCTGAACCTACGCCCAGTTTGTCTGTCAGAATc                                                                                                                                                                  >  1:869813/1‑318 (MQ=255)                                                                    gCCCAGAATC‑GGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑aaaaa                                                                                                                                                                                                                                                                                                   <  1:869812/179‑1 (MQ=255)                                                                         gAATC‑GGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATG‑AAACG‑CAGGCTG‑TTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑‑‑‑AAAAACAGGACC‑‑‑AAGTCCTG‑‑TTTTTCGGCA‑TTT‑AAC‑‑AAGAGGTGTGCTTATTAG‑AACTGGTAAACGATACCCACAGc                                                                                                                                                                                                              >  1:121027/1‑244 (MQ=255)                                                                                                   aaTGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGGA‑AAGGTAT‑‑AAAAAAACAGGACC‑A‑AAGTCCTG‑‑‑TTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCCACAGCAACGGT‑GTCGTCTGAACCTACGCCCAGTTTGT‑TGTCAGAATCGATCTGGTTGatgatg                                                                                                                                                  <  1:869814/282‑1 (MQ=255)                                                                                                                                        aGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑‑‑AAAAAACAGGACC‑a                                                                                                                                                                                                                                                                                       >  1:869819/1‑121 (MQ=255)                                                                                                                                          aaGGGAAGTCCGGCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACGTCAGGCTG‑TTTTGCAAGACGTGAGATTG‑C‑TCTGGA‑AAGGTAT‑‑AAAAAAACAGGACC‑A‑AAGTCCTG‑‑‑TTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGc                                                                                                                                                                                                                                             <  1:869817/157‑1 (MQ=255)                                                                                                                                                       cTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTG‑TTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑‑AAAAAAACAGGACCAA‑AAGTCCTG‑TTTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCca                                                                                                                                                                                                                  <  1:869818/170‑1 (MQ=255)                                                                                                                                                                       agaTTTACAAAATTCCAAC‑TATTAACTGATGAAAACGTCAGGCTG‑TTTTGCAAGACGTGAGATTG‑C‑TCTGGA‑AAGGTAT‑‑AAAAAAACAGGACC‑A‑AAGTCCTG‑‑TTTTTCGGCA‑TTT‑AAc                                                                                                                                                                                                                                                          <  1:869823/118‑1 (MQ=255)                                                                                                                                                                                  aTTCCAAC‑TATTAACTGATGAAAACG‑CAGGCTG‑TTTTGCAAGACGTGAGATTG‑C‑TCTGGA‑AAGGTAT‑‑AAAAAAACAGGACC‑A‑AAGTCCTGTTTTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TAtt                                                                                                                                                                                                                                        <  1:121033/124‑1 (MQ=255)                                                                                                                                                                                                                     tttttGCAAGACGTGAGATTG‑CATCTGGA‑AAGGTAT‑‑AAAAAAACAGGACC‑A‑AAGTCCTG‑‑TTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCCACAGCAACGGT‑GTCGTCTGAACCTACGCCCAg                                                                                                                                                                                  <  1:869826/143‑1 (MQ=255)                                                                                                                                                                                                                                  tGAGATTG‑C‑TCTGGA‑AAGGTAT‑AAAAAAAACAGGACCAA‑AAGTCCTG‑‑TTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCCACAGCAACGGT‑GTCGTCTGAACCTACGCCCAGTTTGT‑TGTCAGAATCGATCTGGTtgatg                                                                                                                                                     <  1:869825/159‑1 (MQ=255)                                                                                                                                                                                                                                    agaTTG‑C‑TCTGGAAAAGGTAT‑AAAAAAAACAGGACCAA‑AAGTCCTG‑TTTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAGAAACTGGTAAACGATACCCACAGCAACGGT‑G‑CGTCTGAACCTACGCCCAGTTTGT‑TGTCAGa                                                                                                                                                                     <  1:869828/143‑1 (MQ=255)                                                                                                                                                                                                                                      aTTG‑C‑TCTGG‑‑AAGGTAT‑AAAAAAAA‑ACGACG‑ACAAGTCCTGTTTTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑a                                                                                                                                                                                                                                    <  1:121037/78‑1 (MQ=255)                                                                                                                                                                                                                                        tG‑C‑TCTGG‑‑AAGGTATAAAAAAAAACAGGACC‑A‑AAGTCCTG‑‑TTTTTCGGCA‑TTTAAACAAAAGAGGTGTGC‑TATTAGAAACTGGTAAACGATACCCACAGCAACGGT‑G‑CGTCTGAACCTACGCCCAGTTTGT‑TGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGAC‑GACGCAGACCGAAATCg  <  1:869827/300‑1 (MQ=255)                                                                                                                                                                                                                                        tG‑C‑TCTGG‑‑AAGGTAT‑‑‑AAAAAACAGGACC‑A‑AAGTCCTG‑TTTTTTCGGCATTTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCCACAGCAACGGT‑GTCGTCTGAACCTACGCCCAGTTTGT‑TGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGcc                                                                                                   >  1:869832/1‑201 (MQ=255)                                                                                                                                                                                                                                                    aaGGTAT‑‑AAAAAAACAGGACC‑A‑AAGTCCTG‑‑‑TTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCCACAGCAACGGT‑GTCGTCTGAACCTACGCCCAGTTTGT‑TGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTgg                                                                                                                                  <  1:869829/160‑1 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                      TTTATAGACACCAATCCCGAGGTTAATTTTGCCGGGACG‑TTCATCGGCACGAAACAGATCGGCCAGGCCCAGAATC‑GGGTCGGCAGGAGCGGCGGTAATGTTCTCAAACATGACGAGGTTCCATTATGGTTACAGAAGGGAAGTCC‑GCTATCAGGGTAACGGGAGATTTACAAAATTCCAAC‑TATT‑ACTGATGAAAACG‑CAGGCTGTTTTTGCAAGACGTGAGATTG‑C‑TCTGG‑‑AAGGTAT‑AAAAAAAACAGGACC‑A‑AAGTCCTG‑TTTTTTCGGCA‑TTT‑AAC‑AAAGAGGTGTGC‑TATTAG‑AACTGGTAAACGATACCCACAGCAACGGT‑GTCGTCTGAACCTACGCCCAGTTTGT‑TGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCG  >  NC_000913/985599‑986122

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