Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,696,861	T→C	L64L (TTA→TTG)	uidR ←	transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,696,861	0	T	C	92.6%	86.7 / ‑0.8	27	L64L (TTA→TTG)	uidR	transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (12/13);  total (13/14)
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.00e+00

CATAATTTCAACCACCAGCGCCCGTTGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGTAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAATCTCTCTGGCAGCATTGAGGAT  >  NC_000913/1696726‑1697007                                                                                                                                        |                                                                                                                                                   cATAATTTCAACCACCAGCGCCCGTTGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTg                                                                                                                                      >  2:368355/1‑150 (MQ=255)     aTTTCAACCACCAGCGCCCGTTGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATc                                                                                                                                  >  2:523887/1‑150 (MQ=255)         cAACCACCAGCGCCCGTTGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAgg                                                                                                                              <  2:180933/150‑1 (MQ=255)          aaCCACCAGCGCCCGTTGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGct                                                                                                                            >  2:511748/1‑151 (MQ=255)          aaCCACCAGCGCCCGTTGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGct                                                                                                                            >  2:269877/1‑151 (MQ=255)          aaCCACCAGCGCCCGTTGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGc                                                                                                                             <  1:388660/150‑1 (MQ=255)                         ttGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGaa                                                                                                              >  2:353579/1‑150 (MQ=255)                          tGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGaaa                                                                                                             <  1:133745/150‑1 (MQ=255)                                                     aCAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATgg                                                                                 <  2:337870/151‑1 (MQ=255)                                                                gACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGAttt                                                                       >  2:597302/1‑150 (MQ=255)                                                                              cGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCAt                                                        <  1:137542/151‑1 (MQ=255)                                                                               gAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATc                                                       <  1:348758/151‑1 (MQ=255)                                                                                         ccTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACt                                               <  1:627505/149‑1 (MQ=255)                                                                                          cTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACtgtg                                            >  2:109424/1‑151 (MQ=255)                                                                                             cAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGaa                                          <  2:36657/150‑1 (MQ=255)                                                                                                  ggTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCAt                                    <  2:63826/151‑1 (MQ=255)                                                                                                     tCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCAtttt                                 >  1:133402/1‑151 (MQ=255)                                                                                                       ccGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGTAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCt                               <  1:16551/151‑1 (MQ=255)                                                                                                              cGGGCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGaaaaaa                        <  2:133399/151‑1 (MQ=255)                                                                                                                 gCCAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAAtct                     <  1:561812/151‑1 (MQ=255)                                                                                                                  ccAGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAAtctc                    <  1:597182/151‑1 (MQ=255)                                                                                                                    aGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAAtctctc                  >  2:631434/1‑151 (MQ=255)                                                                                                                    aGCGCCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAAtctctc                  >  2:387303/1‑151 (MQ=255)                                                                                                                      cgcCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAATCTCTCTg                >  2:68235/1‑151 (MQ=255)                                                                                                                      cgcCCTCTCCTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAATCTCTCTg                >  1:279614/1‑151 (MQ=255)                                                                                                                              ccTGGTCCTGCAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAATCTCTCTGGCAGCAt         <  1:589248/150‑1 (MQ=255)                                                                                                                                    ccTGTAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAATCTCTCTGGCAGCATTGAGGAt  >  2:260539/1‑151 (MQ=255)                                                                                                                                        |                                                                                                                                                   CATAATTTCAACCACCAGCGCCCGTTGTCCAAAGGCTTCATGGGTGAGAGAGACAATGGACTCGACCATATAGTCAACGAAATGAATCCCTTCAATCGGTTCCCGGAAACGGGCCAGCGCCCTCTCCTGGTCCTGTAAGATAATCGCCTGAATCAAGGCTTCTTTGGAGATGAAATGGTGATAGAGCGTCCCGGGACTAATGGCGCAAGATTTACAGATGGCTTTCATCGAGGCACTGTGAAATCCATTTTCTGAAAAAATCTCTCTGGCAGCATTGAGGAT  >  NC_000913/1696726‑1697007

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