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	100.0%	61.9 / NA	19	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);  new base A (6/13);  total (6/13)

CGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGCGTGATATCTGTA  >  NC_000913/41868‑42351                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    cgaTTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGtttttt                                                                                                                                                                                                                                           <  1:26065/251‑1 (MQ=255)      attaACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTa                                                                                                                                                                                                                                      <  1:26833/251‑1 (MQ=255)                     aaaaCCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTaataat                                                                                                                                                                                                                       >  2:40064/1‑251 (MQ=255)                            tCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGttt                                                                                                                                                                                                                 <  1:233382/250‑1 (MQ=255)                                           tttttCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTtgtg                                                                                                                                                                                                 >  1:207716/1‑251 (MQ=255)                                                    tttCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTaaaaa                                                                                                                                                                                         >  2:42540/1‑250 (MQ=255)                                                                  aGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAAAAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACttt                                                                                                                                                                          <  2:207716/251‑1 (MQ=255)                                                                                  aaaaCGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATcc                                                                                                                                                           >  1:147537/1‑250 (MQ=255)                                                                                                 tttGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGa                                                                                                                                            <  1:130119/250‑1 (MQ=255)                                                                                                           cGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCtttt                                                                                                                                 >  1:138484/1‑251 (MQ=255)                                                                                                              ttgtttTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTg                                                                                                                               <  2:262085/250‑1 (MQ=255)                                                                                                                         ttGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTTCGCCgtgt                                                                                                                    <  2:179586/250‑1 (MQ=255)                                                                                                                                    tAATATAATTATTATTAACCTCGTTGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTg                                                                                                         <  2:252999/250‑1 (MQ=255)                                                                                                                                                            ggACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGaa                                                                                  >  2:222551/1‑249 (MQ=255)                                                                                                                                                                            aaCTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCgt                                                                <  1:96320/251‑1 (MQ=255)                                                                                                                                                                                    aTGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATGCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACGATATTTAAAGTTGGATTTATCTGCGTGTGACAtt                                                        <  2:187535/251‑1 (MQ=255)                                                                                                                                                                                    aTGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACAtt                                                        <  1:252224/251‑1 (MQ=255)                                                                                                                                                                                                                        cacaTTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTa                    <  2:67810/251‑1 (MQ=255)                                                                                                                                                                                                                                          gAATAGAAGGCGTTTTGTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGCGTGATATCTGTa  <  2:268336/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    CGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGCGTGATATCTGTA  >  NC_000913/41868‑42351

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