Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	42,108	G→A	intergenic (‑177/‑295)	caiT ← / → fixA	L‑carnitine:gamma‑butyrobetaine antiporter/putative electron transfer flavoprotein FixA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	42,108	0	G	A	94.4%	43.7 / NA	18	intergenic (‑177/‑295)	caiT/fixA	L‑carnitine:gamma‑butyrobetaine antiporter/putative electron transfer flavoprotein FixA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (7/10);  minor base C (1/0);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGCGTGATATCTG  >  NC_000913/41881‑42349                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  cggAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTGAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATGTTTGTTTCGGGTGAATAGAAGGCGTTTTTTc                                                                                                                                                                                                                                        <  1:361698/239‑1 (MQ=255)                  gTTCTATTCCCGTTTTTCTCATTTCATTCTTCATGCGTTAATTTCACTGTAAAAACGAATATTTATTTTGCTCTACCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAAACTCGTGGACGCGTTAATGGCTAAATCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTCTCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGTTTAATAATCAGTTTGttt                                                                                                                                                                                                           >  1:158633/1‑250 (MQ=255)                          cccGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAAGAAGCGGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCgn                                                                                                                                                                                                                                      <  1:82362/216‑2 (MQ=255)                                                      gTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGACGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTaaaaa                                                                                                                                                                                       >  1:185860/1‑235 (MQ=255)                                                                    gAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCtn                                                                                                                                                                             <  1:110612/231‑2 (MQ=255)                                                                                                  tgtttTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGc                                                                                                                            <  1:18833/250‑1 (MQ=255)                                                                                                                        aaTATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCCATAAGATGTATTCTGCGATTGGTATCACATTTTTGATTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATCTTGGt                                                                                                     >  1:242003/1‑251 (MQ=255)                                                                                                                                               ggACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGNTGTNTTCTNTGNTTGGGATCACAGTTTNGTTTCGGGTAAATAGAAGGCTTTTTTTCGTTAATTTTGATTAAGNNTNNNTNTGTTATGc                                                                                                                                                                                                        <  2:312807/129‑1 (MQ=255)                                                                                                                                               ggACGCGCTACTGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGGTTGGTATCAACTTTTTGATTCGGGTACATCGAAGGCGTTTTTTCGTTAATTTTGCTTAATACTCAGTTTCTTATGc                                                                                                                                                                                                        >  1:312807/1‑129 (MQ=255)                                                                                                                                                  cgcgTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTGTCGTTAATTTAGATGAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCGGACTTTCAATATTGGTGAGCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTGTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTg                                                                           <  1:432481/251‑1 (MQ=255)                                                                                                                                                   gcgTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTa                                                                                                                                                                                                           >  1:188059/1‑122 (MQ=255)                                                                                                                                                     gTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACgn                                                                                                                       <  1:38521/204‑2 (MQ=255)                                                                                                                                                         aTGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTAGTATCACATTGTTGTTGCGGGTGAATAGAAGGCGTTTTGTCGGTAATTTTGATTAATAAGCAGGTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTg                                                                                                                             <  1:446933/194‑1 (MQ=255)                                                                                                                                                                  tCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGCATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGaaa                                                                               >  1:177845/1‑231 (MQ=255)                                                                                                                                                                                                  gATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTACCAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTAtt                           >  1:339322/1‑251 (MQ=255)                                                                                                                                                                                                                    gTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCt                                                                                                                                                                              <  1:8384/86‑1 (MQ=255)                                                                                                                                                                                                                    gTTTCGGGAGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCt                                                                                                                                                                              >  2:8384/1‑86 (MQ=255)                                                                                                                                                                                                                           gTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTGATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGCGTGATATCTg  <  1:102896/251‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  CGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGCGTGATATCTG  >  NC_000913/41881‑42349

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