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%	257.9 / NA	84	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/44);  total (40/44)

GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACT  >  NC_000913/4499375‑4499634                                                                                                                              |                                                                                                                                       gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           <  2:695332/139‑1 (MQ=11)      ggCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                      <  2:1120030/139‑1 (MQ=21)          ttCGAGCGTTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATACCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGACGAATCCACCTTGTGCTGa                                                                                                                  >  2:847155/1‑139 (MQ=11)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       <  1:345589/139‑1 (MQ=255)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       <  2:877624/139‑1 (MQ=255)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       <  1:1285144/139‑1 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTta                                                                                               >  1:948261/1‑139 (MQ=255)                               cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGata                                                                                                                >  1:60774/1‑120 (MQ=255)                               cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGata                                                                                                                <  2:60774/120‑1 (MQ=255)                               cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCT‑GAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       >  1:1220049/1‑128 (MQ=255)                               cATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCT‑GAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       <  2:1220049/128‑1 (MQ=255)                                 tCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGTCGAATCCACCTTGTGCTGATATGAATACGGATATTTTATACa                                                                                           <  1:847155/139‑1 (MQ=255)                                    gCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGaa                                                                                        <  1:237181/139‑1 (MQ=255)                                       tAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATAc                                                                                     >  1:935820/1‑139 (MQ=255)                                           gcgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           <  2:107131/97‑1 (MQ=14)                                           gcgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           >  1:107131/1‑97 (MQ=14)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           <  1:972022/96‑1 (MQ=14)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           >  2:972022/1‑96 (MQ=14)                                                    tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAg                                                                        >  1:1014369/1‑139 (MQ=255)                                                      tcgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                      <  1:1081054/91‑1 (MQ=255)                                                      tcgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                      >  2:1081054/1‑91 (MQ=255)                                                       cgctcgtCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     <  1:401512/139‑1 (MQ=255)                                                       cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     <  2:691205/139‑1 (MQ=255)                                                       cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     <  2:1014369/139‑1 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCt                                                                                                                          >  2:256375/1‑82 (MQ=17)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCt                                                                                                                          <  1:256375/82‑1 (MQ=17)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGa                                                                                                      <  2:75833/102‑1 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGa                                                                                                      >  1:75833/1‑102 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGt                                                                 >  1:84648/1‑139 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTAGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       >  2:1293816/1‑101 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTAGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       <  1:1293816/101‑1 (MQ=255)                                                                 gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGCCAAATCCACCTTGTACTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACa                                                           >  2:559667/1‑139 (MQ=255)                                                                        tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                                    >  1:1177458/1‑139 (MQ=255)                                                                        tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                                    <  2:948261/139‑1 (MQ=255)                                                                          tGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCat                                                  <  1:630712/139‑1 (MQ=255)                                                                          tGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCat                                                  <  1:559667/139‑1 (MQ=255)                                                                            cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCatat                                                >  2:503202/1‑139 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCt                                                                                                                    >  2:370728/1‑70 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCt                                                                                                                    <  1:370728/70‑1 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                             <  1:271955/125‑1 (MQ=255)                                                                             ggCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                             >  2:271955/1‑125 (MQ=255)                                                                                         gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCtgtg                                   >  1:1143389/1‑139 (MQ=255)                                                                                            aTGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATc                                <  1:503202/139‑1 (MQ=255)                                                                                              gggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAg                                                          <  1:765548/111‑1 (MQ=255)                                                                                              gggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAACACTTTTCAGTACAAGCGACTGTATTACAg                                                          >  2:765548/1‑111 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           >  1:25087/1‑45 (MQ=12)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                           <  2:25087/45‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATg                                                                                                              >  1:1131690/1‑58 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATg                                                                                                              <  2:1131690/58‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       <  1:794853/65‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       >  2:794853/1‑65 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     <  2:879921/99‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     >  1:1048816/1‑99 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     <  2:1048816/99‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                                         >  2:1161903/1‑111 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                                         <  1:1161903/111‑1 (MQ=255)                                                                                               ggTCAAGTGATAACAGATGTATGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                                     >  1:879921/1‑99 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                      <  1:541051/45‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                      >  2:541051/1‑45 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATg                                                                                                              <  2:930362/53‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATg                                                                                                              >  1:930362/1‑53 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGt                   >  2:390514/1‑139 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCAGCCATATATGATTACCTGAGATCAGGATGAATAAGt                   >  2:963327/1‑139 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                      >  1:347724/1‑39 (MQ=17)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                                      <  2:347724/39‑1 (MQ=17)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACa                                                                                           <  2:585538/66‑1 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACa                                                                                           >  1:585538/1‑66 (MQ=255)                                                                                                              gATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTg              >  2:1236499/1‑139 (MQ=255)                                                                                                              gATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTg              <  2:935820/139‑1 (MQ=255)                                                                                                              gATGTCTGGAAATATAAGGGCAAATCCACCTTGGGCTGATATGAATACGGATATTTTATACAGAGTACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGGGATCAGGAAGAATAAGTGGCTg              <  2:289630/139‑1 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                                 >  2:749084/1‑39 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                                 <  1:749084/39‑1 (MQ=21)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGa             >  1:696108/1‑139 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACa                                                                                           <  2:874209/60‑1 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                       <  2:1163460/80‑1 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                                       >  1:1163460/1‑80 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAg            <  2:1143389/139‑1 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATATTATACa                                                                                           >  1:874209/1‑60 (MQ=255)                                                                                                                  tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                                    <  1:321708/97‑1 (MQ=255)                                                                                                                  tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGcc                                                    >  2:321708/1‑97 (MQ=255)                                                                                                                  tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGt          >  1:609383/1‑139 (MQ=255)                                                                                                                    tGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       >  1:600746/1‑44 (MQ=255)                                                                                                                    tGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                                       <  2:600746/44‑1 (MQ=255)                                                                                                                          tataAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACt  <  1:963327/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                       GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAGGTTTCAAACT  >  NC_000913/4499375‑4499634

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