Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,614,260	G→A	*260* (TAG→TAA)	yjjV →	putative DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,614,260	0	G	A	100.0%	191.0 / NA	60	*260* (TAG→TAA)	yjjV	putative DNase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/24);  total (36/24)

CAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC  >  NC_000913/4614135‑4614391                                                                                                                              |                                                                                                                                    cAGCCTAACCGCCGGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGAGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCg                                                                                                                        >  2:868303/1‑139 (MQ=255)   gCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGtt                                                                                                                      >  2:534430/1‑139 (MQ=255)   gCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGtt                                                                                                                      >  1:714067/1‑139 (MQ=255)              cGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATc                                                                                                           >  2:284905/1‑139 (MQ=255)                        tGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                                 <  1:158032/139‑1 (MQ=255)                        tGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                                 <  2:847974/139‑1 (MQ=255)                          cccGTGTGTTTGCCGTGCCTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                                 >  1:1047606/1‑137 (MQ=255)                          cccGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                                 <  2:1047606/137‑1 (MQ=255)                          cccGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                               >  1:30812/1‑139 (MQ=255)                           ccGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcgc                                                                                              >  1:142805/1‑139 (MQ=255)                           ccGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcgc                                                                                              >  1:489144/1‑139 (MQ=255)                                    cGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCt                                                                                     >  1:10826/1‑139 (MQ=255)                                        gTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATg                                                                                 <  2:489144/139‑1 (MQ=255)                                            tttGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGAc                                                                             >  1:808429/1‑139 (MQ=255)                                                   gTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                      <  2:142805/139‑1 (MQ=255)                                                             gAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                            <  2:57819/139‑1 (MQ=255)                                                               aCCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                          <  2:30812/139‑1 (MQ=255)                                                                 cggcggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGc                                                        >  2:896496/1‑139 (MQ=255)                                                                                cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                          <  2:995238/122‑1 (MQ=255)                                                                                cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                          >  1:995238/1‑122 (MQ=255)                                                                                     gCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                      >  1:519182/1‑105 (MQ=255)                                                                                     gCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                      <  2:519182/105‑1 (MQ=255)                                                                                           cTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                          <  2:846324/111‑1 (MQ=255)                                                                                           cTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                          >  1:846324/1‑111 (MQ=255)                                                                                               ataaCACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGa                          >  1:369724/1‑139 (MQ=255)                                                                                                taaCACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGCTTGTGCCGAACCGTACACCGGaa                         >  2:341725/1‑138 (MQ=255)                                                                                                   cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAAcc                                    >  1:569153/1‑125 (MQ=255)                                                                                                   cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAAcc                                    <  2:569153/125‑1 (MQ=255)                                                                                                   cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                      >  2:111006/1‑139 (MQ=255)                                                                                                     cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCAcgccgc                                                                                                              >  1:1251533/1‑49 (MQ=255)                                                                                                     cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCAcgccgc                                                                                                              <  2:1251533/49‑1 (MQ=255)                                                                                                     cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                                    >  1:83327/1‑59 (MQ=255)                                                                                                     cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                                    <  2:83327/59‑1 (MQ=255)                                                                                                     cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTAt                                                             >  2:1063726/1‑98 (MQ=255)                                                                                                     cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTAt                                                             <  1:1063726/98‑1 (MQ=255)                                                                                                        ataCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtt                 >  1:632530/1‑139 (MQ=255)                                                                                                        ataCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtt                 >  1:1244394/1‑139 (MQ=255)                                                                                                              tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAg                                                                                                     <  2:1007889/49‑1 (MQ=255)                                                                                                              tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAg                                                                                                     >  1:1007889/1‑49 (MQ=255)                                                                                                              tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATg                                                                                 <  1:531478/69‑1 (MQ=255)                                                                                                              tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATg                                                                                 >  2:531478/1‑69 (MQ=255)                                                                                                              tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc           <  1:284905/139‑1 (MQ=255)                                                                                                              tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc           >  2:557191/1‑139 (MQ=255)                                                                                                              tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc           <  1:557191/139‑1 (MQ=255)                                                                                                              tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc           >  2:51643/1‑139 (MQ=255)                                                                                                                tttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                                >  2:1187078/1‑100 (MQ=255)                                                                                                                tttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                                <  1:1187078/100‑1 (MQ=255)                                                                                                                tttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgccgc         >  2:310836/1‑139 (MQ=255)                                                                                                                tttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgccgc         >  1:125464/1‑139 (MQ=255)                                                                                                                  tAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCAt       >  1:245120/1‑139 (MQ=255)                                                                                                                     cGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg    >  1:1043419/1‑139 (MQ=255)                                                                                                                       tGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc  <  2:125464/139‑1 (MQ=255)                                                                                                                        gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                      <  1:88668/70‑1 (MQ=255)                                                                                                                        gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                      >  2:88668/1‑70 (MQ=255)                                                                                                                        gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                                        >  2:431202/1‑100 (MQ=255)                                                                                                                        gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                                        <  1:431202/100‑1 (MQ=255)                                                                                                                        gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGc                              >  2:316638/1‑110 (MQ=255)                                                                                                                        gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGc                              <  1:316638/110‑1 (MQ=255)                                                                                                                          cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                      >  1:74253/1‑68 (MQ=255)                                                                                                                          cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                      <  2:74253/68‑1 (MQ=25)                                                                                                                              |                                                                                                                                    CAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGC  >  NC_000913/4614135‑4614391

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