Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	317,653	A→G	intergenic (‑517/+73)	ykgA ← / ← rclC	putative DNA‑binding transcriptional regulator YkgA/reactive chlorine species resistance protein C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	317,653	0	A	G	100.0%	122.2 / NA	35	intergenic (‑517/+73)	ykgA/rclC	putative DNA‑binding transcriptional regulator YkgA/reactive chlorine species resistance protein C
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/20);  total (15/20)

TATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAA  >  NC_000913/317579‑317729                                                                           |                                                                             tATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCg                                                                         <  2:66113/80‑1 (MQ=255)   ttttACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTc                                                                       <  1:66114/80‑1 (MQ=255)   ttttACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTc                                                                       <  2:66115/80‑1 (MQ=255)     ttACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAgg                                                                    >  2:66148/1‑81 (MQ=255)      tACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGc                                                                   >  1:66159/1‑81 (MQ=255)       aCCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCt                                                                  >  1:66149/1‑81 (MQ=255)            aGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTtcat                                                              <  2:66116/80‑1 (MQ=255)             gacgacAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTtcatca                                                            >  2:66157/1‑81 (MQ=255)             gacgacAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTtcatca                                                            >  1:66150/1‑81 (MQ=255)             gacgacAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCACATCTCCGGTAATGTGGCGTCAGGCTtcatca                                                            >  2:66162/1‑81 (MQ=255)                gacAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTtcatcat                                                           >  2:66156/1‑79 (MQ=255)                   aaaaaCCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGaa                                                      <  1:66117/81‑1 (MQ=255)                     aaaCCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAg                                                    <  1:66118/81‑1 (MQ=255)                        ccTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAgggg                                                 >  1:66158/1‑81 (MQ=255)                            ttAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGcc                                              <  2:66119/80‑1 (MQ=255)                             tAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCt                                             >  1:66163/1‑80 (MQ=255)                              aaCAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTg                                            <  2:66120/80‑1 (MQ=255)                                   ggAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCt                                      <  2:66121/81‑1 (MQ=255)                                    gAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTg                                     <  2:66122/81‑1 (MQ=255)                                     aGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTg                                     >  1:66161/1‑80 (MQ=255)                                         atatTGATCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGaaaa                                 <  1:66123/80‑1 (MQ=255)                                                tCTTCACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaac                          >  1:66160/1‑80 (MQ=255)                                                    cACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaacaac                       <  1:66126/79‑1 (MQ=255)                                                    cACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaacaa                        >  2:66165/1‑78 (MQ=255)                                                    cACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaacaa                        <  1:66125/78‑1 (MQ=255)                                                    cACTCTGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaac                          <  2:66124/76‑1 (MQ=255)                                                      ctctGTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagac                   <  2:66127/81‑1 (MQ=255)                                                          gTCATATCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAgtg               >  2:66164/1‑81 (MQ=255)                                                             atatCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGtt             <  1:66128/80‑1 (MQ=255)                                                             atatCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTc            <  1:66129/81‑1 (MQ=255)                                                              tatCTCCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCg           <  2:66130/81‑1 (MQ=255)                                                                 ctcCGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAATTAGTGAACAACAGACAGTGTTCGGAt        >  2:66172/1‑81 (MQ=255)                                                                    cGGTAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGAtt       <  1:66131/79‑1 (MQ=255)                                                                       tAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCaa  >  2:66154/1‑81 (MQ=255)                                                                       tAATGTGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCaa  <  1:66132/81‑1 (MQ=255)                                                                           |                                                                             TATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAA  >  NC_000913/317579‑317729

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