Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,018,586	G→A	L418L (CTG→CTA)	rmuC →	putative recombination limiting protein RmuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,018,586	0	G	A	100.0%	139.0 / NA	41	L418L (CTG→CTA)	rmuC	putative recombination limiting protein RmuC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/25);  total (16/25)

ACAACATCATGCTGGTTAGCCCGACTACGCTGCTGGTGGCGCTGCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTGGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGA  >  NC_000913/4018345‑4018827                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   acaacaTCATGCTGGTTAGCCCGACTACGCTGCTGGTGGCGCTTCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCgcaggc                                                                                                                                                                                                                                            >  1:220270/1‑250 (MQ=255)           gCTGGTTAGCCCGACTACGCTGCTGGTGGCGCTGCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTc                                                                                                                                                                                                                                 <  1:217610/251‑1 (MQ=255)                   gccagtcTACGCTGCTGGTGGCGCTGCGCTCTATCGCCAACCTGGGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTa                                                                                                                                                                                                                         <  2:54644/245‑1 (MQ=255)                      cGACTACGCTGCTGGTGGCGCTGCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAggg                                                                                                                                                                                                                      >  2:351931/1‑250 (MQ=255)                                   gttggctcttcgcACTTTCGCCAACCTGTTGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCGGTACGACAAGATGCGTTTGTTCATCGATGACAGGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGTCAGGCAATG‑AAAAATCTCTCTTCAGGGCGCGTAAATGTGCTAGCGCATGCAGAAGCTTTTCGCGGTTTAGGAGTAGAAATTaa                                                                                                                                                                                                           <  2:418075/239‑1 (MQ=255)                                   ggtggCGCTGCGCACTCTCGCCAACCTTTTGCGTTATGAGAAGCAAAGCCGCAACGCCCAGCAAATCGCCGTTCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAAc                                                                                                                                                                                                         <  2:349876/251‑1 (MQ=255)                                   ggtggCGCTGCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTaaa                                                                                                                                                                                                          <  2:228371/250‑1 (MQ=255)                                        ggcTGCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTAGCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCgag                                                                                                                                                                                                    <  2:287785/250‑1 (MQ=255)                                          cTGCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAAcgcg                                                                                                                                                                                                      >  2:136334/1‑247 (MQ=255)                                           tGCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAACTCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGACATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGAtt                                                                                                                                                                                                 >  1:372503/1‑251 (MQ=255)                                            gcgcACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTa                                                                                                                                                                                                <  1:232595/251‑1 (MQ=255)                                                   aTCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGa                                                                                                                                                                                         <  1:310320/251‑1 (MQ=255)                                                        cAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTgg                                                                                                                                                                                    >  1:140258/1‑251 (MQ=255)                                                          aCCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAACTAGCCGATCGTGACAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGCAATGGGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCt                                                                                                                                                                                  >  1:65386/1‑251 (MQ=255)                                                            cggtgGCGGTAGGAGCATCAAAGCCGCAACGCCCAGCAAATAGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGa                                                                                                                                                                                <  2:435735/249‑1 (MQ=255)                                                                             tCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATAGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAgg                                                                                                                                                               >  1:161917/1‑251 (MQ=255)                                                                                aaGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTGCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATg                                                                                                                                                            <  1:360896/251‑1 (MQ=255)                                                                                         cGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATc                                                                                                                                                  <  1:166788/251‑1 (MQ=255)                                                                                              aGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACCTGTCCGCGATTGGTCAAAGTCTCGACAACGCGCAGGATAATTATCGGCAGCCCATG‑AAAAAACTCTCTTCCGGGCGCGGAAATGGGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTTGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCCGGCTGAA‑GGGTATCGACtt                                                                                                                                             >  2:137866/1‑251 (MQ=255)                                                                                                                ccAGCAAGCTGTACGACAAGATGAGTTTGTTAATCGATACCATGTCCGCGATTGGCCAAAGTCTCGCAAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAAATATCTTCAGGGCGCGGAAACGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTCAACGCGAAATTAATCCGGATTTGGCTGAACAGACGGTGAGACAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGagc                                                                                                                             >  1:271554/1‑249 (MQ=255)                                                                                                                    cAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGTGGCAGGATGAA‑GATTATCGACTTCGGTCGGTTCCGGAGCAGCCGa                                                                                                                       <  1:177791/251‑1 (MQ=255)                                                                                                                     aaGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGaa                                                                                                                      <  1:195203/251‑1 (MQ=255)                                                                                                                                 aaGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATGGAAAAAACCATCTCCAGGGCGCGGAATTGTGTTAGCGAAGCCGGAAGGGTTTCGGGGTTTAGGAGAAGAAATTAACGGCGAGATAATTCCGGATTTGCCTGAACAGGCGGTAGGCCAGG‑TGAAAGGGTAACGACTTCGGCCGGTTCCGGAGCAGCCGatgaataaagctt                                                                                                           >  1:83067/1‑239 (MQ=255)                                                                                                                                      gCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAAGCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCaa                                                                                                      <  1:157796/250‑1 (MQ=255)                                                                                                                                        gTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAAcgc                                                                                                   <  1:225663/251‑1 (MQ=255)                                                                                                                                                     tGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAatat                                                                                       >  1:426084/1‑250 (MQ=255)                                                                                                                                                            tCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAat                                                                                         >  1:75533/1‑241 (MQ=255)                                                                                                                                                            tCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAat                                                                                         <  2:75533/241‑1 (MQ=255)                                                                                                                                                            tCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGACGCATGTTGAAGCTTATCAACGCGTTGATGAATATAATcagc                                                                                >  2:7190/1‑250 (MQ=255)                                                                                                                                                             ccGCGGTTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATcagcag                                                                              <  2:256089/251‑1 (MQ=255)                                                                                                                                                                               ccGACAAAGCGCAGGATAATTCTCGGCAGGCAATG‑AAAAAACTCTCTGCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTGTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCTGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTggg                                                            <  2:267336/250‑1 (MQ=255)                                                                                                                                                                               ccGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCTAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTggg                                                            <  2:77070/250‑1 (MQ=255)                                                                                                                                                                                cGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTgg                                                             <  1:427251/249‑1 (MQ=255)                                                                                                                                                                                          gaGGATAATTATCGGCAAGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGc                                                 <  1:242777/250‑1 (MQ=255)                                                                                                                                                                                                 aTTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGt                                            >  1:321986/1‑249 (MQ=255)                                                                                                                                                                                                        gcaggcaATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCAGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGGAACGACTTCGGGCGGTTCCGGAGCAGCCGAATGATGAAGCTTAACAACGCGAGGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGAAGTTAAGCCGGGTAGAAATCt                                    >  1:302545/1‑250 (MQ=255)                                                                                                                                                                                                             caATG‑ATAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGc                               >  1:116274/1‑250 (MQ=255)                                                                                                                                                                                                              aagg‑AAAAAACTCTCTTCAGTGCGCGGAAAGGGGCTAGCGCATGCAGAAGTGTTTCGCGGTTTAGGATTATAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTTAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGAGGATGAATATAAGCAGCAGTCGCGCTAGTCCATTGTGTGGAGTGAAGCCGGGTAGAAATCTAGGGCAt                             <  2:169797/248‑1 (MQ=255)                                                                                                                                                                                                               aTG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCAGTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATc                            <  1:402913/251‑1 (MQ=255)                                                                                                                                                                                                                             cgtCAGGGCGCGGAAATGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAAGATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAAACTAGGGCATCGACGCCCAATCTg               <  2:87518/249‑1 (MQ=255)                                                                                                                                                                                                                                          aaaTGTGCTAGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGa  <  1:136334/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   ACAACATCATGCTGGTTAGCCCGACTACGCTGCTGGTGGCGCTGCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATCGTGCCAGCAAGCTGTACGACAAGATGCGTTTGTTCATCGATGACATGTCCGCGATTGGTCAAAGTCTCGACAAAGCGCAGGATAATTATCGGCAGGCAATG‑AAAAAACTCTCTTCAGGGCGCGGAAATGTGCTGGCGCAGGCAGAAGCGTTTCGCGGTTTAGGAGTAGAAATTAAACGCGAGATTAATCCGGATTTGGCTGAACAGGCGGTGAGCCAGGATGAA‑GAGTATCGACTTCGGTCGGTTCCGGAGCAGCCGAATGATGAAGCTTATCAACGCGATGATGAATATAATCAGCAGTCGCGCTAGCCCATTGGGAGTAGTTAAGCCGGGTAGAAATCTAGGGCATCGACGCCCAATCTGTTACACTTCTGGA  >  NC_000913/4018345‑4018827

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