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

GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAG  >  NC_000913/4499375‑4499624                                                                                                                              |                                                                                                                             gCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  2:503609/139‑1 (MQ=11)       gCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGc                                                                                                           >  2:493935/1‑139 (MQ=25)       gCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGc                                                                                                           <  2:489778/139‑1 (MQ=25)           tCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       <  2:45218/139‑1 (MQ=255)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:862318/139‑1 (MQ=255)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:1306018/139‑1 (MQ=255)                     taACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:37815/139‑1 (MQ=255)                           aTGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATAttt                                                                                       >  2:959269/1‑139 (MQ=255)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            >  1:1011891/1‑116 (MQ=25)                             ggCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            <  2:1011891/116‑1 (MQ=25)                                          tGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACttt                                                                        >  2:190304/1‑139 (MQ=255)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTc                                                                      >  1:59905/1‑139 (MQ=255)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTc                                                                      >  2:603435/1‑139 (MQ=255)                                            cgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTc                                                                      >  1:142897/1‑139 (MQ=255)                                               tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATa                                                                                                >  1:245147/1‑110 (MQ=255)                                               tGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATa                                                                                                <  2:245147/110‑1 (MQ=255)                                                    tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       <  1:1005723/98‑1 (MQ=255)                                                    tATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       >  2:1005723/1‑98 (MQ=255)                                                     aTCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             <  2:142897/139‑1 (MQ=255)                                                     aTCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             <  1:603435/139‑1 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 >  2:683651/1‑81 (MQ=14)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  1:683651/81‑1 (MQ=14)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            <  2:900116/86‑1 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             >  1:375277/1‑133 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGc                                                             <  2:375277/133‑1 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGt                                                       >  2:379180/1‑139 (MQ=255)                                                           cgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGAGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            >  1:900116/1‑86 (MQ=255)                                                            gcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            <  2:673270/85‑1 (MQ=255)                                                            gcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            >  1:673270/1‑85 (MQ=255)                                                                 gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACa                                                 >  1:1324522/1‑139 (MQ=255)                                                                  ggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAg                                                <  2:1324522/139‑1 (MQ=255)                                                                   gAATATCTGCGGCAGCGGACTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTTTATTACAGt                                               <  1:903943/139‑1 (MQ=255)                                                                    aaTATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                            >  1:1254421/1‑109 (MQ=255)                                                                    aaTATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                            <  2:1254421/109‑1 (MQ=255)                                                                      tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                            <  2:27218/139‑1 (MQ=255)                                                                      tatCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTg                                            <  2:59905/139‑1 (MQ=255)                                                                        tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTc                                              <  2:204699/135‑1 (MQ=255)                                                                        tCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTc                                              >  1:204699/1‑135 (MQ=255)                                                                           gCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCata                                       >  1:1083706/1‑139 (MQ=255)                                                                            cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGata                                                                                                      >  1:57076/1‑75 (MQ=255)                                                                            cGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGata                                                                                                      <  2:57076/75‑1 (MQ=255)                                                                            cGGCAGCGGACTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAg                                                <  2:841237/129‑1 (MQ=255)                                                                            cGGCAGCGGACTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAg                                                >  1:841237/1‑129 (MQ=255)                                                                               cAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       <  1:642347/71‑1 (MQ=255)                                                                               cAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGat                                                                                                       >  2:642347/1‑71 (MQ=255)                                                                                 gCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGAt                                 <  2:1180318/139‑1 (MQ=255)                                                                                    ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:595171/76‑1 (MQ=255)                                                                                    ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:595171/1‑76 (MQ=255)                                                                                    ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           <  1:166611/110‑1 (MQ=255)                                                                                    ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGa                                                           >  2:166611/1‑110 (MQ=255)                                                                                         gTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCtgtg                         <  1:493935/139‑1 (MQ=255)                                                                                           aaTGGGTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGAt                       <  1:190304/139‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            <  2:890337/50‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            >  1:890337/1‑50 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACt                                                         <  2:1234781/101‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACt                                                         >  1:1234781/1‑101 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                   <  2:1255316/107‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTa                                                   >  1:1255316/1‑107 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                               >  2:455874/1‑111 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                               <  1:455874/111‑1 (MQ=255)                                                                                               ggTTAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAgg                   >  1:658210/1‑139 (MQ=255)                                                                                                 ttAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                   <  1:181215/73‑1 (MQ=255)                                                                                                 ttAAGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                   >  2:181215/1‑73 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 >  1:819748/1‑40 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  2:819748/40‑1 (MQ=11)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 <  2:349290/40‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCAcc                                                                                                                 >  1:349290/1‑40 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:627628/60‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:627628/1‑60 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                               <  1:442168/106‑1 (MQ=255)                                                                                                    aGTGATAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGt                                               >  2:442168/1‑106 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            >  1:1302350/1‑40 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTtgtg                                                                                                            <  2:1302350/40‑1 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                   >  1:660145/1‑65 (MQ=255)                                                                                                         tAACAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTtata                                                                                   <  2:660145/65‑1 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTg        >  1:380075/1‑139 (MQ=255)                                                                                                          aaCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTg        >  1:774413/1‑139 (MQ=255)                                                                                                           aCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:745615/53‑1 (MQ=255)                                                                                                           aCAGATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:745615/1‑53 (MQ=255)                                                                                                              gATGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTg    <  1:959269/139‑1 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  2:426618/1‑49 (MQ=255)                                                                                                               aTGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:426618/49‑1 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  1:353702/48‑1 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  2:353702/1‑48 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                            >  2:649786/1‑65 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATa                                                                            <  1:649786/65‑1 (MQ=255)                                                                                                                tGTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAg  >  2:122812/1‑139 (MQ=255)                                                                                                                 gTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCa                     >  1:496353/1‑119 (MQ=255)                                                                                                                 gTCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCa                     <  2:496353/119‑1 (MQ=255)                                                                                                                  tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             <  2:78345/46‑1 (MQ=255)                                                                                                                  tCTGGAAATATAAGGGCAAATCCACCTTGTGCTGATATGAATACgg                                                                                             >  1:78345/1‑46 (MQ=255)                                                                                                                              |                                                                                                                             GCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCACCTTGTGCTGATATGAATACGGATATTTTATACAGAATACTTTTCAGTACAAGCGACTGTATTACAGTCTGCCATATATGATTACCTGTGATCAGGAAGAATAAGTGGCTGAG  >  NC_000913/4499375‑4499624

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: ‑