Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,458,745	(G)5→4	coding (787/1353 nt)	pmbA →	metalloprotease PmbA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,458,741	0	G	.	100.0%	282.6 / NA	66	coding (783/1353 nt)	pmbA	metalloprotease PmbA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (27/39);  total (27/39)

ATCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGTGGGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACGCCATTC  >  NC_000913/4458602‑4458876                                                                                                                                            |                                                                                                                                        aTCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCtggtg                                                                                                                                         >  2:1017919/1‑140 (MQ=255) aTCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTATGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCtggtg                                                                                                                                         >  2:1436504/1‑140 (MQ=255)    tGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGTgggg                                                                                                                                      <  2:241320/140‑1 (MQ=255)    tGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGTgggg                                                                                                                                      <  1:1430090/140‑1 (MQ=255)       aaaCGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGTggggcga                                                                                                                                   >  1:521777/1‑137 (MQ=255)        aaCGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGTggggcgat                                                                                                                                  <  1:145299/140‑5 (MQ=255)            ccAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCg                                                                                                                             >  2:322282/1‑140 (MQ=255)                agTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTg                                                                                                                         >  2:464798/1‑140 (MQ=255)                agTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTg                                                                                                                         >  1:701437/1‑140 (MQ=255)                  tGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGa                                                                                                                       <  2:1213824/140‑1 (MQ=255)                  tGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGa                                                                                                                       <  1:54459/140‑1 (MQ=255)                   gggTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGAt                                                                                                                      <  1:160316/140‑1 (MQ=255)                   gggTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGAt                                                                                                                      <  1:464798/140‑1 (MQ=255)                   gggTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGAt                                                                                                                      <  1:789757/140‑1 (MQ=255)                   gggTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGAt                                                                                                                      <  2:701437/140‑1 (MQ=255)                   gggTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGAt                                                                                                                      <  2:1242397/140‑1 (MQ=255)                       tGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGt                                                                                                                  <  2:124575/140‑1 (MQ=255)                       tGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGt                                                                                                                  <  1:1337816/140‑1 (MQ=255)                         gggCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGttt                                                                                                                >  1:288395/1‑140 (MQ=255)                          ggCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGtt                                                                                                                 <  1:1124955/138‑1 (MQ=255)                          ggCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGtt                                                                                                                 >  2:1124955/1‑138 (MQ=255)                           gCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTAt                                                                                                              <  1:1017919/140‑1 (MQ=255)                            ccGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATc                                                                                                             <  2:1448650/140‑1 (MQ=255)                               aCTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTa                                                                                                          <  2:1441590/140‑1 (MQ=255)                                    gCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCt                                                                                                     >  2:553417/1‑140 (MQ=255)                                     cTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTa                                                                                                    <  1:701924/140‑1 (MQ=255)                                       cGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTAcc                                                                                                  <  1:35283/140‑1 (MQ=255)                                        gCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCt                                                                                                 <  1:672554/140‑1 (MQ=255)                                          cGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCA‑TTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTcc                                                                                              >  2:1089477/1‑140 (MQ=255)                                          cGTACTTCATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTc                                                                                               >  2:1407363/1‑140 (MQ=255)                                                ttATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCctgctg                                                                                         >  1:787365/1‑140 (MQ=255)                                                    cgcgTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGAtt                                                                                     >  2:1066724/1‑140 (MQ=255)                                                      cgccTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCg                                                                                   <  2:1289004/137‑1 (MQ=255)                                                      cgTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCg                                                                                   <  2:912266/140‑1 (MQ=255)                                                      cgTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCg                                                                                   >  1:254345/1‑140 (MQ=255)                                                      cgTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCg                                                                                   >  1:285429/1‑140 (MQ=255)                                                      cgTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCg                                                                                   >  1:861135/1‑140 (MQ=255)                                                        tCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCt                                                                                 <  2:285429/140‑1 (MQ=255)                                                        tCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCt                                                                                 <  1:628203/140‑1 (MQ=255)                                                        tCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCt                                                                                 >  1:1180582/1‑140 (MQ=255)                                                              aCCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTaa                                                                           <  1:1066724/140‑1 (MQ=255)                                                                   gTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACaaa                                                                      <  1:460262/140‑1 (MQ=255)                                                                           tcCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCg                                                              >  1:927237/1‑140 (MQ=255)                                                                                 aTGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTgg                                                        <  2:1419192/140‑1 (MQ=255)                                                                                  tGAAAGCGCCAGTCA‑TTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCt                                                      <  1:1089477/140‑1 (MQ=255)                                                                                   gAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCt                                                      <  2:1180582/140‑1 (MQ=255)                                                                                   gAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCt                                                      <  1:1407363/140‑1 (MQ=255)                                                                                     aaGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGa                                                    <  2:521777/140‑1 (MQ=255)                                                                                     aaGCGCCAGTCATTTTTGCCAACGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGa                                                    >  1:469943/1‑140 (MQ=255)                                                                                      aGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGAc                                                   <  2:658971/140‑1 (MQ=255)                                                                                      aGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGAc                                                   <  2:920893/140‑1 (MQ=255)                                                                                         gcCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCAt                                                <  1:665049/140‑1 (MQ=255)                                                                                          ccAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCAtt                                               >  2:1326900/1‑140 (MQ=255)                                                                                          ccAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCAtt                                               >  2:1013524/1‑140 (MQ=255)                                                                                          ccAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCAtt                                               >  2:274880/1‑140 (MQ=255)                                                                                              tCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAag                                           <  2:787365/140‑1 (MQ=255)                                                                                               cATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAaga                                          <  2:809777/140‑1 (MQ=255)                                                                                               cATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAaga                                          <  1:553417/140‑1 (MQ=255)                                                                                                       ccAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGc                                  <  1:1436504/140‑1 (MQ=255)                                                                                                       ccAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGc                                  >  1:1438744/1‑140 (MQ=255)                                                                                                       ccAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGc                                  >  1:575965/1‑140 (MQ=255)                                                                                                          aTGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATc                               <  1:322282/140‑1 (MQ=255)                                                                                                          aTGAAGTGGCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATc                               >  2:225362/1‑140 (MQ=255)                                                                                                                  gCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGaaa                       >  2:39459/1‑140 (MQ=255)                                                                                                                  gCAACCGGGCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGaaa                       >  2:245660/1‑140 (MQ=255)                                                                                                                        gggCTTTTTGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTg                 <  1:39459/140‑1 (MQ=255)                                                                                                                              tttGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCg           >  1:153256/1‑140 (MQ=255)                                                                                                                               ttGGGCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGa          <  2:861135/140‑1 (MQ=255)                                                                                                                                 gggCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACg        >  2:310758/1‑140 (MQ=255)                                                                                                                                 gggCATCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACg        <  2:254345/140‑1 (MQ=255)                                                                                                                                      tCTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACGCCAtt   >  1:1348254/1‑140 (MQ=255)                                                                                                                                       cTGGT‑GGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACGCCATTc  <  2:153256/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                        ATCTGCAAACGCCAGAGTGGGTTGGGGCCGACTGTGCTCGCCGTACTTTATCGCGTCTGTCACCGCGTAAACTCTCCACCATGAAAGCGCCAGTCATTTTTGCCAATGAAGTGGCAACCGGGCTTTTTGGGCATCTGGTGGGGGCGATAGCGGGTGGATCGGTTTATCGTAAATCTACCTTCCTGCTGGATTCGCTGGGTAAACAAATTCTGCCGGACTGGCTGACCATTGAAGAGCATCCGCATCTGCTGAAAGGGCTGGCGTCGACGCCATTC  >  NC_000913/4458602‑4458876

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