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

GTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCC  >  NC_000913/4614124‑4614382                                                                                                                                         |                                                                                                                           gTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGc                                                                                                                          >  1:769474/1‑139 (MQ=255)      cAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGAt                                                                                                                     <  2:770611/139‑1 (MQ=255)                   aCCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCAcgc                                                                                                        >  1:1097596/1‑139 (MQ=255)                                 gCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGtt                                                                                          <  2:1113834/139‑1 (MQ=255)                                     cccGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                      <  2:1097596/139‑1 (MQ=255)                                      ccGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcgc                                                                                     >  2:781054/1‑139 (MQ=255)                                        gtgtgtTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGcac                                                                                   <  2:769474/139‑1 (MQ=255)                                        gcgtgtTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGcac                                                                                   <  2:607942/137‑1 (MQ=255)                                                              gTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAAc                                                             <  1:1415404/139‑1 (MQ=255)                                                               ttGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCAcg                                                                                                         >  1:1115591/1‑94 (MQ=255)                                                               ttGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCAcg                                                                                                         <  2:1115591/94‑1 (MQ=255)                                                                 gcgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATgtg                                                                      <  1:395820/127‑3 (MQ=255)                                                                 gcgcCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATgtg                                                                      >  2:395820/1‑125 (MQ=255)                                                                                    gagaTTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                      <  2:748030/92‑1 (MQ=255)                                                                                    gagaTTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGgcg                                                                                      >  1:748030/1‑92 (MQ=255)                                                                                        ttGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGc                                                                                                                          >  1:743982/1‑52 (MQ=255)                                                                                            gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTaa                                                              <  1:386472/108‑1 (MQ=255)                                                                                            gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTaa                                                              >  2:386472/1‑108 (MQ=255)                                                                                            gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                       >  2:1121765/1‑131 (MQ=255)                                                                                            gcAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGt                                       <  1:1121765/131‑1 (MQ=255)                                                                                                gCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAAcc                           <  2:603721/139‑1 (MQ=255)                                                                                                 cGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCAcgccg                                                                                                      <  2:432777/63‑1 (MQ=255)                                                                                                 cGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCAcgccg                                                                                                      >  1:432777/1‑63 (MQ=255)                                                                                                          ataaCACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                               >  1:991296/1‑125 (MQ=255)                                                                                                          ataaCACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCg                               <  2:991296/125‑1 (MQ=255)                                                                                                           taaCACGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGAt                >  1:1195222/1‑139 (MQ=255)                                                                                                              cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                   <  2:1486780/101‑1 (MQ=255)                                                                                                              cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                   >  1:1486780/1‑101 (MQ=255)                                                                                                              cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg             >  2:577271/1‑139 (MQ=255)                                                                                                              cacGTATACGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg             >  1:982571/1‑139 (MQ=255)                                                                                                                     aCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGcaca                                                                                  <  1:246484/63‑1 (MQ=255)                                                                                                                     aCGTTGTTTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGcaca                                                                                  >  2:246484/1‑63 (MQ=255)                                                                                                                       gttgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg    <  2:1195222/139‑1 (MQ=255)                                                                                                                         tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                           <  1:358556/50‑1 (MQ=255)                                                                                                                         tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGt                                                                                           >  2:358556/1‑50 (MQ=255)                                                                                                                         tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTaa                                                              >  2:1358634/1‑79 (MQ=255)                                                                                                                         tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTaa                                                              <  1:1358634/79‑1 (MQ=255)                                                                                                                         tgttTAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc  <  2:982571/139‑1 (MQ=255)                                                                                                                           tttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGca                                                                                    <  2:351451/55‑1 (MQ=255)                                                                                                                           tttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGca                                                                                    >  1:351451/1‑55 (MQ=255)                                                                                                                           tttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg    <  2:1237097/135‑1 (MQ=255)                                                                                                                           tttAACGTGCCGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg    >  1:1237097/1‑135 (MQ=255)                                                                                                                                     cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                   >  2:1254892/1‑78 (MQ=255)                                                                                                                                     cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                   <  1:1254892/78‑1 (MQ=255)                                                                                                                                     cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                   <  2:729175/78‑1 (MQ=255)                                                                                                                                     cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                   >  1:729175/1‑78 (MQ=255)                                                                                                                                     cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGc                     >  2:646582/1‑108 (MQ=255)                                                                                                                                     cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGc                     <  1:646582/108‑1 (MQ=255)                                                                                                                                     cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc  >  2:1480625/1‑127 (MQ=255)                                                                                                                                     cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgcc  <  1:1480625/127‑1 (MQ=255)                                                                                                                                         |                                                                                                                           GTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACACGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCC  >  NC_000913/4614124‑4614382

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