Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,049,827	T→G	intergenic (‑83/‑6)	gfcA ← / → insA‑4	threonine‑rich inner membrane protein GfcA/IS1 protein InsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,049,827	0	T	G	100.0%	151.2 / NA	43	intergenic (‑83/‑6)	gfcA/insA‑4	threonine‑rich inner membrane protein GfcA/IS1 protein InsA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/26);  total (17/26)

CTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAG  >  NC_000913/1049693‑1049900                                                                                                                                       |                                                                          cTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCa                                                                      <  2:798945/140‑1 (MQ=255)  tGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAg                                                                     >  2:373314/1‑140 (MQ=255)   ggCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGt                                                                    <  1:1101661/140‑1 (MQ=255)   ggCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGt                                                                    <  2:1190473/140‑1 (MQ=255)     cgtcgtCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTgg                                                                  >  2:541313/1‑140 (MQ=255)      gtcgtcAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGc                                                                 <  2:943661/140‑1 (MQ=255)       tcgtcAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                                                                >  1:117285/1‑140 (MQ=255)        cgtcAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                               <  1:1038397/140‑1 (MQ=255)        cgtcAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                               <  1:541313/140‑1 (MQ=255)        cgtcAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGATTTTGAGGTGCTCCAGTGGCtt                                                               <  2:1086255/140‑1 (MQ=255)                tGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCt                                                       <  2:347572/140‑1 (MQ=255)                  aaGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGTCAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTAt                                                     <  2:771168/140‑1 (MQ=255)                    ggCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa                                                   <  1:373314/140‑1 (MQ=255)                     gCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAg                                                  <  1:358708/140‑1 (MQ=255)                      ccATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGc                                                 <  1:162901/140‑1 (MQ=255)                       cATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTAAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCt                                                >  2:843977/1‑140 (MQ=255)                        aTAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                                            >  1:761367/1‑127 (MQ=255)                        aTAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                                            <  2:761367/127‑1 (MQ=255)                        aTAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTg                                               >  1:20724/1‑140 (MQ=255)                         tAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                                              >  1:1306612/1‑140 (MQ=255)                             gATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcct                                          <  1:2347/140‑1 (MQ=255)                              aTTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctc                                         <  2:1306612/140‑1 (MQ=255)                              aTTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctc                                         <  2:1204004/140‑1 (MQ=255)                              aTTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctc                                         <  2:444712/140‑1 (MQ=255)                               ttGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcc                                        >  1:998409/1‑140 (MQ=255)                               ttGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcc                                        >  1:1286967/1‑140 (MQ=255)                               ttGCAGAAAGCTTGTGGTTCATAACTTTTCCTTTATTCATCGCATGGACAATCCGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTCTTGAGGTGCTCCAGTGGCTACTGTTTCTATCAGCTGTCcctcc                                        >  2:620367/1‑140 (MQ=21)                                tGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                                              <  2:504886/133‑1 (MQ=255)                                tGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                                              >  1:504886/1‑133 (MQ=255)                                     aaaGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCa                                  >  1:1117890/1‑140 (MQ=255)                                       aGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGc                                >  1:886256/1‑140 (MQ=255)                                       aGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGc                                <  2:1286967/140‑1 (MQ=255)                                        gCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCt                               <  2:1199301/140‑1 (MQ=255)                                             tgtTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGa                          <  1:570184/140‑1 (MQ=255)                                               tttCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg                        <  1:203784/140‑1 (MQ=255)                                                  cATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgggg                     <  2:328883/140‑1 (MQ=255)                                                   aTAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTATTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggt                    >  2:865461/1‑140 (MQ=255)                                                          ttCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGATGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGt             >  2:932707/1‑140 (MQ=255)                                                          ttCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGt             >  2:864041/1‑140 (MQ=255)                                                             cTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAAc          <  1:367883/140‑1 (MQ=255)                                                             cTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTAAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAAc          <  1:843977/140‑1 (MQ=255)                                                               ttATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACgg        >  1:969997/1‑140 (MQ=255)                                                                  ttCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaa     >  1:675679/1‑140 (MQ=255)                                                                     aTCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAg  <  2:1117890/140‑1 (MQ=255)                                                                                                                                       |                                                                          CTGGCGTCGTCAACATGAAGGCCATAAGGATTGCAGAAAGCTTGTGTTTCATAACTTTTCCTTTATTCATCGCATGGACAATACGGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAG  >  NC_000913/1049693‑1049900

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