Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,532,727	T→C	H304H (CAT→CAC)	ydcC →	H repeat‑associated putative transposase YdcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,532,727	0	T	C	100.0%	99.2 / NA	35	H304H (CAT→CAC)	ydcC	H repeat‑associated putative transposase YdcC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/16);  total (19/16)

AATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAA  >  NC_000913/1532591‑1532784                                                                                                                                         |                                                          aaTTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAcgtg                                                        <  2:646338/140‑1 (MQ=37)  aTTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTgg                                                       >  2:1060593/1‑140 (MQ=38)   ttATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAcgt                                                         >  1:1383571/1‑137 (MQ=37)     aTGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGgaga                                                    <  1:1075216/140‑1 (MQ=37)      tGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGaa                                                   <  1:776375/140‑1 (MQ=37)        cGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATAGCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATa                                                 >  2:1057231/1‑140 (MQ=255)                 ctcCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCAca                                        <  1:453819/140‑2 (MQ=37)                 ctcCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCAca                                        <  2:868901/140‑2 (MQ=37)                 ctcCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCAca                                        >  2:454189/1‑139 (MQ=255)                    cTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                                     <  2:1193503/140‑1 (MQ=255)                     tttCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCg                                    >  2:1273005/1‑140 (MQ=255)                          gTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTgg                               <  2:1243852/140‑1 (MQ=255)                            ccATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGAc                             >  1:131109/1‑140 (MQ=255)                             cATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGAc                             >  1:1078223/1‑139 (MQ=255)                             cATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGAc                             <  2:1078223/139‑1 (MQ=255)                             cATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                            <  1:1226186/140‑1 (MQ=255)                              ataataGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                           >  2:414221/1‑140 (MQ=255)                                  taGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTa                       <  2:838171/140‑1 (MQ=16)                                  taGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTa                       <  2:1339058/140‑1 (MQ=21)                                     cAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatg                    >  1:908107/1‑140 (MQ=25)                                     cAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatg                    >  2:973314/1‑140 (MQ=255)                                         aCAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaatg                <  1:454189/140‑1 (MQ=16)                                          cAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAatgaatga               >  1:913450/1‑140 (MQ=255)                                          cAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAatgaatga               >  2:761426/1‑140 (MQ=255)                                              agaaagaGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac           <  1:1273005/140‑1 (MQ=16)                                              agaaagaGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac           <  1:1060593/140‑1 (MQ=16)                                                  agagCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACt       >  2:496773/1‑140 (MQ=25)                                                    agCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACtac     >  1:947400/1‑138 (MQ=16)                                                     gCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAGGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa    >  1:581659/1‑140 (MQ=255)                                                     gCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACtaca    <  1:414221/140‑4 (MQ=255)                                                     gCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa    >  1:1067293/1‑140 (MQ=16)                                                     gCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa    >  1:706981/1‑140 (MQ=16)                                                       tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac           <  1:993173/131‑1 (MQ=2)                                                       tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac           >  2:993173/1‑131 (MQ=17)                                                       tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACta      >  2:752146/1‑135 (MQ=12)                                                       tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACta      <  1:752146/136‑2 (MQ=2)                                                       tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaa  >  1:245375/1‑140 (MQ=12)                                                                                                                                         |                                                          AATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAA  >  NC_000913/1532591‑1532784

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