Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,623,680	G→A	R146R (CGC→CGT)	lplA ←	lipoate‑‑protein ligase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,623,680	0	G	A	86.2%	58.5 / 3.4	30	R146R (CGC→CGT)	lplA	lipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (14/11);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

TGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTG‑GAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTT  >  NC_000913/4623437‑4623925                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      tGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACACGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTG‑GAAGCCGCGATCtt                                                                                                                                                                                                                                                >  1:188561/1‑251 (MQ=255)                         ttCCGCTTCTCCCCGCTCGCCATAAGGGGCGAAAAAGGCCTGGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGGGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGACAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGcctatcct                                                                                                                                                                                                                         <  2:130535/250‑7 (MQ=255)                                tctacacgcTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCttt                                                                                                                                                                                                                                               >  1:313807/7‑221 (MQ=255)                                tctacacgcTCGCCATAAGGGGCGAAAAAGGCCTTGGTTATGGCCGCGCAAACCTGCTCATGTGTGATCCCCGGCAGCAGCTCGGTTAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCttt                                                                                                                                                                                                                                               <  2:313807/215‑1 (MQ=255)                                       gcTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCg                                                                                                                                                                                                          >  1:348907/1‑251 (MQ=255)                                                  ggggCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGGTGGTCACGCGGGAACGCACTGATGTAATGCCTTGCGCCGCCAGTTGCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTTAGCAGCAAGGTGCCGTGGTG‑GAATCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGt                                                                                                                                                                                               <  1:36970/250‑1 (MQ=255)                                                         aaaaGGCCTTGGTTATGGCCACGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCCCGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGAAAGGTTGCCAGGCGCCTAGGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GGAGCCCCGATCTTTGGATTCGCGATAGGCCGCTCCTGAGACTTAGCGTTCGCCTTCGGCGGgttttatc                                                                                                                                                                                        >  2:128382/1‑247 (MQ=255)                                                              gCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTGaccacc                                                                                                                                                                                     >  1:201457/1‑249 (MQ=255)                                                                               gCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACTCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTTACCCCCCGATCGTTACGCCCGGa                                                                                                                                                                    >  1:124009/1‑249 (MQ=255)                                                                                 aaaCCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAACAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGCCTTTGCGGTCGCCTTCGGCGGTTTTCCCCACCAGATCGTATCGTCCGGCCGcc                                                                                                                                                                >  2:9489/1‑250 (MQ=255)                                                                                  aaCCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACACGATCGTTACGCCCGGACGCtt                                                                                                                                                               >  1:5892/1‑251 (MQ=255)                                                                                          cATGGGTGATCCCCGGCAGCAGATCGGAGAGGTTGGTCACGCGGGAACGCACTGAGGTACTGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGACGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTTCCCACCAGATCGTTGCGCCCGGACGCTTCCGCGCTg                                                                                                                                                       >  2:191468/1‑251 (MQ=255)                                                                                          cATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCTGTGGTG‑GAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTg                                                                                                                                                       >  1:72714/1‑251 (MQ=255)                                                                                                      ccGGCAGCAGCTCTGTGAGGTTGGTCACGCGTGCACGGACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGAGATGTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GACGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACGTTGCGGTCGCCTTCGGCGGTTGTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGTTGACGCCAAGcgcg                                                                                                                                           <  2:205750/251‑1 (MQ=255)                                                                                                                               cACGCGGGAACGCACTGAGGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACg                                                                                                                  <  1:363643/251‑1 (MQ=255)                                                                                                                                       aacgcactgaggTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCAACCCCAGAACGTTACGCCCGGGCGCTTCGGCGCTTGCGCCAAGAGCGGTCAGCGCATTGAGCACCATCGGCGttgaagtg                                                                                                          >  1:270265/12‑244 (MQ=255)                                                                                                                                            actgaggTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGtttt                                                                                                     <  1:97883/245‑1 (MQ=255)                                                                                                                                            aCCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTG‑GAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGACGCCTTCGACGGTTTTCACCACCAGCTCGTTACGTCAGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTTAGCCAAATCGACGAGGAGATAGtttt                                                                                                     >  2:162065/1‑251 (MQ=255)                                                                                                                                              tggggTAATGCCTTTCTCCTCCAGGTGCTTTTTATCCGGATTGCGATAGTTTGCCCGGCTGCTAAGGTCGGCATTGAGCATACAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGTTCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTTGAGATGGtttt                                                                                                     <  2:226441/245‑1 (MQ=255)                                                                                                                                              tgaggTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGCCGCTTCGGCGCTGACGCCAAGCGCGGTCAGCGCATTGAGCACAATCGACGTGGAGAtgggttttt                                                                                                   >  1:27255/5‑243 (MQ=255)                                                                                                                                                                  ccAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTgccagcc                                                                               >  1:39410/1‑251 (MQ=255)                                                                                                                                                                    aGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTG‑GAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCACGTTCAGCGCATTGAGCACAATCGCCGTGGAGATAGCTTTATCGTACTCCGGCTTgccagccc                                                                              >  2:141610/1‑249 (MQ=255)                                                                                                                                                                                   ggATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTTAGCACCATCGACGTGGAGATGGGTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAggg                                                              >  1:113825/1‑250 (MQ=255)                                                                                                                                                                                                ttGCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTTTCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGccc                                                       >  1:25167/1‑244 (MQ=255)                                                                                                                                                                                                  gCCAGGCGGCTAAGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGTCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTAATCCGGCTTGCCAGACATAAAGGTAAAGCAGGTATTTCCGAGATCGt                                                <  1:90585/250‑1 (MQ=255)                                                                                                                                                                                                             aaGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACAcc                                        <  1:273479/246‑1 (MQ=255)                                                                                                                                                                                                             aaGGTCGGCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGCTTAGGGCCTTGAGCACAAACGCCGTGGAGATGGGTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAACACGTATTGCCggagtcgtggaccaacagct                                    >  2:273479/2‑232 (MQ=255)                                                                                                                                                                                                                   ggCATTGAGCAGCAAGGTGCCGTGGTG‑GAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGCTTTTCACCACCAGATCGTTACGCCCGGACGCTTCCGCGCTGACGCCAAGCGCGTTCAGCGCCTTCAGAACCATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTTCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACAACGCGCCGCCAc                              >  1:381936/1‑251 (MQ=255)                                                                                                                                                                                                                    gaaTTGAGCACCAAGGTGGCG‑GGTGTGAAGCCACGTTATTTGGTTTCGCGATAGGCCTATCCTGAGACTTTGCGGTCGCCGTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACTCTTCGGCGCGTACGACAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTATTCTGTCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGTGATCGTGGAACACCCCGCCGCCAcc                             <  2:251072/249‑1 (MQ=255)                                                                                                                                                                                                                                                aaGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGtt  <  1:290664/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      TGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAGGTGCCGTGGTG‑GAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTT  >  NC_000913/4623437‑4623925

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