Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,466,996	T→C	pseudogene (1605/2513 nt)	ydbA →	pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,466,996	0	T	C	100.0%	328.6 / NA	92	pseudogene (1605/2513 nt)	ydbA	pseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (47/45);  total (47/45)

CAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTC  >  NC_000913/1466868‑1467124                                                                                                                                 |                                                                                                                                 cAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                        <  2:454198/139‑1 (MQ=255)              aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  2:862003/139‑1 (MQ=255)              aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  1:1245213/139‑1 (MQ=255)              aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  2:1245213/1‑139 (MQ=255)               tAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGc                                                                                                          >  1:57721/1‑139 (MQ=255)                  cGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCAtt                                                                                                       >  2:1684072/1‑139 (MQ=255)                            agCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTgg                                                                                             >  1:758369/1‑139 (MQ=255)                              cACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGAt                                                                                           <  1:591035/139‑1 (MQ=255)                                 tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGc                                                                                                                   <  2:1750306/112‑1 (MQ=255)                                 tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGc                                                                                                                   >  1:1750306/1‑112 (MQ=255)                                      ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCt                                                                                   >  2:1478533/1‑139 (MQ=255)                                       cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGc                                                                                                                   <  2:1711814/106‑1 (MQ=255)                                       cAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGc                                                                                                                   >  1:1711814/1‑106 (MQ=255)                                        aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTaa                                                                                                   <  2:1318347/121‑1 (MQ=255)                                        aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTaa                                                                                                   >  1:1318347/1‑121 (MQ=255)                                         aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGAc                                                                                >  2:173680/1‑139 (MQ=255)                                           ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGcc                                                                                    <  2:427952/133‑1 (MQ=255)                                           ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGcc                                                                                    >  1:427952/1‑133 (MQ=255)                                           ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTg                                                                              >  2:228167/1‑139 (MQ=255)                                           ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTg                                                                              >  2:497490/1‑139 (MQ=255)                                             cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGATTAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAAACTGGATGGTAGCCTGACTGtt                                                                            >  1:1157683/1‑139 (MQ=255)                                             cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCgg                                                                                                                     >  1:524973/1‑98 (MQ=255)                                             cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCgg                                                                                                                     <  2:524973/98‑1 (MQ=255)                                                 accggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTa                                                                                                    >  2:361528/1‑111 (MQ=255)                                                 accggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTa                                                                                                    <  1:361528/111‑1 (MQ=255)                                                     gcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTg                                                                                  <  2:721257/125‑1 (MQ=255)                                                     gcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTg                                                                                  >  1:721257/1‑125 (MQ=255)                                                     gcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCgg                                                                    <  1:1478533/139‑1 (MQ=255)                                                           aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  1:228167/139‑1 (MQ=255)                                                           aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  2:758369/139‑1 (MQ=255)                                                           aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  1:497490/139‑1 (MQ=255)                                                             aTCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATg                                                            <  1:1684072/139‑1 (MQ=255)                                                               cAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGa                                                                                 >  1:1374750/1‑116 (MQ=255)                                                               cAACGGTAACGACGCCACCGCGAATAAAAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGa                                                                                 <  2:1374750/116‑1 (MQ=255)                                                                  cGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACt                                                                               <  1:918737/115‑1 (MQ=255)                                                                  cGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACt                                                                               >  2:918737/1‑115 (MQ=255)                                                                     tAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              >  2:373786/1‑129 (MQ=255)                                                                     tAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  1:373786/129‑1 (MQ=255)                                                                       gcgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATAGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                      >  2:1153876/2‑71 (MQ=255)                                                                       acgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                      <  1:1153876/71‑1 (MQ=255)                                                                       acgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGaa                                                  >  1:1541482/1‑139 (MQ=255)                                                                           cgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              >  2:218911/1‑123 (MQ=255)                                                                           cgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  1:218911/123‑1 (MQ=255)                                                                             ccACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                        >  1:574872/1‑63 (MQ=255)                                                                             ccACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                        <  2:574872/63‑1 (MQ=255)                                                                             ccACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTgg                                            >  2:18191/1‑139 (MQ=255)                                                                             ccACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTgg                                            >  1:227753/1‑139 (MQ=255)                                                                                 cgcgAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGAc                                        <  2:57721/139‑1 (MQ=255)                                                                                     aaTAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  2:896155/1‑68 (MQ=255)                                                                                     aaTAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  1:896155/68‑1 (MQ=255)                                                                                             tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                      <  2:1370858/49‑1 (MQ=255)                                                                                             tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                      >  1:1370858/1‑49 (MQ=255)                                                                                             tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTgc                                                                 <  2:797571/102‑1 (MQ=255)                                                                                             tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTgc                                                                 >  1:797571/1‑102 (MQ=255)                                                                                             tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGtt                            <  2:227753/139‑1 (MQ=255)                                                                                             tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGtt                            <  1:173680/139‑1 (MQ=255)                                                                                              ggAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCggg                                                                                                                    >  2:26310/1‑50 (MQ=255)                                                                                              ggAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCggg                                                                                                                    <  1:26310/50‑1 (MQ=255)                                                                                              ggAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  2:1318125/1‑59 (MQ=255)                                                                                              ggAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  1:1318125/59‑1 (MQ=255)                                                                                               gAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  1:1182138/58‑1 (MQ=255)                                                                                               gAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  2:1182138/1‑58 (MQ=255)                                                                                               gAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaa                          >  1:757300/1‑139 (MQ=255)                                                                                                  aaaCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                        <  1:913177/42‑1 (MQ=255)                                                                                                  aaaCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                        >  2:913177/1‑42 (MQ=255)                                                                                                    aCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAacaaca                     >  2:1731038/1‑139 (MQ=255)                                                                                                     ccACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAAcaacaa                    <  2:757300/139‑1 (MQ=255)                                                                                                       aCTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGACAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgt                                                          >  2:571277/1‑99 (MQ=255)                                                                                                       aCTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGACAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgt                                                          <  1:571277/99‑1 (MQ=255)                                                                                                        cTGTTGATGGAAAAGATGCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTCGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAgg                 <  2:1157683/139‑1 (MQ=255)                                                                                                         tGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCAt                                                                                                        <  2:34731/51‑1 (MQ=255)                                                                                                         tGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCAt                                                                                                        >  1:34731/1‑51 (MQ=255)                                                                                                         tGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAAt                                                                                                 >  1:131585/1‑58 (MQ=255)                                                                                                         tGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAAt                                                                                                 <  2:131585/58‑1 (MQ=255)                                                                                                           ttGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGtt                                                       <  2:1328859/98‑1 (MQ=255)                                                                                                           ttGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGgaga              >  1:290455/1‑139 (MQ=255)                                                                                                           ttGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCAGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGtt                                                       >  1:1328859/1‑98 (MQ=255)                                                                                                              aTGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGatat           >  1:1871441/1‑139 (MQ=255)                                                                                                                   aaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  1:565705/83‑1 (MQ=255)                                                                                                                   aaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              >  2:565705/1‑83 (MQ=255)                                                                                                                     aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGc                                                                      <  2:1794461/73‑1 (MQ=255)                                                                                                                     aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGc                                                                      >  1:1794461/1‑73 (MQ=255)                                                                                                                     aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  2:783942/81‑1 (MQ=255)                                                                                                                     aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              >  1:783942/1‑81 (MQ=255)                                                                                                                     aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGtt                            >  2:1073483/1‑115 (MQ=255)                                                                                                                     aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGtt                            <  1:1073483/115‑1 (MQ=255)                                                                                                                     aGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTgttgtt    >  1:750949/1‑139 (MQ=255)                                                                                                                       aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc  <  1:1731038/139‑1 (MQ=255)                                                                                                                        ttCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAgg                 >  1:33079/1‑123 (MQ=255)                                                                                                                        ttCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAgg                 <  2:33079/123‑1 (MQ=255)                                                                                                                            aCGGGCACGAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATg                                                                                          <  1:485731/46‑1 (MQ=255)                                                                                                                            aCGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATg                                                                                          >  2:485731/1‑46 (MQ=255)                                                                                                                                 |                                                                                                                                 CAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTC  >  NC_000913/1466868‑1467124

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