Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,499,393	T→G	intergenic (+192/‑55)	cysG → / → yhfL	siroheme synthase/DUF4223 domain‑containing lipoprotein YhfL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,499,393	0	T	G	97.4%	136.1 / NA	38	intergenic (+192/‑55)	cysG/yhfL	siroheme synthase/DUF4223 domain‑containing lipoprotein YhfL
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (18/19);  total (18/20)
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

CGGGTTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTGCTTATATTCATATAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGC  >  NC_000913/3499154‑3499634                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  cGGGTTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcg                                                                                                                                                                                                                                        >  2:390762‑M1/1‑247 (MQ=255)       gttgGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgacc                                                                                                                                                                                                                                     >  1:228048‑M1/1‑241 (MQ=255)          ggCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataa                                                                                                                                                                                                                               <  1:175080‑M1/251‑14 (MQ=255)                                           cATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatc                                                                                                                                                                                              <  2:144235‑M1/251‑47 (MQ=255)                                                    aaCCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTTctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgct                                                                                                                                                                                      <  2:95553‑M1/250‑55 (MQ=255)                                                    aaCCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgctt                                                                                                                                                                                     <  1:36122‑M1/251‑56 (MQ=255)                                                    aaCCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgct                                                                                                                                                                                      <  2:108385‑M1/250‑55 (MQ=255)                                                                             tcttAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcct                                                                                                                                                            >  1:426928‑M1/1‑171 (MQ=255)                                                                               ttAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgc                                                                                                                                                          >  1:491628‑M1/1‑169 (MQ=255)                                                                                tAATGACTCTTTTTATATATTCAGCAAANAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcct                                                                                                                                                            <  2:228048‑M1/248‑81 (MQ=255)                                                                                   tGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtc                                                                                                                                                       >  2:268270‑M1/1‑165 (MQ=255)                                                                                        ctTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgc                                                                                                                                                                 <  2:15196‑M1/235‑76 (MQ=255)                                                                                        ctTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgc                                                                                                                                                                 >  1:15196‑M1/1‑160 (MQ=255)                                                                                                  atTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaac                                                                                                                                       >  1:206800‑M1/1‑150 (MQ=255)                                                                                                       gCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgnggcctgcgtcaggccaacggatgaacggcat                                                                                                                                  >  2:242655‑M1/1‑145 (MQ=255)                                                                                                         aaaTAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattc                                                                                                                                <  2:31896‑M1/251‑109 (MQ=255)                                                                                                                  aTAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgcc                                                                                                                         >  1:311877‑M1/1‑134 (MQ=255)                                                                                                                    aGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagct                                                                                                                     >  1:23086‑M1/1‑132 (MQ=255)                                                                                                                                      atCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaat                                                                                                    <  2:108556‑M1/250‑137 (MQ=255)                                                                                                                                          gCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgatt                                                                                               <  1:307536‑M1/251‑142 (MQ=255)                                                                                                                                             aattgtttatattcCCCGTTACAGCGACTTACAAAATTAATTCATTCGCTTACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGTGCAGTCTTCGTATATTCctcgaccagataactgccattgcggccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggtcaacggaagaacggcattcatttgccagcgcatctatggttttcaatcgattac                                                                                             <  2:427288‑M1/244‑144 (MQ=255)                                                                                                                                                gggTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgc                                                                                          >  1:135501‑M1/1‑104 (MQ=255)                                                                                                                                                  gTTATATTCACCGTTACAGCGAATTTCAAAATTAAGTCTTTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcact                                                                                       <  2:299715‑M1/251‑150 (MQ=255)                                                                                                                                                       atTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaccggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttc                                                                                   >  2:337368‑M1/1‑97 (MQ=255)                                                                                                                                                         tCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcac                                                                                 <  1:268270‑M1/250‑156 (MQ=255)                                                                                                                                                                          tACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATTTTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgacccgataactgccattgcgcccccggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggcc                                                                <  2:93088‑M1/250‑173 (MQ=255)                                                                                                                                                                            cAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagc                                                             >  2:82521‑M1/1‑76 (MQ=255)                                                                                                                                                                                       cATTCGCTCACTGACGCCAATTAAATATTCCTTTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgac                                                  <  2:135501‑M1/251‑187 (MQ=255)                                                                                                                                                                                        aTTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacg                                                 <  1:390762‑M1/251‑188 (MQ=255)                                                                                                                                                                                                 cTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaacgcacttttcacccaacaaataacgggccagcctcgtgatgacgacgccttct                                        >  2:79095‑M1/1‑55 (MQ=255)                                                                                                                                                                                                             aaaTATTCCTGTATCGTTCTTAACGACCAGTCTTCGTAAATTCctcgaccagaaaactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacacggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgctcttttcacccaacaaataacggaccagcgtcgtgatgacgacgccttcttttttgataaac                            <  2:151079‑M1/251‑209 (MQ=255)                                                                                                                                                                                                                       gTATCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgc                   <  2:488141‑M1/250‑218 (MQ=255)                                                                                                                                                                                                                         aTCGTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggctgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgtcgacgccttcttttttgataaacgttcgacgcat                 >  2:365531‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                            gTTCTTAACGAGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataat              >  1:140050‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                             gtcTTAATGAGCAGTCTTCTTTTATTCctcgaccaggtaactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgctgcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgtttaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataatctttcgacgcataat              <  2:214461‑M1/248‑223 (MQ=255)                                                                                                                                                                                                                                      aGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataataaattttcagt   >  2:369706‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                      aGCAGTCTTCGTATATTCctcgaccagataactgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatgggtttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataataaattttcagt   >  2:384924‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                       gCAGTCTTCGTATATTCctcgaccatatacctgccattgcgcccacggcgttcaatccgtatcgctccgcttgactccagcgttttcaacgcggcctgcgtcaggccaacggatgaacggcattcatttgccagctcatctatggttttcaatcgattaccgcacttttcacccaacaaataacgggccagcgtcgtgatgacgacgccttcttttttgataaacgttcgacgcataataaattttcagta  <  2:23086‑M1/251‑235 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  CGGGTTGTTGGCCTGCGCGATAAACTGAACTGGTTCTCCAACCATTAATTTAACCCGGGCCAGAGAATTCTGGCCTTCTTAATGACTCTTTTTATATATTCAGCAAATAAAACATAGCCCCTTAATAAATATTATCTGCTAATGGGTTATATTCACCGTTACAGCGACTTACAAAATTAAGTCATTCGCTCACTGACGCCAATTAAATATTCCTGTATCGTTCTTAACGAGCAGTCTGCTTATATTCATATAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGC  >  NC_000913/3499154‑3499634

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