Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	730,736	A→G	D385G (GAC→GGC)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	730,736	0	A	G	100.0%	111.5 / NA	31	D385G (GAC→GGC)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/5);  total (26/5)

TAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCA  >  NC_000913/730602‑730870                                                                                                                                        |                                                                                                                                       tAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGca                                                                                                                                   >  1:33439/1‑140 (MQ=255)  aaaTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcat                                                                                                                                  <  1:654625/140‑1 (MQ=255)      aCCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcatcacc                                                                                                                              >  1:491761/1‑140 (MQ=255)        cGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcatcaccat                                                                                                                            >  1:542725/1‑140 (MQ=255)        cGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAAAAGGGCCGcatcaccat                                                                                                                            >  1:227589/1‑140 (MQ=255)          gggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcatcaccatca                                                                                                                          >  1:596800/1‑140 (MQ=255)          gggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcatcaccat                                                                                                                            <  1:748387/138‑1 (MQ=255)          gggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcatcaccat                                                                                                                            >  2:748387/1‑138 (MQ=255)           ggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcatcaccat                                                                                                                            <  1:659862/137‑1 (MQ=255)           ggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcatcaccat                                                                                                                            >  2:659862/1‑137 (MQ=255)             ccGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCg                                                                                                                       <  1:1107919/140‑1 (MQ=255)                gATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcatca                                                                                                                                <  1:261179/128‑1 (MQ=255)                gATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGcatca                                                                                                                                >  2:261179/1‑128 (MQ=255)                     tggCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTg                                                                                                               >  1:980722/1‑140 (MQ=255)                       gcgcACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCTACAGCCTGaa                                                                                                             >  1:35754/1‑140 (MQ=255)                               cccccggggCCGACGGAAAA‑TTTGTTCCCGTTCCGACGGCGAGGGGGGGGGGACAGAACAGCAAACCCGGGCGGGTTTAAGTTACACGTACCAGTAGGAGAAGGGCCGCCCCCCCACCCCCGACAGCCTGAACCGCCGTg                                                                                                     <  2:1043151/134‑1 (MQ=255)                                   gggggCCG‑CCGGAAATTCTTTTACCGTTACGACAGCGAGGGGGGGGTGAAAAAACAGCAAACCCCGGCAGGTTTAATTTACCCGTATAAGTAGGAAAGGGCCCGCACCCCCACCCCCGACACCCTGAACCGCCGTGAAGt                                                                                                 <  2:542725/138‑1 (MQ=9)                                       gagagagggAAA‑TTTGTTCCCGTAACGACAGGGAGGGGGGGGGGAAAAAACAGCAAACCCCGGCGGTTTTAAGTTACACGTATCAGTATGAGAGGGCCCGCACCCCCATCCCCGACACCCTGAACCGCCGTGAAGTCCTg                                                                                             <  2:244104/133‑1 (MQ=9)                                       gCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTg                                                                                             >  1:1043151/1‑140 (MQ=255)                                       gCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTg                                                                                             >  2:194482/1‑140 (MQ=255)                                       gCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTg                                                                                             >  1:1212010/1‑140 (MQ=255)                                             ggaaaa‑TCTTTTCCCGTTACGACAGGGAGGGGGGGGTGAAAAAACAGCAAACCCCGGCGGTTTAAATTTCCACGTATCAGTATAAAAAGGCCCCCACCCCCACCCCCGACACCTTAAACCGCCGTGAAGTCCTGCACACg                                                                                       <  2:227589/137‑1 (MQ=255)                                                        ttccggTTACGACGGGGGGGGGGGGGGAACAAAACAGCAAACCCCGGCAGGCTTAATTTACACGTATCAGTAGGAGAGGGGCCGCACCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAg                                                                             <  2:360091/135‑1 (MQ=255)                                                        ttACCGTTACGACAGCGAGGGGGGGGTGAAAGAACAGCTAACCCCGGCAGGTTTAAGTTACACGTATCAGTATGAGAGGGGCCGCATCCCCATCCCCTACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAg                                                                             <  2:35754/140‑1 (MQ=14)                                                          aCCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAgggctgtctc                                                                           >  1:244104/1‑133 (MQ=255)                                                             gTTACGACGGCGAGGGGGGGGTGAAAAACCAGCTAACCCGGGCGGGTTAAAGTTACACATATCAGTATGAGAGGGGCCCCACCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCggg                                                                        <  2:421637/140‑1 (MQ=255)                                                             gTTACGACAGCGAGGGGGGGGTGAAAAAACACCCAACCCCGGCAGACTTAAGTTCCTCGTATCATTAGGAGAGGGGCCCCCCCCCCACCCCCGACGCCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCggg                                                                        <  2:160505/140‑1 (MQ=255)                                                                aCGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTg                                                                     >  1:1070444/1‑140 (MQ=255)                                                                  gACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGaa                                                                   >  1:1104551/1‑140 (MQ=255)                                                                       cGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCggg                                                              >  1:160505/1‑140 (MQ=255)                                                                        ggggggggggTGACAGAACAGCAAACCCCGGCGGGTTTAAGTTACAGGTATCAGTATGAGAGGGGCCCCACCCCCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGggt                                                             <  2:180193/137‑1 (MQ=255)                                                                            ggcgggTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGgtggtg                                                         >  1:810065/1‑140 (MQ=255)                                                                                             cTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAacactgtct                                        >  1:36900/1‑135 (MQ=255)                                                                                                accccGGCGGGTTAAAGTAACACGTATAAGTAGGAGAGGGGCCCCACCCCCACCCCCGACACCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACgg                                     <  2:1212010/138‑1 (MQ=14)                                                                                                aaCCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACgg                                     >  1:1321124/1‑140 (MQ=255)                                                                                                     ggCAGGCTTAAGTTACACGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCg                                >  1:782639/1‑140 (MQ=255)                                                                                                          gCTAAAGTTACACGTATCAGTAGGAGAAGGCCCGCACCCCCACCCCCGACAGCCTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACg                           <  2:805520/140‑1 (MQ=14)                                                                                                                    cacGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGt                 >  1:144602/1‑140 (MQ=255)                                                                                                                    cacGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGt                 >  1:1122123/1‑140 (MQ=255)                                                                                                                      cGTATCAGTATGAGAAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGttt               >  1:547314/1‑140 (MQ=255)                                                                                                                            gaaaggagaAGGCCCGCACCACCACCCCCAACAGCAGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGAcgcg         <  2:590847/135‑1 (MQ=14)                                                                                                                                 gagaAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTggg    >  2:1116660/1‑140 (MQ=255)                                                                                                                                 gagaAGGCCCGCACCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTggg    <  2:1104551/140‑1 (MQ=14)                                                                                                                                   gaAGGGCCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCa  >  1:571727/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       TAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAA‑TCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCA  >  NC_000913/730602‑730870

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: ‑