Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,105,994	T→C	T12T (ACT→ACC)	pdhD →	dihydrolipoyl dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,105,994	0	T	C	96.5%	102.6 / ‑6.3	29	T12T (ACT→ACC)	pdhD	dihydrolipoyl dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/13);  minor base G (0/1);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATTAAGCTTTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCATTACT  >  CP000730/1105848‑1106134                                                                                                                                                   |                                                                                                                                             tCATTAAGCTTTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACcat                                                                                                                                            <  1:42271/149‑1 (MQ=255)      aaGCTTTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTaa                                                                                                                                       >  1:22924/1‑149 (MQ=255)      aaGCTTTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGATTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGCAACAGATACCATAGTaa                                                                                                                                       >  1:225281/1‑149 (MQ=255)          tttGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCgg                                                                                                                                   >  1:134619/1‑149 (MQ=255)          tttGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCgg                                                                                                                                   >  2:156946/1‑149 (MQ=255)                 acaGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGt                                                                                                                            >  2:118299/1‑149 (MQ=255)                      cAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCtgg                                                                                                                       <  2:274423/149‑1 (MQ=255)                         aTTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCtggtgg                                                                                                                    <  2:193320/149‑1 (MQ=255)                             aTGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATAc                                                                                                                <  1:204741/149‑1 (MQ=255)                                                  cAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCt                                                                                           <  2:58243/149‑1 (MQ=255)                                                               tAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaa                                                                              >  2:31377/1‑149 (MQ=255)                                                                 aaCGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACaaaaa                                                                            >  1:164228/1‑149 (MQ=255)                                                                  tcGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAg                                                                           <  1:141157/148‑1 (MQ=255)                                                                    gTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTa                                                                         <  2:270923/149‑1 (MQ=255)                                                                            aataaaCCAGAAATATTATTAATGGAGAAAAAAAAAATAGTAGTTGGAGATTTCCCAATTGAAAAAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGAAAAAAAAAAAAAAAAAGTAACAATCGt                                                                 <  2:214452/149‑1 (MQ=255)                                                                             ataatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGtt                                                                >  1:260336/1‑149 (MQ=255)                                                                             ataatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGtt                                                                >  1:93212/1‑149 (MQ=255)                                                                             ataatCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGtt                                                                >  1:73935/1‑149 (MQ=255)                                                                                   cAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCACAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGAAAAAAAAAAACAAAAAGTAACAATCGTTGAGaaa                                                          <  2:94462/149‑1 (MQ=255)                                                                                            tattaATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATNAGGACAAAAAGTNACAATCGTTGAGAAAGGTAATCtt                                                 <  1:155891/149‑1 (MQ=255)                                                                                               taATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTtggt                                              >  2:77459/1‑149 (MQ=255)                                                                                                 aTGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTtggtgg                                            >  1:135175/1‑149 (MQ=255)                                                                                                             aaCGTGGTAGTTGGAGTTTTCGCAAGTGAAACAGATAGCATAGGAATTGGTGCAGGTCCTGGTGGATACGTTGCAGGGATTCGTGCGGCTCGATTAGTACAAAGAGTAACTGTCGTGGAGAAAGGTAATCTTGGTGGTGTTTGCTTaaa                                <  2:225281/149‑1 (MQ=255)                                                                                                             aaCATGGTAGTTGGAGATTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTaaa                                <  2:22924/149‑1 (MQ=255)                                                                                                                           agaTTTCCCAATTGAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTc                  >  1:111009/1‑149 (MQ=255)                                                                                                                              tttCCCAATTGAAACAGATACGATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGAAATTCGTGCAGATCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCtt               <  2:110724/149‑1 (MQ=255)                                                                                                                                        gAAACAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCAtt     >  1:33473/1‑149 (MQ=255)                                                                                                                                         aaaCAGATACCAGAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATGAGGACAAAAAGTTACAATCGTGTATAAATGTTATCGTGGTGGTGTTTGCGTAAACGTATGATGTAGTCCGTCAAAAGCATta    >  1:140029/1‑149 (MQ=255)                                                                                                                                           aCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCATtact  <  1:7239/149‑1 (MQ=255)                                                                                                                                           aCAGATACCATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCATtact  <  1:195816/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             TCATTAAGCTTTGACCACAGACAAATTGATGGTGCAACTGGCCAAAATGCAATGAATCACATTAAACGTTTATTAAATAATCCAGAATTATTATTAATGGAGGGGTAAAACATGGTAGTTGGAGATTTCCCAATTGAAACAGATACTATAGTAATCGGAGCAGGTCCTGGTGGATACGTTGCAGCAATTCGTGCAGCTCAATTAGGACAAAAAGTAACAATCGTTGAGAAAGGTAATCTTGGTGGTGTTTGCTTAAACGTAGGATGTATTCCTTCAAAAGCATTACT  >  CP000730/1105848‑1106134

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