Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,171,417	G→A	L850L (CTG→TTG)	mfd ←	transcription‑repair coupling factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,171,417	0	G	A	97.0%	93.5 / NA	34	L850L (CTG→TTG)	mfd	transcription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (19/13);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCCG  >  NC_000913/1171178‑1171661                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     aCGTGTTTGTGCCTCGGTATTCATCGCTTTCTGATGCGTTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCTTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTTCAAAACAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATtcgcgctcg                                                                                                                                                                                                                                            <  2:288469/250‑1 (MQ=255)          ggCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGcc                                                                                                                                                                                                                                  <  1:178680/250‑1 (MQ=255)                      cagcGCTTTTGGAAGCGGTGTCAGCAACCATGCAGATGCCGGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTGCAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATg                                                                                                                                                                                                                       <  2:36717/247‑1 (MQ=255)                              ttgtacgctgtggccgCAACCATGCATCTGTCTGATGATGCGAACGTCCGACGCGACCGCGTAACTCGGGCATCTGCGCCAGAGCGAAGTGATCCTCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCTCTGTTTCGATAATGGTTGTACAAACCAGCACGGTGAAACGTTGATGATGGAAATCATTCATCACCCGTAGCAAGTT‑GCGCTCGCGCTGCTGCCCGTGAACGATGGCGATcag                                                                                                                                                                                                               <  2:259130/235‑3 (MQ=255)                                    gCGGTGTCAGCAACCATTCTTATTCCCGATGATGCGAACGGCCGACGCGTCTGCGTAACTGGTGCAGCTGGTCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCTATCCTGGTTTCGATAATGGTTGTCCAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTTCCCGTGACCGATGGCGATCCGTGCTTCt                                                                                                                                                                                                        <  2:35469/251‑1 (MQ=255)                                          tCAGCAACCATACATATGCCTGATGATGCGAACGTCCGAAGCGACCGCGCAACTGAAGCAACTGCACCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATATCGATCCCGGTTTCAATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGCTGTACATCATTCATCACCAGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACAGTTGGCGATCCGCGCATCTggaccc                                                                                                                                                                                                  >  2:19978/1‑247 (MQ=255)                                                  cATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGc                                                                                                                                                                                          >  1:30635/1‑251 (MQ=255)                                                       atttGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGGGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGCCCGATGGCGATCCGCGCTTCTGGCACCAGTTCGGCCAGcc                                                                                                                                                                                     <  2:30635/248‑1 (MQ=255)                                                                    tgCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCttt                                                                                                                                                                        >  1:243015/1‑251 (MQ=255)                                                                    ggCGAACGTCCGCCGCGACCGCGTTACTGGTGCAGCGGCGCCAGTCCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGCCGGGATGTCGAGCCCGGTTTCGATAACGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATGCGCGCTCGCGAATCGGCCCTTGACCGGTGTCGATCCGCGCTTCTGGCACCAGTTCGGCCAGCCGTGCGGCGGCttt                                                                                                                                                                        <  2:174598/250‑1 (MQ=255)                                                                                aCGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGCTGGAAATCATTCATCACCCGT‑TCCACTTCGCGCTCGCGCATCTGCCAGTGACCGATGGCGATCCGCGCTTCTGGCACCCGTACTGCCAGCCGTTCGTCGGCTTTCTGAATGTTTTc                                                                                                                                                            >  2:209564/1‑251 (MQ=255)                                                                                aCGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTc                                                                                                                                                            >  1:77340/1‑251 (MQ=255)                                                                                   cgACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCcac                                                                                                                                                         >  1:374104/1‑251 (MQ=255)                                                                                     aCCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGCTGGCGATCCGCGCTTCTGGCCCCCGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCCcat                                                                                                                                                       >  2:324097/1‑251 (MQ=255)                                                                                         cgTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCAtt                                                                                                                                                   >  1:87037/1‑251 (MQ=255)                                                                                                     aGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGataata                                                                                                                                       <  1:377458/251‑1 (MQ=255)                                                                                                         gcgcCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATCCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAAc                                                                                                                                    <  1:19978/250‑1 (MQ=255)                                                                                                              aGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGAGCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGtcc                                                                                                                              <  1:324097/251‑1 (MQ=255)                                                                                                                   gAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGGTTTCCACATCATTGGAGAGATAATAAACCTGTCCTCcgc                                                                                                                         >  1:262804/1‑251 (MQ=255)                                                                                                                                           aataGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGAGCATCTGCCCGTGACCGATGGCGATCCGCGATTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCt                                                                                                 >  1:153915/1‑251 (MQ=255)                                                                                                                                                tgtTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCCCGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACTGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCgg                                                                                            >  1:364196/1‑251 (MQ=255)                                                                                                                                                             gATGTCGCTCTCTGTTTCGATAATGGTTGTACAACCAAGCAAGTTGAAACGTTGATGATGGAAATCATTCATCACCAGT‑TCCAATTCGCTCTCTCGCATCTGCCCGTGACCGATGTCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTATGAATGTTTTCCACATCATTGTAGCGATAATAAACCTGACCTCCGCGCAAAATTTCACGCAGGATCCCGTCCAGGACCTCCAttcttt                                                                               >  2:16029/1‑246 (MQ=255)                                                                                                                                                              aTGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCATCACGGACCACCATGCAGTc                                                                              >  1:278291/1‑251 (MQ=255)                                                                                                                                                                    aTCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGCAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTc                                                                        >  2:174064/1‑251 (MQ=255)                                                                                                                                                                      cccGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCAc                                                                      >  1:97579/1‑251 (MQ=255)                                                                                                                                                                                aTAATGGTTGTACAAACCAGCACGATGAAACGTTGATGATGGACATAATTCATCACCCGA‑TCCAATCCGCGATCGCGCATCTGCCCGTGACCGATGGCGATCCGCGATTCTGGCACCAGTTCTGCCCGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGAGAATAAACCTGTCCTCAGCGCAAAATTTCACGCAGGATCGACTCCCGGCCCACCATGATGTCATaatccaggaccacggt                                                            >  2:47636/1‑236 (MQ=255)                                                                                                                                                                                            cAAACCCGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGcc                                                  >  1:55503/1‑249 (MQ=255)                                                                                                                                                                                                 cAGCACGTTGAAACGTTGATGATGGAAATCCTTCATCCCCCGT‑CCCAATTCGCGCTCGCGCACCTGCCCGTGTCCGATGGCGATCCGCGCTTCTGTCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAgacgac                                           <  2:104129/251‑1 (MQ=255)                                                                                                                                                                                                    cACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGAcggg                                        >  1:12737/1‑251 (MQ=255)                                                                                                                                                                                                    cACGTGGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAAGTCGCGCTCGCGCATCTGCCCTTGACCGATGGCGATCCGCGCTTCTGGTACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATGGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGAcggg                                        <  2:87037/251‑1 (MQ=255)                                                                                                                                                                                                                  gatgatGGAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGc                           <  1:117788/250‑1 (MQ=255)                                                                                                                                                                                                                        tgAAATCATTCATCACCCGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTAACCGATAGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCCTTGTCGAGATCATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGACATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGTCGGGCGGCGTGGCGATAAt                     >  2:38641/2‑250 (MQ=255)                                                                                                                                                                                                                                       cccGT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCAtt     <  1:209564/251‑1 (MQ=255)                                                                                                                                                                                                                                          gT‑TCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCCg  >  1:360869/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     ACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGT‑TCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCCG  >  NC_000913/1171178‑1171661

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