Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,823,537	A→G	L78P (CTG→CCG)	recA ←	DNA recombination/repair protein RecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,823,537	0	A	G	100.0%	138.5 / NA	44	L78P (CTG→CCG)	recA	DNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/27);  total (17/27)

TCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTG  >  NC_000913/2823294‑2823780                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    tCGCCTTCGATTTCCGCTTTCGGCGCCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCGTCTCGCTGCGCTGCGGCGAGCACCTGcggcgtca                                                                                                                                                                                                                                               <  1:210252/250‑1 (MQ=255)              ccGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTAcc                                                                                                                                                                                                                                  <  2:5336/250‑1 (MQ=255)                   ttCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAg                                                                                                                                                                                                                             <  1:96175/250‑1 (MQ=255)                    tCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGa                                                                                                                                                                                                                            <  1:172310/250‑1 (MQ=255)                     cGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAgctt                                                                                                                                                                                                                          >  2:245469/1‑248 (MQ=255)                         gTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCCGCGTGGTTTTACCGGAAGATTCCgg                                                                                                                                                                                                                      >  2:223239/1‑251 (MQ=255)                           cAGTGCCGCCACGGAGTCAACGACGATAACTTCTACTGCGCCAGGACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTTCGTAGATTGGGTCCAGCGCGTGGTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTc                                                                                                                                                                                                                    <  1:46548/251‑1 (MQ=255)                           cAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTc                                                                                                                                                                                                                    <  1:217480/251‑1 (MQ=255)                             gTGCCGCCACGGAGTCAACGACGCTAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTCCCGGAAGATTCCGGTCCg                                                                                                                                                                                                                  >  2:42997/1‑251 (MQ=255)                              ggccgccACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGt                                                                                                                                                                                                                 <  2:46244/250‑1 (MQ=255)                                    cACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGAttt                                                                                                                                                                                                           <  2:260523/251‑1 (MQ=255)                                            cAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGAt                                                                                                                                                                                                    <  2:203103/250‑1 (MQ=255)                                             aaCGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATAc                                                                                                                                                                                                  <  1:234747/251‑1 (MQ=255)                                               cgacgaTAACGTCTACTGCGCCAGAACGCGCCAGGGCGGCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCGGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGGTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACgg                                                                                                                                                                                                <  2:144384/251‑1 (MQ=255)                                               cgacgaTAACGTCTACCGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCTACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCCCCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACgg                                                                                                                                                                                                <  1:152516/251‑1 (MQ=255)                                                          tCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCCGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGc                                                                                                                                                                                       >  1:262149/1‑249 (MQ=255)                                                           cTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCACGCCAGTGTCCGGCTGGGAGCACAGCAGCTTGTCGCTATCGACGCCCCGTTTACGTGCGTCGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGCTCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGc                                                                                                                                                                                       >  1:260529/1‑248 (MQ=255)                                                                gcgcCAGAACGCGCCAGTGCGGCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACTCTTCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCTACGATTCGGCCCCTCGGCAGaccac                                                                                                                                                                                <  2:142191/250‑1 (MQ=255)                                                                              cAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGcgcgg                                                                                                                                                                 >  2:260952/1‑250 (MQ=255)                                                                                    gTCACAGATTGCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTGCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTGTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCAGCGGCAGACCACCTGCTCCAAGCGCGATATCCa                                                                                                                                                           <  2:31248/251‑1 (MQ=255)                                                                                        caGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTcgac                                                                                                                                                                                                       <  1:104286/203‑1 (MQ=255)                                                                                        caGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTcgac                                                                                                                                                                                                       >  2:104286/1‑203 (MQ=255)                                                                                             ttCCAATGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGCTGTCGATATCGCCGCCCCGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCACGCGCGATCTCCAGTGAAAGCg                                                                                                                                                  >  2:139824/1‑251 (MQ=255)                                                                                                 aGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAAcc                                                                                                                                              <  1:76440/251‑1 (MQ=255)                                                                                                    gCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGt                                                                                                                                           <  1:245469/251‑1 (MQ=255)                                                                                                    gCCTGCGCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGt                                                                                                                                           <  1:223239/251‑1 (MQ=255)                                                                                                                  gtCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTcc                                                                                                                              <  1:65973/250‑1 (MQ=255)                                                                                                                           ggAGCACAGCAGGTTGTCGATATCGACGCCCATTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTGTGACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCAGCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATgg                                                                                                                    <  2:113083/251‑1 (MQ=255)                                                                                                                                                          aGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGcc                                                                                     <  1:260952/251‑1 (MQ=255)                                                                                                                                                                 gTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTAcc                                                                               >  1:18775/1‑250 (MQ=255)                                                                                                                                                                     gTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAAtt                                                                          <  1:125218/251‑1 (MQ=255)                                                                                                                                                                            gggTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTtctc                                                                   <  1:42997/251‑1 (MQ=255)                                                                                                                                                                             ggTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTGTCTCa                                                                  <  1:74989/251‑1 (MQ=255)                                                                                                                                                                                      cgTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGccc                                                         >  1:149621/1‑251 (MQ=255)                                                                                                                                                                                            cAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGtgc                                                    <  2:149621/250‑1 (MQ=255)                                                                                                                                                                                               cATCGATAAGCGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCCCCTTCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCTCTCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCGgcc                                                <  2:220667/251‑1 (MQ=255)                                                                                                                                                                                                  cGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTgccgcc                                             >  1:58874/1‑251 (MQ=255)                                                                                                                                                                                                       aaCGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGc                                        >  1:146456/1‑251 (MQ=255)                                                                                                                                                                                                          gCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCttt                                     >  1:40557/1‑251 (MQ=255)                                                                                                                                                                                                              aGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGACCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCtgt                                 >  1:142087/1‑251 (MQ=255)                                                                                                                                                                                                                  ttttccTTCACGCTGCGCTGCGGCGATCCCCTTCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCtgtttg                              <  2:146456/246‑1 (MQ=255)                                                                                                                                                                                                                               tgcgctgcgGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGCTTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTAGTCGATAGc                >  1:251022/1‑251 (MQ=255)                                                                                                                                                                                                                               tgcgctgcgGCGATCACCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaa                                                                                                                                                      >  2:106024/1‑117 (MQ=255)                                                                                                                                                                                                                               tgcgctgcgGCGATCACCGGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGaa                                                                                                                                                      <  1:106024/117‑1 (MQ=255)                                                                                                                                                                                                                                             aaCCTGCGGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCAGTAGATTTCGACGCTACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCCGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGATTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTCCTCCTg  >  2:249433/2‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    TCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGCGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTCCTG  >  NC_000913/2823294‑2823780

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