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%	150.4 / NA	48	*260* (TAG→TAA)	yjjV	putative DNase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/26);  total (22/26)

CAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCC  >  NC_000913/4614129‑4614382                                                                                                                                    |                                                                                                                           cAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGAt                                                                                                                     <  1:1325055/139‑1 (MQ=255)         gCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGtt                                                                                                             >  1:1192053/1‑139 (MQ=255)               ccGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCAcgcc                                                                                                       <  1:84035/139‑1 (MQ=255)                                          cGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCt                                                                            >  2:230415/1‑139 (MQ=255)                                               tGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGc                                                                       <  2:1353034/139‑1 (MQ=255)                                                         gTTGCGCCGGGAACCGGCGGATTAGATTGCGCAAGCGTTGTTTAATATCACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                             <  1:1377739/139‑1 (MQ=255)                                                         gTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                        <  2:276607/112‑1 (MQ=255)                                                         gTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                        >  1:276607/1‑112 (MQ=255)                                                         gTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                             <  1:230415/139‑1 (MQ=255)                                                           tGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAcgc                                                           >  1:1792735/1‑139 (MQ=255)                                                             cgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcgc                                                                                     <  1:1596244/111‑1 (MQ=255)                                                             cgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcgc                                                                                     >  2:1596244/1‑111 (MQ=255)                                                                   gAACCGGCGGATGAGATTGCGCGAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                   <  2:681431/139‑1 (MQ=255)                                                                     aCCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                 <  2:190866/139‑1 (MQ=255)                                                                       cggcggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                   <  2:927709/135‑1 (MQ=255)                                                                       cggcggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                   >  1:927709/1‑135 (MQ=255)                                                                              tGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAgg                                        >  2:1315305/1‑139 (MQ=255)                                                                                       gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATc                                                                                                  <  2:1357859/72‑1 (MQ=255)                                                                                       gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATc                                                                                                  <  2:864314/72‑1 (MQ=255)                                                                                       gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATc                                                                                                  >  1:1357859/1‑72 (MQ=255)                                                                                       gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATc                                                                                                  >  1:864314/1‑72 (MQ=255)                                                                                       gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                        >  1:620664/1‑82 (MQ=255)                                                                                       gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                        <  2:620664/82‑1 (MQ=255)                                                                                                gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                      >  1:1590285/1‑75 (MQ=255)                                                                                                gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                      <  2:1590285/75‑1 (MQ=255)                                                                                                gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTg                                                                           >  2:215906/1‑86 (MQ=255)                                                                                                gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCCTTCACGCCTCATCCGGCAGTTGGCGCACAATGCCTg                                                                           <  1:215906/86‑1 (MQ=21)                                                                                                 cTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                           <  2:760597/69‑1 (MQ=255)                                                                                                 cTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                           >  1:760597/1‑69 (MQ=255)                                                                                                    aataaCACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCgg                  <  2:1792735/139‑1 (MQ=255)                                                                                                              ataCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGc                                                                       >  1:1688391/1‑76 (MQ=255)                                                                                                              ataCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGc                                                                       <  2:1688391/76‑1 (MQ=255)                                                                                                                aCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGt                                                         >  1:1330744/1‑88 (MQ=255)                                                                                                                aCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGt                                                         <  2:1330744/88‑1 (MQ=255)                                                                                                                       ttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                      <  2:823199/52‑1 (MQ=255)                                                                                                                       ttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                      >  1:823199/1‑52 (MQ=255)                                                                                                                              gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                      <  2:444499/45‑1 (MQ=12)                                                                                                                              gCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                      >  1:444499/1‑45 (MQ=255)                                                                                                                               ccGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAcg                                                            >  2:1692048/1‑70 (MQ=255)                                                                                                                               ccGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAcg                                                            <  1:1692048/70‑1 (MQ=255)                                                                                                                               ccGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGc                     >  1:854030/1‑109 (MQ=255)                                                                                                                               ccGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGc                     <  2:854030/109‑1 (MQ=255)                                                                                                                                cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                       >  2:1663814/1‑90 (MQ=255)                                                                                                                                cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                       <  1:1663814/90‑1 (MQ=255)                                                                                                                                cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc  <  2:32237/127‑1 (MQ=255)                                                                                                                                cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc  >  1:32237/1‑127 (MQ=255)                                                                                                                                cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg    <  2:202522/125‑1 (MQ=255)                                                                                                                                cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg    >  1:202522/1‑125 (MQ=255)                                                                                                                                    |                                                                                                                           CAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCC  >  NC_000913/4614129‑4614382

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