Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,012,789	G→A	L10L (CTA→TTA)	metR ←	DNA‑binding transcriptional dual regulator MetR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,012,789	0	G	A	96.9%	97.3 / NA	32	L10L (CTA→TTA)	metR	DNA‑binding transcriptional dual regulator MetR
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (11/20);  total (11/21)
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

GTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACATAT  >  NC_000913/4012551‑4013035                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       gTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAggcg                                                                                                                                                                                                                                            >  1:353263/1‑249 (MQ=255)            ttCATTGCAGGCTTGCAGGGCCTGGCTACATTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGCTTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGATGTTCCAGCTCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGGTGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCg                                                                                                                                                                                                                                  >  2:114229/1‑250 (MQ=255)             tCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGa                                                                                                                                                                                                                                 <  1:323546/250‑1 (MQ=255)              cATTGCAGGCTTGCAGGGCCTGGTTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGAt                                                                                                                                                                                                                                <  2:71326/250‑1 (MQ=255)                                   gggtAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCg                                                                                                                                                                                                            <  2:199186/246‑1 (MQ=255)                                           tGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGaata                                                                                                                                                                                                   <  1:224172/250‑1 (MQ=255)                                                           ttGCCAGTTGCAACATGATTTCTCCCTGCGGTGTAAAGCGTAGCTGCTGGCTCTTACGCTCAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCTAGCGAGCCGCAGTTCCGCAACGCTTTTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGaa                                                                                                                                                                                   <  2:380547/250‑1 (MQ=255)                                                                  tagCAAAAGGATTTCTCCCTGCGGTGTAAAGCGTAGAGGCTGGCTCTTACGAACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAAAGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTACGACACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATAATGAAAGTCCTTCACTTAGGCATGAATAATTTGCGCTTGAGGAATATACAGt                                                                                                                                                                           >  2:43397/3‑251 (MQ=255)                                                                                       gggTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATgtgt                                                                                                                                                      <  1:70438/250‑1 (MQ=255)                                                                                       cGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATgtgt                                                                                                                                                      >  1:162814/1‑251 (MQ=255)                                                                                                      gggctggctCTTACGCACAAATAGCCGGAAGCCACGGCGTTGTTCCAGATCGCTAACCTGGTGAGACAGGGCGGTTTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGGGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGc                                                                                                                                       <  2:353263/250‑1 (MQ=255)                                                                                                                  aCGCACAAATAGCCGGAAGCCATGGCGTTGTTCCAGATCGCTAAACGGGTGAGACAGGGCGGATTG‑CTTCTGAGGCACCGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCg                                                                                                                           <  2:305604/251‑1 (MQ=255)                                                                                                                       caAATAGCGGGAAGCCAAGGCGTTTTCCCGGATCGCGAAGTTGTTGAGACTGGGCGG‑TTGTCGTCTGATGCACCGTCGCCGCAGCGTCTGCTGTCGAGCCGCAGTTCCGCAACGCTTTTAACTTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCttt                                                                                                                      <  2:28341/251‑1 (MQ=255)                                                                                                                                         ggCGTTTTTCCAGATCGTTACCCTGGGGAGAGAGGGCGTATTG‑CGTCTGGTGCAACGTGGCCGCAGCGTGTGCCAGCTTGCCGCAGTTCCGCAACGCTTGTAACGTTTTCATTTTTTTTCCTTTTATCATAAAAGTCCTTCACTTCGGCATGATTAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATTGATGTGTAAGCATCTGGACGGCTAAAAGCCTTCGTCTTTTAAATTTATGGTGCGtcg                                                                                                    <  2:27772/251‑3 (MQ=255)                                                                                                                                                 tCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGtt                                                                                            >  1:349097/1‑251 (MQ=255)                                                                                                                                                    aGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCGGCGGCTGGGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTtct                                                                                         <  2:38840/251‑1 (MQ=255)                                                                                                                                                        cGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCAc                                                                                     >  2:314035/1‑251 (MQ=255)                                                                                                                                                                       aCAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCCCCCCGGTCacttactt                                                                      >  1:194727/1‑251 (MQ=255)                                                                                                                                                                                        cTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAg                                                                                                                                                                            <  1:301069/133‑1 (MQ=255)                                                                                                                                                                                        cTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAg                                                                                                                                                                            >  2:301069/1‑133 (MQ=255)                                                                                                                                                                                         tGATGCACCGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCgggg                                                      <  1:159459/250‑1 (MQ=255)                                                                                                                                                                                          aaTGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGAAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGAt                                                    <  1:162915/250‑1 (MQ=255)                                                                                                                                                                                                          aGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAAtt                                                                                                                                                                <  1:159888/127‑1 (MQ=255)                                                                                                                                                                                                          aGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAAtt                                                                                                                                                                >  2:159888/1‑127 (MQ=255)                                                                                                                                                                                                                  cgatcgtgcCGCAGTTCCGCATCGCTTGTAACGTTTTCAGGTTTTTTTCTTCGATCATGAAAGTCTTTCTCTTCGGCATTAATAATTTGCGCTTGAGGAATATACAGTATCCGGCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGACTTCCCTa                             <  2:76153/243‑1 (MQ=255)                                                                                                                                                                                                                    ggcgagcAGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCGATAATGAAAGTCCTTCCCTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAAt                          <  1:309455/250‑1 (MQ=255)                                                                                                                                                                                                                     gcgagcGGCAGTTCCGCAACGCTTGTAACGTTTTCAGGTGTTTTACTTCAATCATAAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTGGAGGAAAATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAAtt                         <  1:308873/251‑1 (MQ=255)                                                                                                                                                                                                                          ccgcagcggcggagctcTTGTAGCGTTTTCAGTTGTTGTACTTTGATCACGAAAGTCGTGCATTTCTGCATGAATAGTTTGCACTTGAGGAATGTGCGCTATCCGCCAATTACGGATTTGTAACCATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGTCGGCGTTTCTTCGCCCCTGTCACATACTTCAGTAAGCTCCCGGGGATGAATAAACGTGCCGGCTTCCCTAAATTCaaaa                    <  2:222798/235‑1 (MQ=255)                                                                                                                                                                                                                                  ccGCAACGCTTGTAACGTTTTAAGGTGTTTTACTTCGATCATGAAATTCCTTCACTGCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAg             <  2:9327/250‑1 (MQ=255)                                                                                                                                                                                                                                    ccAACGCTTGTTACGTTTTGAGGTGTTTTGCTTCGATCATAAAAGTACTTCACTGCGGAATGAAGAAGTTGCGTTGGAGGAATATACAGTATCCGCCAATTATGGATGTGTAAACAGCTGGACGGCTAAAATCCTTCGTCTTTGAAATTTATGGTGCGTGGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGCGATGAATAAACTTGCCGCCTTCCCTAGATTCAAAATCCATAGGAt          <  2:253068/250‑1 (MQ=255)                                                                                                                                                                                                                                      aaCGCTTGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGNttt        >  2:365854/1‑251 (MQ=255)                                                                                                                                                                                                                                            tGTAACGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACatat  >  1:276086/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       GTCTGCTGCGGTTCATTGCAGGCTTGCAGGGCCTGGCTAATTTGCGGCAGTACCTGGTTTGCCAGTTGCAACAGGATTTCTCCCTGCGGTGTAAAGCGTAGCGGCTGGCTCTTACGCACAAATAGCCGGAAGCCAAGGCGTTGTTCCAGATCGCTAAACTGGTGAGACAGGGCGGATTG‑CGTCTGATGCAACGTCGCCGCAGCGGCTGCGAGCGAGCCGCAGTTCCGCAACGCTTGTAGCGTTTTCAGGTGTTTTACTTCGATCATGAAAGTCCTTCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGATGAATAAACTTGCCGCCTTCCCTAAATTCAAAATCCATAGGATTTACATAT  >  NC_000913/4012551‑4013035

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: ‑