Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,015,515	(T)7→6	intergenic (+56/‑200)	pqiC → / → rmf	intermembrane transport lipoprotein PqiC/ribosome modulation factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,015,509	0	T	.	84.6%	46.1 / 0.9	26	intergenic (+50/‑206)	pqiC/rmf	intermembrane transport lipoprotein PqiC/ribosome modulation factor
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base . (7/15);  total (8/18)
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.33e-01

GAACCACCAGGGACAACTGATCAAACGTCCGTTCCGTCTGG‑AAGGAGTGCAAACTCAGGATGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGT‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑TTTTCTTTTCCACCAGAAACCAGTATGAGGGAAACGAGGCATGAAGAGACAAAAACGAGATCGCCTGGAACGGGCACATCAACGTGGTTATCAGGCCGGC  >  NC_000913/1015312‑1015774                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                  gAACCACCAGGGACAACTGATCAAACGTCCGTTCCGTCTGG‑AAGGAGTGCAAACTCAGGATGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGt‑t                                                                                                                                                                                                                                                                                 >  1:874146/1‑199 (MQ=255)       ccaGGGACAACTGATCAAACGTCCGTTCCGTCTGGAAAGGAGTGCAAACTCAGGATGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGATAAAAGACGT‑CTACC‑TTAATTAT‑AAGA‑TTGTAAATATAACCGTCTCCGGTATGTTGCCTGA‑GCGG‑‑TGTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑TATTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGa                                                                                                                                                    <  1:124316/313‑1 (MQ=255)       ccaGGGACAACTGATCAAACGTCCGTTCCGTCTGG‑AAGGAGTGCAAACTCAGGATGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGA‑GCGGT‑GTTTTTGTCTC‑TAACGTGCGG‑AAAAATTTTGTTCCTCTTCACA‑TATTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCgg                                                                                                                                             <  1:124317/321‑1 (MQ=255)                          cGTCCGTTCCGTCTGG‑AAGGAGTGCAAACTCAGGATGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCG‑‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAatat                                                                                                                                                                                           <  1:874147/253‑1 (MQ=255)                          cGTCCGTTCCGTCTGG‑AAGGAGTGCAAACTCAGGATGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCG‑‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAatat                                                                                                                                                                                           <  1:874148/253‑1 (MQ=255)                                  ccgtCTGG‑AAGGAGTGCAAACTCAGGATGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTCTTAACGTGCGG‑‑AAAA‑TTTGTTcc                                                                                                                                                                                                                                             >  1:874152/1‑198 (MQ=255)                                                        tCAGGATGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCG‑‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACtgtg                                                                                                                             <  1:874151/285‑1 (MQ=255)                                                              tGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑‑AAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAatat                                                                                                                                                                                           >  1:874156/1‑218 (MQ=255)                                                                gTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTcaca                                                                                                                                                                                                <  1:874153/212‑1 (MQ=255)                                                                 ttACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑‑AAAA‑TTTGTTCctactt                                                                                                                                                                                                                                        >  1:874157/1‑168 (MQ=255)                                                                               ttAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑ttt                                                                                                   >  1:874161/1‑288 (MQ=255)                                                                                tAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGtt                                                                                                                                                                                                                                               >  1:874162/1‑151 (MQ=255)                                                                                       cTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCG‑‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑tt                                                                                                    <  1:874155/278‑1 (MQ=255)                                                                                                 tCTGGAGTCAGG‑AGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGT‑AATATAACCGTCTCCGGTATGTTGCCTGAGGCG‑‑‑‑‑TTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑‑‑‑TTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGAc                                                                                                                                                                                        <  1:124327/180‑1 (MQ=255)                                                                                                       gTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGT‑AATATAACCGTCTCCGGTATGTTGCCTGA‑GCGGTGTTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCtt                                                                                                                                     <  1:874160/232‑1 (MQ=255)                                                                                                           ggAAGCCGCTTCTATTGCACAAAGAGAT‑AAAG‑CGTACTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACATT‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTTAAACTGTGATTGCACATTTAGTAATCACTG‑TTTTCTTTTCCACCAGAAACCAGTATGAGGGAAACGAGGCATg                                                           <  1:874164/304‑1 (MQ=255)                                                                                                                         ttGCACAAAGAGAT‑AAAG‑CGTACTACCTTTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTcac                                                                                                                                                                                                 <  1:874167/157‑1 (MQ=255)                                                                                                                                   agaT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCG‑‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑tt                                                                                                    <  1:874166/235‑1 (MQ=255)                                                                                                                                             cGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGaa                                                                                                                                                                             <  1:874168/157‑1 (MQ=255)                                                                                                                                                    cc‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAATTTTGTTCCTCTTCACA‑T‑TTTTTGTAc                                                                                                                                                                                                                         >  1:874170/1‑108 (MQ=255)                                                                                                                                                            ataaaaGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGG‑‑‑TTTTTGTCTC‑TAACGTGCGG‑‑AAAT‑TTGTTTCCTCTTCACATT‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTGTTTTTCTTTTccac                                                                                         >  1:124331/4‑226 (MQ=255)                                                                                                                                                                    tttGTAAATATAACCGTCTCCGGTATGTTGCCTGA‑GCGG‑‑TGTTTTGTCTC‑TAACGTGCGG‑AAAAA‑ATTGTTCCTCTTCACA‑T‑ATTTTGTACAACCGACATG‑CCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑TTTTCTTTTCCACCAGAAACCAGTATGAGGGAAACGAGGCATGAAGAGACAAAAACgag                                           <  1:124330/262‑1 (MQ=255)                                                                                                                                                                            tataACCGTCTCCGGTATGTTGCCTGAGGCG‑‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑TTTTCTTTTCCACCAGAAACCAg                                                                               <  1:874169/219‑1 (MQ=255)                                                                                                                                                                                         ggTATGTTGCCTGAGGCGG‑‑TTTTTTGTCTC‑TAACGTGCGGAAAAAA‑TTTGTTCCTCTTCACA‑‑‑TTTTTGTACAACCGACATG‑CCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑TTTTC‑TTTCCACCAG‑AACCAGTATGAGGGAAACGAGGCATGAAGAGACaaaacgag                                            <  1:874171/239‑5 (MQ=255)                                                                                                                                                                                         ggTATGTTGCCTGA‑GCGGT‑GTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑TTTTCTTTTCCACCAGAAACCAGTATGAGGGAAACGAGGCATGAAGAGA‑ACAAACGAGATCGCCTGGAACGGGCACATCAACGTGGTTATCAGGCCGGc  <  1:874172/283‑1 (MQ=255)                                                                                                                                                                                                       tttggT‑TGTTTTGCCTC‑TGACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACATT‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑TTTTCTTTTCCACCAGAAACCAGTATGAGGGAAACGAGGCATGAAGAGACAAAAACGAGATCGCCTGGAACGGGCACATCAACGTGGTTATCAGGcc     >  1:874177/4‑269 (MQ=255)                                                                                                                                                                                                       ggCGG‑‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACtgt                                                                                                                              <  1:874173/147‑1 (MQ=255)                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                  GAACCACCAGGGACAACTGATCAAACGTCCGTTCCGTCTGG‑AAGGAGTGCAAACTCAGGATGGTTACGATGAGATGGTTAAAGTGCTGGCCGGTGTCTGGAGTCAGGAAGCCGCTTCTATTGCAC‑AAGAGAT‑AAAG‑CGT‑CTACC‑TTAATTATAAAGATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGT‑TTTTTTGTCTC‑TAACGTGCGG‑AAAAA‑TTTGTTCCTCTTCACA‑T‑TTTTTGTACAACCGACATGCCCGTGTAGCTCACAAATATGACAGTGGCGTGAATTTTGCGCATTGACGGCAGTTATGATTCGCGGTATTGCTT‑AACTGTGATTGCACATTTAGTAATCACTG‑TTTTCTTTTCCACCAGAAACCAGTATGAGGGAAACGAGGCATGAAGAGACAAAAACGAGATCGCCTGGAACGGGCACATCAACGTGGTTATCAGGCCGGC  >  NC_000913/1015312‑1015774

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