Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,695,724	G→A	A116A (GCC→GCT)	uidA ←	beta‑glucuronidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,695,724	0	G	A	100.0%	128.3 / NA	38	A116A (GCC→GCT)	uidA	beta‑glucuronidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/21);  total (17/21)

TAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCTGCCCAACCTTTCGGTATAAAGACTTC  >  NC_000913/1695593‑1695858                                                                                                                                    |                                                                                                                                       tAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCaaa                                                                                                                                <  2:42789/140‑1 (MQ=255) tAAGACTGCTTTTTCATGCCGTTTTGGTCCGTAATCAACATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCaaa                                                                                                                                <  1:229645/140‑1 (MQ=255)      cTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCg                                                                                                                           <  1:430746/140‑1 (MQ=255)                   ccGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTg                                                                                                              >  2:138237/1‑140 (MQ=255)                        ttCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACACACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCt                                                                                                         <  2:421710/140‑1 (MQ=255)                            tcgGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCAt                                                                                                     <  1:389624/140‑1 (MQ=255)                            tcgGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCAt                                                                                                     <  1:88841/140‑1 (MQ=255)                               gTAATCACAATTCCCGATGGGAAAGTCTGCCTGCTCAGTCCGTTGTTCTCACAAACGGTGATAAGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCAc                                                                                                  <  2:329622/140‑1 (MQ=255)                                         ttCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGatta                                                                                        <  2:488430/140‑1 (MQ=255)                                          tCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGattat                                                                                       >  2:588952/1‑140 (MQ=255)                                              ggcggGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTGAAATGGCGGATACCCGACCTGATGCTCCATAACTTCCTGATTATTGAc                                                                                   >  2:227115/1‑140 (MQ=255)                                               gcggGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGatta                                                                                        >  1:28844/1‑134 (MQ=255)                                               gcggGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGatta                                                                                        <  2:28844/134‑1 (MQ=255)                                                       gTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACttt                                                                          <  2:552273/140‑1 (MQ=255)                                                             cagttcagttCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTa                                                                    >  1:701383/1‑140 (MQ=255)                                                                ttcagttcGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGAc                                                                                   <  2:490406/122‑1 (MQ=255)                                                                ttcagttcGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGAc                                                                                   >  1:490406/1‑122 (MQ=255)                                                                   agttcGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGt                                                              >  1:611315/1‑140 (MQ=255)                                                                   agttcGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGt                                                              <  2:36834/140‑1 (MQ=255)                                                                        gttgttCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCg                                                         <  2:382018/140‑1 (MQ=255)                                                                          tgttCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCa                                                       >  2:444601/1‑140 (MQ=255)                                                                                  caAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACg                                               <  2:611315/140‑1 (MQ=255)                                                                                    aaCGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACTCCTTGCCGTAATGAGAGACCGCATCGAAACgca                                             >  2:315956/1‑140 (MQ=255)                                                                                    aaCGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgca                                             <  2:660349/140‑1 (MQ=255)                                                                                      cGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACgcagc                                           >  2:610927/1‑140 (MQ=255)                                                                                        gTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCAc                                         >  1:601431/1‑140 (MQ=255)                                                                                                tACACTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCt                                 >  1:36021/1‑140 (MQ=255)                                                                                                   acTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGc                              >  1:489340/1‑140 (MQ=255)                                                                                                   acTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGc                              >  1:412389/1‑140 (MQ=255)                                                                                                   acTTTTCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGc                              <  2:601431/140‑1 (MQ=255)                                                                                                     ttttCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCt                            >  2:116643/1‑140 (MQ=255)                                                                                                        tCCCGGCAATAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCTGcc                         <  2:539852/140‑1 (MQ=255)                                                                                                                 tAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCTGCCCAACCTTTc                <  1:737173/140‑1 (MQ=255)                                                                                                                 tAACATACGGCGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCTGCCCAACCTTTc                <  2:678620/140‑1 (MQ=255)                                                                                                                 tAACATACGGCGTGAAATCAGCTTCCAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCTGCCCAACCTTTc                <  1:315956/140‑1 (MQ=255)                                                                                                                           cGTGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCTGCCCAACCTTTCGGTATAAAGa      >  2:427724/1‑140 (MQ=255)                                                                                                                             tGACATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCTGCCCAACCTTTCGGTATAAAGACt    >  2:583507/1‑140 (MQ=255)                                                                                                                               aCATCAGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCTGCCCAACCTTTCGGTATAAAGACTTc  <  2:36021/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       TAAGACTGCTTTTTCTTGCCGTTTTCGTCGGTAATCACCATTCCCGGCGGGATAGTCTGCCAGTTCAGTTCGTTGTTCACACAAACGGTGATACGTACACTTTTCCCGGCAATAACATACGGCGTGACATCGGCTTCAAATGGCGTATAGCCGCCCTGATGCTCCATCACTTCCTGATTATTGACCCACACTTTGCCGTAATGAGTGACCGCATCGAAACGCAGCACGATACGCTGGCCTGCCCAACCTTTCGGTATAAAGACTTC  >  NC_000913/1695593‑1695858

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