Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	646,838	G→A	intergenic (+192/‑395)	lplA2 → / → mvaK1	lipoate‑‑protein ligase A/mevalonate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	646,838	0	G	A	96.4%	95.1 / ‑6.1	28	intergenic (+192/‑395)	lplA2/mvaK1	lipoate‑‑protein ligase A/mevalonate kinase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (11/16);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGAAAGTATCCATAAAAAGCATAAATATGTTATTGACCTAATATTTAAATTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAGGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAATTATTATTTGTAGTATTATAGTTAGGGTGGAAAATTATAAATTGT  >  CP000730/646693‑646968                                                                                                                                                  |                                                                                                                                   ttttGAAAGTATCCATAAAAAGCATAAATATGTTATTGACCTAATATTTAAATTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAagg                                                                                                                                  <  1:56624/148‑1 (MQ=255)       aaGTATCCATAAAAAGCATAAATATGGTATTGACCTAATATTTAAAGTGAACATAATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATATCATGTTTTACGTATTAAATAGATTGAGAATAAGGtatata                                                                                                                            >  2:263462/1‑148 (MQ=255)        aGTATCCATAAAAAGCATAAATATGTTATTGACCTAATATTTAAATTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGtatatat                                                                                                                           <  1:217409/148‑1 (MQ=255)                  aaaaGCATAAATATGTTATTGACCTAATATTTAAATTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATg                                                                                                                 <  2:56680/148‑1 (MQ=255)                                     tGACCTAATATTTAAATTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCt                                                                                              <  1:7971/148‑1 (MQ=255)                                            atatTTAAATTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTg                                                                                       >  2:46747/1‑148 (MQ=255)                                                  aaaTTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGt                                                                                 <  1:3246/148‑1 (MQ=255)                                                  aaaTTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGt                                                                                 <  1:214833/148‑1 (MQ=255)                                                  aaaTTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGt                                                                                 <  2:352107/148‑1 (MQ=255)                                                  aaaTTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGt                                                                                 <  2:275045/148‑1 (MQ=255)                                                          catcatTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATAc                                                                         <  1:207728/148‑1 (MQ=255)                                                          catcatTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATAc                                                                         <  2:240635/148‑1 (MQ=255)                                                            tcatTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACtt                                                                       <  2:210928/148‑1 (MQ=255)                                                                                atatgatGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATAtt                                                   <  2:94254/148‑1 (MQ=255)                                                                                atatgatGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATAtt                                                   <  2:24857/148‑1 (MQ=255)                                                                                atatgatGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATAtt                                                   <  2:198470/148‑1 (MQ=255)                                                                                          aataatAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAattatta                                         >  2:18658/1‑148 (MQ=255)                                                                                          aataatAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAattatta                                         >  2:202124/1‑148 (MQ=255)                                                                                          aataatAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAattatta                                         >  2:155980/1‑148 (MQ=255)                                                                                          aataatAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAattatta                                         >  2:32118/1‑148 (MQ=255)                                                                                          aataatAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAattatta                                         >  1:94841/1‑148 (MQ=255)                                                                                          aataatAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAattatta                                         >  1:71767/1‑148 (MQ=255)                                                                                            taataaATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAattattatt                                       <  1:2535/148‑1 (MQ=255)                                                                                                aaaTTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCTAATTATTATTTgta                                   >  2:16800/1‑148 (MQ=255)                                                                                                aaaTTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAATTATTATTTgta                                   >  1:5931/1‑148 (MQ=255)                                                                                                                        tttACGTATTAAATAGATTGAGAATAAGAAATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAATTATTATTTGTAGTATTATAGTTAGGGTGGAAAAtt           <  2:137263/148‑1 (MQ=255)                                                                                                                                 tcgatcgcttGAGAGTCGGGTATATATAGTGCGTATGACTATGAAAGGATCGACAGTCTATTGAAGTGTTCTAATACTTCAGTAGTAAATCTATATATTTGATTTATATTTTGTAGTATGACGGTGAGGATGAAAAGGAATAAATTGt  >  2:329592/9‑148 (MQ=255)                                                                                                                                 tAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAATTATTATTTGTAGTATTATAGTTAGGGTGGAAAATTATAAATTGt  <  1:46747/148‑1 (MQ=255)                                                                                                                                 tAAATAGATTGAGAATAAGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAATTATTATTTGTAGTATTATAGTTAGGGTGGAAAATTATAAATTGt  >  2:19432/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                   TTTTGAAAGTATCCATAAAAAGCATAAATATGTTATTGACCTAATATTTAAATTGAACATCATTAAATTAAAAGTAATGATATGATGGAAATAATAAATTTAATAAGAAAATTTCATGTTTTACGTATTAAATAGATTGAGAATAGGGTATATATAATGCGTATGACTATGAAAGAAACGACCTTCTATTGAAATGTTCTAATACTTTAGTAGTAAATCTATATATTCCAATTATTATTTGTAGTATTATAGTTAGGGTGGAAAATTATAAATTGT  >  CP000730/646693‑646968

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