Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,499,500	G→A	*302* (TAG→TAA)	insD1 →	IS2 transposase TnpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,499,500	0	G	A	100.0%	263.4 / NA	86	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/46);  total (40/46)

GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCA  >  NC_000913/4499375‑4499630                                                                                                                              |                                                                                                                                   gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                       <  1:1530973/139‑1 (MQ=11)      ggCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  2:846199/139‑1 (MQ=21)      ggCGTTCGAGCACTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTAGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  1:1091793/139‑1 (MQ=255)        cGTTCGAGCATTATGACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAAGATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCt                                                                                                                <  1:847633/139‑1 (MQ=255)          ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                       <  1:730032/130‑1 (MQ=11)          ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                       >  2:730032/1‑130 (MQ=11)          ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  >  2:1544597/1‑135 (MQ=21)          ttCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  1:1544597/135‑1 (MQ=21)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                   <  2:1044262/139‑1 (MQ=255)                     aagggAAAGGCATCGGAATAGTGGGTTGGGTTATGGCTCGCAAGGGGAATTTCTGCGGCAGGGGGTTTGAAATGGTTTAAGTAATAAAAGATGTTTGGAAATATAAGGGCAAACCCACCTTGTGCTGATATGAATACgg                                                                                                   <  1:3397/135‑1 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATGTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  >  1:1120405/1‑116 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATGTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGAAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  2:1120405/116‑1 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  >  2:1274078/1‑116 (MQ=25)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  >  2:354116/1‑116 (MQ=25)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  1:1274078/116‑1 (MQ=25)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  1:354116/116‑1 (MQ=25)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                   <  2:866264/131‑1 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                   >  1:866264/1‑131 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTta                                                                                           >  1:346002/1‑139 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTta                                                                                           >  2:661285/1‑139 (MQ=255)                                         gTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtaaaacctgctttt                                                                               >  2:3397/1‑127 (MQ=17)                                          tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACttt                                                                              >  1:172531/1‑139 (MQ=255)                                              cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                       >  1:886103/1‑94 (MQ=14)                                              cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                       <  2:886103/94‑1 (MQ=14)                                              cTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAg                                                                          >  1:211261/1‑139 (MQ=255)                                                  gTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACa                                                                      <  2:211261/139‑1 (MQ=255)                                                  gTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACa                                                                      >  1:649471/1‑139 (MQ=255)                                                    tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  >  2:536288/1‑93 (MQ=255)                                                    tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  1:536288/93‑1 (MQ=37)                                                      tcgctcgcCACGGGAATATCTGCGGCAGCGGACTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGata                                                                                                            <  1:862709/97‑1 (MQ=255)                                                      tcgctcgcCACGGGAATATCTGCGGCAGCGGACTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGata                                                                                                            >  2:862709/1‑97 (MQ=255)                                                       cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                 <  1:661285/139‑1 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCt                                                                                                                >  2:1527347/1‑88 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCt                                                                                                                <  1:1527347/88‑1 (MQ=255)                                                             ccACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  >  2:735426/1‑84 (MQ=255)                                                             ccACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  1:735426/84‑1 (MQ=255)                                                               aCGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                         <  1:479953/139‑1 (MQ=255)                                                                   gAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAACTCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                                     <  1:599015/139‑1 (MQ=255)                                                                   gAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                                     <  2:172531/139‑1 (MQ=255)                                                                      tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                  <  1:503735/139‑1 (MQ=255)                                                                      tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                  <  1:383264/139‑1 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                  >  2:1678327/1‑132 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                  <  1:1678327/132‑1 (MQ=255)                                                                                 gCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                       <  2:649471/139‑1 (MQ=255)                                                                                 gCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                       >  1:1578518/1‑139 (MQ=255)                                                                                  cGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAtt                                      >  2:1176165/1‑139 (MQ=255)                                                                                    ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  >  2:673874/1‑61 (MQ=255)                                                                                    ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  1:673874/61‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATAc                                                                                                     >  1:1492350/1‑63 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATAc                                                                                                     <  2:1492350/63‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTa                                                                        <  2:326522/92‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTa                                                                        >  1:326522/1‑92 (MQ=255)                                                                                                gTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                   <  2:95368/64‑1 (MQ=255)                                                                                                gTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                   >  1:95368/1‑64 (MQ=255)                                                                                                   aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGTTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGa                              >  2:173004/1‑130 (MQ=255)                                                                                                   aaGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGTTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGa                              <  1:173004/130‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                                    <  2:365385/91‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                                    >  1:365385/1‑91 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                 <  2:567906/94‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                 >  1:567906/1‑94 (MQ=255)                                                                                                      tGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATa                  <  1:1176165/139‑1 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  <  2:601302/40‑1 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                  >  1:601302/1‑40 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                   >  1:207499/1‑55 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                   <  2:207499/55‑1 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTATTCAGTACAAGCGACTGTATTACAGTCTGCTATATATGATTACCTGTGATCAGGAAGAATAAGt               >  1:1408166/1‑139 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                             <  1:932304/44‑1 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                             >  2:932304/1‑44 (MQ=255)                                                                                                           aCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                  <  2:1276883/102‑1 (MQ=255)                                                                                                           aCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                                  >  1:1276883/1‑102 (MQ=255)                                                                                                              gATGTCTGGAAATATAAGGGCAAAGCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTg          <  1:1797969/139‑1 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGAt                             <  2:484760/119‑1 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGAt                             >  1:484760/1‑119 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGa         >  1:19720/1‑139 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                 <  1:342840/82‑1 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                 >  2:342840/1‑82 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAg        <  2:748310/139‑1 (MQ=255)                                                                                                                 gTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                         >  2:1246028/1‑57 (MQ=255)                                                                                                                 gTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                         <  1:1246028/57‑1 (MQ=255)                                                                                                                  tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACtttt                                                                             >  2:71240/1‑68 (MQ=255)                                                                                                                  tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACtttt                                                                             <  1:71240/68‑1 (MQ=255)                                                                                                                   cTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                   >  1:1449205/1‑45 (MQ=255)                                                                                                                   cTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                   <  2:1449205/45‑1 (MQ=255)                                                                                                                      gAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACa                                                                                       >  2:1262383/1‑54 (MQ=255)                                                                                                                      gAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACa                                                                                       <  1:1262383/54‑1 (MQ=255)                                                                                                                      gAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCa  >  2:1265056/1‑139 (MQ=255)                                                                                                                              |                                                                                                                                   GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCA  >  NC_000913/4499375‑4499630

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