Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,333,691	A→G	pseudogene (712/4605 nt)	yfaS ←	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family, putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,333,691	0	A	G	100.0%	288.6 / NA	82	pseudogene (712/4605 nt)	yfaS	pseudogene, bacterial alpha2‑macroglobulin YfaS variant family, putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/42);  total (40/42)

CTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGA  >  NC_000913/2333569‑2333825                                                                                                                           |                                                                                                                                       cTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCttgt                                                                                                                        <  2:1174706/139‑1 (MQ=255) cTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCttgt                                                                                                                        <  2:1653229/139‑1 (MQ=255)     tttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGc                                                                                                                    >  2:825964/1‑139 (MQ=255)         tACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                <  1:825964/139‑1 (MQ=255)         tACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                <  2:1663610/139‑1 (MQ=255)         tACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                <  1:1948706/139‑1 (MQ=255)         tACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                <  2:52720/139‑1 (MQ=255)          aCCCGCGGTCCACACCAGAAGCTCTTTGCCTTACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGc                                                                                                                              <  2:1764918/124‑1 (MQ=255)          aCCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGAGCAGTCCCGGc                                                                                                                              >  1:1764918/1‑124 (MQ=255)                ggTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAata                                                                                                         >  1:1211105/1‑139 (MQ=255)                 gTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAatat                                                                                                        <  2:1211105/139‑1 (MQ=255)                            gaGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTc                                                                                                                                   >  1:1286504/2‑101 (MQ=255)                            aaGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGTCTCGACGAGGTGCAGTc                                                                                                                                   <  2:1286504/101‑1 (MQ=255)                            aaGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGagag                                                                                                              <  2:1324529/122‑1 (MQ=255)                            aaGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGagag                                                                                                              >  1:1324529/1‑122 (MQ=255)                              gCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCgg                                                                                                                               <  1:985213/103‑1 (MQ=255)                              gCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCgg                                                                                                                               >  2:985213/1‑103 (MQ=255)                                 ctTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTcc                                                                                                                                  <  2:495227/97‑1 (MQ=255)                                 ctTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTcc                                                                                                                                  >  1:495227/1‑97 (MQ=255)                                 ctTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTcc                                                                                                                                  <  2:1521460/97‑1 (MQ=255)                                 ctTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTcc                                                                                                                                  >  1:1521460/1‑97 (MQ=255)                                 ctTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTg                                                                                        >  2:1989663/1‑139 (MQ=255)                                   ttGCCTTACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcg                                                                                      >  1:860040/1‑139 (MQ=255)                                   ttGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGcc                                                                                                                   <  2:829849/110‑1 (MQ=255)                                   ttGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGcc                                                                                                                   >  1:829849/1‑110 (MQ=255)                                   ttGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCt                                                                                         >  2:799641/1‑136 (MQ=255)                                   ttGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCt                                                                                         <  1:799641/136‑1 (MQ=255)                                   ttGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcg                                                                                      >  2:670667/1‑139 (MQ=255)                                    tGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcgg                                                                                     >  2:1132404/1‑139 (MQ=255)                                              tttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGcc                                                                                                                   >  2:2009046/1‑99 (MQ=255)                                              tttGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGcc                                                                                                                   <  1:2009046/99‑1 (MQ=255)                                                tGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGagag                                                                                                              >  1:733520/1‑102 (MQ=255)                                                tGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGagag                                                                                                              <  2:733520/102‑1 (MQ=255)                                                      gcgcCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGa                                                                   >  2:777431/1‑139 (MQ=255)                                                        gcCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATg                                                                 <  1:586672/139‑1 (MQ=255)                                                          cACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgcg                                                               >  1:1364387/1‑139 (MQ=255)                                                          cACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgcg                                                               <  1:1932061/139‑1 (MQ=255)                                                          cACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgcg                                                               <  2:860040/139‑1 (MQ=255)                                                            cGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGcc                                                             >  2:927561/1‑139 (MQ=255)                                                                 tCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcggcg                                                                                   >  1:670615/1‑112 (MQ=255)                                                                 tCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcggcg                                                                                   <  2:670615/112‑1 (MQ=255)                                                                 tCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGACCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGAGCCTTCCa                                                        >  2:1223280/1‑139 (MQ=255)                                                                   ggAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGt                                                      >  1:671599/1‑139 (MQ=255)                                                                    gAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCga                                                                                  >  1:994866/1‑110 (MQ=255)                                                                    gAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCga                                                                                  <  2:994866/110‑1 (MQ=255)                                                                      aacaaacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGtttt                                                   <  1:777431/139‑1 (MQ=255)                                                                        caaacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTc                                                                                                                                   <  1:748104/57‑1 (MQ=255)                                                                        caaacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTc                                                                                                                                   >  2:748104/1‑57 (MQ=255)                                                                         aaacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTc                                                          >  2:653825/1‑129 (MQ=255)                                                                         aaacaCCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTc                                                          <  1:653825/129‑1 (MQ=255)                                                                               caccGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGagag                                                                                                              <  1:1758745/71‑1 (MQ=255)                                                                               caccGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGagag                                                                                                              >  2:1758745/1‑71 (MQ=255)                                                                               caccGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAatat                                                                                                        >  2:1620396/1‑77 (MQ=255)                                                                               caccGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAatat                                                                                                        <  1:1620396/77‑1 (MQ=255)                                                                                  cgtcgtCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgcg                                                               <  2:558029/115‑1 (MQ=255)                                                                                  cgtcgtCGCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATgcg                                                               >  1:558029/1‑115 (MQ=255)                                                                                    tcgtcgCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCgg                                     >  2:995665/1‑139 (MQ=255)                                                                                      gtcgCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                <  1:720956/62‑1 (MQ=255)                                                                                      gtcgCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                >  2:720956/1‑62 (MQ=255)                                                                                         gCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcgg                                                                                     >  1:252273/1‑86 (MQ=255)                                                                                         gCCCGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcgg                                                                                     <  2:252273/86‑1 (MQ=255)                                                                                            cGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGagag                                                                                                              >  2:1117023/1‑58 (MQ=255)                                                                                            cGATACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGagag                                                                                                              <  1:1117023/58‑1 (MQ=255)                                                                                              aTACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                <  1:1555452/54‑1 (MQ=255)                                                                                              aTACCCACCAACCATCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                >  2:1555452/1‑54 (MQ=255)                                                                                                     ccaaccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCt                    >  2:211397/1‑139 (MQ=255)                                                                                                      caaccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTg                   <  2:671599/139‑1 (MQ=255)                                                                                                        accaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCTAGTTTTACCCCTTGCTGCGGCTCGAATGGCTTAGCCTGcc                 >  2:1074738/1‑139 (MQ=255)                                                                                                         ccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                >  1:665169/1‑43 (MQ=255)                                                                                                         ccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCgag                                                                                                                <  2:665169/43‑1 (MQ=255)                                                                                                         ccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcgg                                                                                     <  2:1341854/70‑1 (MQ=255)                                                                                                         ccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTgcgg                                                                                     >  1:1341854/1‑70 (MQ=255)                                                                                                         ccaTCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCa                <  1:211397/139‑1 (MQ=255)                                                                                                            tCGCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAg             <  1:927561/139‑1 (MQ=255)                                                                                                              gCCTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTg           <  1:670667/139‑1 (MQ=255)                                                                                                                cTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGaa                                                                                                           >  2:1516202/1‑41 (MQ=255)                                                                                                                cTCGACGAGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGaa                                                                                                           <  1:1516202/41‑1 (MQ=255)                                                                                                                  cgacgaGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCt                                                            >  2:2025764/1‑86 (MQ=255)                                                                                                                  cgacgaGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCt                                                            <  1:2025764/86‑1 (MQ=255)                                                                                                                   gacgaGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCAGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCa                <  2:1216479/129‑1 (MQ=255)                                                                                                                   gacgaGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCAGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCa                >  1:1216479/1‑129 (MQ=255)                                                                                                                       aGGTGCAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGa  <  1:995665/139‑1 (MQ=255)                                                                                                                           |                                                                                                                                       CTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGAGGTACAGTCCCGGCTCTTGTTGGCCGAGAGGAATATAAATGTTACCCGGCTGCGGCGAGATGAAATTGCTGGATGCGCCTTCCAGTTTTACCCCTTGCTGCGGCTCGAACGGTTTAGCCTGCCATAGTGGATAACGGA  >  NC_000913/2333569‑2333825

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