Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,866,420	C→T	G18D (GGC→GAC)	mipA ←	MltA‑interacting protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,866,420	0	C	T	95.2%	57.0 / ‑5.9	21	G18D (GGC→GAC)	mipA	MltA‑interacting protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (8/12);  total (8/13)
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.95e-01

ATTTGGTGATCGCCACTGTCTTTCGCTTTGAAGTAAAGCGGCGACCAGTAAGCGGTAATTGAAAGTTTATCCGTTGCGTCATTCCACAGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGCCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACAAAGTTTACCCGTCAATACATCGAAACCCAACCTTTTTACGTTTTCATTTTTTAAAGTAACTGTTTAATTTTCCTGACGCGGATGACACCGCGCTTACAGCCAAATGAATTTAGCGTTACTGGCGAGCCTGGTCTTTACAT  >  NC_000913/1866177‑1866650                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       aTTTGGTGATCGCCACTGTCTTTCGCTTTGAAGTAAAGCGGCGACCAGTAAGCGGTAATTGAAAGTTTATCCGTTGCGTCATTCCACAGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGa                                                                                                                                                                                                                                 <  1:356380/251‑1 (MQ=255)               actttCTTTCGCTTGGAAGTAAAGCGGCGACCAGTAAGCGGTAATTGAAAGTTTACCCGTTGCGTCATTCCACAGGTAGTAACCACCACCGAAGCCACGGAACCAGAAGTTATCGCCTTCATAGGTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCa                                                                                                                                                                                                                   <  2:338767/247‑1 (MQ=255)                                    aaGCGGCGACCAGTAAGCGGTAATTGAAAGTTTATCCGTTGCGTCATTCCACAGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAagtttga                                                                                                                                                                                              <  2:342114/251‑1 (MQ=255)                                               aGTAAGCGGTAATTGAAAGTTTATCCGTTGCGTCATTCCACAGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTcac                                                                                                                                                                                    >  1:236512/1‑250 (MQ=255)                                                    gCGGTAATTGAAAGTTTATCCGTTGCGTCATTCCACAGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTaa                                                                                                                                                                              <  1:159523/251‑1 (MQ=255)                                                                  tttgtacGTTGCGTCATTCCACCGGTAGTAGCCACCACATAATCTACGTAACCAGAATTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTATTCAACGACCCCTAAGCTTGCGCCCAGGGGAATTTTTCCTTCAGAGTGCGCTGAGTCTGACGACGTTACGTTAAGAACTCCAATTACCAGTAGTTTGAGTTTTGTCACAATTACTAATTCCTTAAAc                                                                                                                                                                 >  2:109187/7‑250 (MQ=255)                                                                              gTCATTCCACAGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACGCCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCgg                                                                                                                                                    <  1:80559/251‑1 (MQ=255)                                                                                cATTCCACAGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGc                                                                                                                                                  >  1:269859/1‑251 (MQ=255)                                                                                 aTTCCACAGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCg                                                                                                                                                 >  1:479545/1‑251 (MQ=255)                                                                                       caGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACaa                                                                                                                                             >  1:286021/1‑249 (MQ=255)                                                                                          ggagtaACCACCACCTAAGCCACGGAACCAGAAGTTATCGCATTCATAGTTGATTTCCGGTACGGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACAAAGttt                                                                                                                                        <  2:560115/249‑1 (MQ=255)                                                                                                                       aGAAGTTATCGCCTTCATAGGTGAGTACCGGGACTGGGTAAACATCGGTATCGTAAGCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCGTCAGCGTGCGCTACGCCTGCAGACGTTGTGATAAGCACTCCAAGTGCCAGAAGTTTGAGTGTGGTCCCAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACAAAGTTTACCCGTCAATACATCGAAACCCAACCttt                                                                                                           <  2:144357/251‑1 (MQ=255)                                                                                                                          aGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACAAAGTTTACCCGTCAATACATCGAAACCCAACCTTTTTa                                                                                                        <  2:202219/251‑1 (MQ=255)                                                                                                                                     tCATAGTTGATTTCCTGTCCTTGGTAACCTTCGGTTTCGTAATCTTTATATGTGTGTTCAACGACTCCTTCGCCTGCTCCCTGGTACATTTTGCCTTCCGCGTGCTCTTCGTCTTCATACGTTTCGATAAGCACTCCAAGTTCCAGAAGTTTGATTTTGTTCACATTTCATCATTCCTTAATCAAATGTTTATCGTGCGACAAAGTTTACCCGTCAATACATCGAAACCCATCCTTTTTACGTTTTCAt                                                                                               <  2:224736/249‑1 (MQ=255)                                                                                                                                               ggaCCGGTACTGCGTAATCATCGGTTGCTTAATTTTTATATTGGTGTTAAACGACACCTCCGCCATCTCCCAGGGAAAATTTTCCTTCAGCGTGCGCTTCGTTTGCTGATGTTGCGATAAGCCCTCCAAGTGCCTGAAGGTTGAGTTTGGTCACAATTAATCATTCCTAAAACAAATGTTTAGCGGGCGACAAAGTTTACCCGTCAATACATCTCAACCCATCCTTTGTACGTTTGCAGTTTTTAAAGTaa                                                                                   <  2:402560/249‑1 (MQ=255)                                                                                                                                                   cGGTACTGGGGAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACCCCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACCATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACAAAGTTTACCAGTCAATACATCGAAACCCAACCTTTTTACGTTTTCATTTTTTAAAGTAACTGt                                                                               >  2:83096/1‑251 (MQ=255)                                                                                                                                                                                                      acCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACAAAGTTTACCCGTCAATACATCGAAACCCAACCTTTTTACGTTTTCATTTTTTAAAGTAACTGTTTAATTTTCCTGACGCGGATGACACCGCGCTTACAGCCAAATGAATTTAGc                            <  1:556515/251‑1 (MQ=255)                                                                                                                                                                                                               tGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACAAAGTTTACCCGTCAATACATCGAAACCCAACCTTTTTACGTTTTCATTTTTTAAAGTAACTGTTTAATTTTCCTGACGCGGATGACACCGCGCTTACAGCCAAATGAATTTAGCGTTACTGGc                   >  1:314733/1‑251 (MQ=255)                                                                                                                                                                                                                            aatgtACCGTCAGCGTGCGCTACGTCTGCAGACTTTGCGATAAGCTCTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATTTTTAGCGGGCGACAAAGTTTACCCGGCAATACATCGAAACCCACCCTTTTTACGTTTTCATTTTTTAAAGTAACTGTTTAATTTTCCTGACGCGGATGACACCGCGCTTACAGCCAAATGAATTTAGCGTTACTGGCGAGCCTGGTCttt      <  2:479545/247‑1 (MQ=255)                                                                                                                                                                                                                            aaTTTGCCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGAAAAAGTTTTCCCGTCAATTCATCGAAACCCAACCTTTTTACGTTTTCATTTTTTAAAGTAACTGTTTAATTTTCCTGACGCGGATGACACCGCGCTTACAGCCAAATGAATTTAGCGTTACTGGCGAGCCTGGTCttt      <  1:275319/251‑1 (MQ=255)                                                                                                                                                                                                                                tACCTTCAGCGTGCGCTACGTCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACAAAGTTTACCCGTCAATACATCGAAACCCAACCTTTTTACGTTTTCATTTTTTAAAGTAACTGTTTAATTTTCCTGACGCGGATGACACCGCGCTTACAGCCAAATGAATTTAGCGTTACTGGCGAGCCTGGTCTTTACat  >  1:570209/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       ATTTGGTGATCGCCACTGTCTTTCGCTTTGAAGTAAAGCGGCGACCAGTAAGCGGTAATTGAAAGTTTATCCGTTGCGTCATTCCACAGGTAGTAACCACCACCTAAGCCACGGAACCAGAAGTTATCGCCTTCATAGTTGATTACCGGTACTGGGTAAACATCGGTATCGTAATCTTTATATGGGTGTTCAACGACACCTACGCCTGCGCCCAGGGAAAATTTACCTTCAGCGTGCGCTACGCCTGCAGACGTTGCGATAAGCACTCCAAGTGCCAGAAGTTTGAGTTTGGTCACAATTAATCATTCCTTAAACAAATGTTTAGCGGGCGACAAAGTTTACCCGTCAATACATCGAAACCCAACCTTTTTACGTTTTCATTTTTTAAAGTAACTGTTTAATTTTCCTGACGCGGATGACACCGCGCTTACAGCCAAATGAATTTAGCGTTACTGGCGAGCCTGGTCTTTACAT  >  NC_000913/1866177‑1866650

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