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

GGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGT  >  NC_000913/4614132‑4614394                                                                                                                                 |                                                                                                                                       gggCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAg                                                                                                                              <  1:140307/139‑1 (MQ=255)    cAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCg                                                                                                                           >  1:1448362/1‑139 (MQ=255)       ccTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTc                                                                                                                        >  1:578786/1‑139 (MQ=255)                     gCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                          <  1:631682/139‑1 (MQ=255)                                 tgtgtTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGcaca                                                                                              <  2:1448362/139‑1 (MQ=255)                                                         gcgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCgg                                                                                                                                   <  1:59680/78‑1 (MQ=255)                                                         gcgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCgg                                                                                                                                   >  2:59680/1‑78 (MQ=255)                                                          cgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGt                                                                     <  1:978156/139‑1 (MQ=255)                                                               ggAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTAt                                                                <  2:1729875/139‑1 (MQ=255)                                                                      gcggATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCt                                                         >  1:1658442/1‑139 (MQ=255)                                                                                   cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGa                                                                                                                                  >  2:1131111/1‑53 (MQ=255)                                                                                   cgcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGa                                                                                                                                  <  1:1131111/53‑1 (MQ=255)                                                                                    gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGAt                                                                                                                                 <  2:616909/53‑1 (MQ=255)                                                                                    gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGAt                                                                                                                                 >  1:616909/1‑53 (MQ=255)                                                                                    gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAcgcg                                                                      <  2:1280342/112‑1 (MQ=255)                                                                                    gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAcgcg                                                                      >  1:1280342/1‑112 (MQ=255)                                                                                        gCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                                    <  2:1431294/78‑1 (MQ=255)                                                                                        gCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                                                                                    >  1:1431294/1‑78 (MQ=255)                                                                                        gCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAAcc                                       <  2:1658442/139‑1 (MQ=255)                                                                                           ttGCTTAATAAGATGGATACGTTGTTTAACCTGCCGTAAGCCGGTTGAGGGGTTCACGCCGCATCCGGCAGTTGGCGCACTATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTa                                    >  1:1847823/1‑139 (MQ=18)                                                                                             gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAAcc                                       <  2:465549/134‑1 (MQ=255)                                                                                             gCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAAcc                                       >  1:465549/1‑134 (MQ=255)                                                                                              cTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGc                                 <  2:578786/139‑1 (MQ=255)                                                                                                      cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                         <  2:1849972/91‑1 (MQ=255)                                                                                                      cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                                         >  1:1849972/1‑91 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgc                                                                                                   <  1:748042/63‑1 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgc                                                                                                   >  2:748042/1‑63 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCtt                                                                            >  1:987029/1‑86 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCtt                                                                            <  2:987029/86‑1 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                                   <  2:564567/111‑1 (MQ=255)                                                                                                        cGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                                   >  1:564567/1‑111 (MQ=255)                                                                                                         gTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                             >  1:225595/1‑100 (MQ=255)                                                                                                         gTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                             <  2:225595/100‑1 (MQ=255)                                                                                                           ataCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtt                    >  2:1207660/1‑139 (MQ=255)                                                                                                               gttgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg                <  1:1207660/139‑1 (MQ=255)                                                                                                                 tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                                  >  1:763261/1‑55 (MQ=255)                                                                                                                 tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                                  <  2:763261/55‑1 (MQ=255)                                                                                                                        cGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                         >  1:699498/1‑121 (MQ=255)                                                                                                                        cGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                         <  2:699498/121‑1 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGt                                                                     >  2:419225/1‑72 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGt                                                                     <  1:419225/72‑1 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                               >  2:1894984/1‑78 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                               <  2:1584488/78‑1 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                               <  1:1894984/78‑1 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                               >  1:1584488/1‑78 (MQ=255)                                                                                                                             cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt  >  2:1694287/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                       GGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGT  >  NC_000913/4614132‑4614394

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