Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,984,193	G→A	*412* (TAG→TAA)	aslB →	putative AslA‑specific sulfatase‑maturating enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,984,193	0	G	A	100.0%	242.4 / NA	74	*412* (TAG→TAA)	aslB	putative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/36);  total (38/36)

AATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGG  >  NC_000913/3984059‑3984328                                                                                                                                       |                                                                                                                                        aaTTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAAt                                                                                                                                     >  1:1242977/1‑139 (MQ=255)           gtgCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGgccgc                                                                                                                           >  2:1331730/1‑139 (MQ=255)            tgCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGgccgcc                                                                                                                          >  1:1549885/1‑139 (MQ=255)              ccGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTg                                                                                                                        <  2:862393/139‑1 (MQ=255)                           ttATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTg                                                                                                                        <  2:248331/126‑1 (MQ=255)                           ttATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTg                                                                                                                        >  1:248331/1‑126 (MQ=255)                             aTTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGgccgcc                                                                                                                          >  1:359283/1‑122 (MQ=255)                             aTTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGgccgcc                                                                                                                          <  2:359283/122‑1 (MQ=255)                                              ccATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGAATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCtt                                                                                        >  2:114347/1‑139 (MQ=255)                                                 tatCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                     <  2:1242977/139‑1 (MQ=255)                                                   tCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  2:793190/139‑1 (MQ=255)                                                   tCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGAATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  1:114347/139‑1 (MQ=255)                                                       aaaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                 >  2:893427/1‑137 (MQ=255)                                                       aaaGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                 <  1:893427/137‑1 (MQ=255)                                                           cAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                      <  2:62095/128‑1 (MQ=255)                                                           cAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                      >  1:62095/1‑128 (MQ=255)                                                             aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTa                                                                                                >  1:1331337/1‑116 (MQ=255)                                                             aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTa                                                                                                <  2:1331337/116‑1 (MQ=255)                                                             aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGc                                                                                           >  2:1255896/1‑121 (MQ=255)                                                             aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGc                                                                                           <  1:1255896/121‑1 (MQ=255)                                                             aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                      <  1:947897/126‑1 (MQ=255)                                                             aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTAt                                                                                      >  2:947897/1‑126 (MQ=255)                                                             aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  2:930602/129‑1 (MQ=255)                                                             aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   >  1:930602/1‑129 (MQ=255)                                                             aTGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGa                                                                         >  2:1421760/1‑139 (MQ=255)                                                              tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                     <  1:398080/126‑1 (MQ=255)                                                              tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                     >  2:398080/1‑126 (MQ=255)                                                              tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  2:342545/128‑1 (MQ=255)                                                              tGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   >  1:342545/1‑128 (MQ=255)                                                               ggCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCg                                                                                                                      >  2:79872/1‑92 (MQ=255)                                                               ggCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCg                                                                                                                      <  1:79872/92‑1 (MQ=255)                                                                     tttGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCaa                                                                 >  1:1644366/1‑139 (MQ=255)                                                                         cTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   >  1:5068/1‑117 (MQ=255)                                                                         cTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  2:5068/117‑1 (MQ=255)                                                                         cTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTc                                                             >  2:195794/1‑139 (MQ=255)                                                                         cTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGG‑TGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  2:1073778/116‑1 (MQ=255)                                                                         cTGCCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGG‑TGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   >  1:1073778/1‑116 (MQ=39)                                                                          tGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATcggcg                                                                                                                                <  1:1503726/71‑1 (MQ=255)                                                                          tGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATcggcg                                                                                                                                >  2:1503726/1‑71 (MQ=255)                                                                             cgcACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATa                                                         <  1:195794/139‑1 (MQ=255)                                                                              gcACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCaa                                                                 >  2:586846/1‑130 (MQ=255)                                                                              gcACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCaa                                                                 <  1:586846/130‑1 (MQ=255)                                                                                aCGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcg                                                                                                     <  2:485539/92‑1 (MQ=255)                                                                                aCGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcg                                                                                                     >  1:485539/1‑92 (MQ=255)                                                                                 cGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTaa                                                                                               <  2:1598511/97‑1 (MQ=255)                                                                                 cGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTaa                                                                                               >  1:1598511/1‑97 (MQ=255)                                                                                   gTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGc                                                   >  2:1107630/1‑139 (MQ=255)                                                                                   gTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTACGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                     >  2:186469/1‑105 (MQ=255)                                                                                   gTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTACGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                     <  1:186469/105‑1 (MQ=255)                                                                                       cccGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCg                                               <  2:1644366/139‑1 (MQ=255)                                                                                                          catgcGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   >  1:870344/1‑84 (MQ=255)                                                                                                          catgcGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  2:870344/84‑1 (MQ=255)                                                                                                          catgcGCATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGa                            <  1:1421760/139‑1 (MQ=255)                                                                                                             gcgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  2:747267/81‑1 (MQ=255)                                                                                                             gcgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   >  1:747267/1‑81 (MQ=255)                                                                                                             gcgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                 >  2:997111/1‑83 (MQ=255)                                                                                                             gcgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                 <  1:997111/83‑1 (MQ=255)                                                                                                             gcgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGc                         >  2:467206/1‑139 (MQ=255)                                                                                                             gcgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGc                         >  1:639044/1‑139 (MQ=255)                                                                                                              cgcATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCg                        >  1:569320/1‑139 (MQ=255)                                                                                                                cATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   >  2:1616864/1‑78 (MQ=255)                                                                                                                cATTTGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  1:1616864/78‑1 (MQ=255)                                                                                                                  tttGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAg                                       <  2:594416/120‑1 (MQ=255)                                                                                                                  tttGCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAg                                       >  1:594416/1‑120 (MQ=255)                                                                                                                     gCTGGTGGTGAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcat                 <  1:1107630/139‑1 (MQ=255)                                                                                                                           ggtgAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                     <  2:153527/65‑1 (MQ=255)                                                                                                                           ggtgAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                     >  1:153527/1‑65 (MQ=255)                                                                                                                             tgAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAg                                                 >  1:99588/1‑99 (MQ=255)                                                                                                                             tgAGTAAGTAAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAg                                                 <  2:99588/99‑1 (MQ=255)                                                                                                                                gtaagtaaAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGt                       <  2:737577/122‑1 (MQ=255)                                                                                                                                gtaagtaaAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGt                       >  1:737577/1‑122 (MQ=255)                                                                                                                                   agtaAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATg                                                                                                       >  1:1282611/1‑39 (MQ=38)                                                                                                                                   agtaAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATg                                                                                                       <  2:1282611/39‑1 (MQ=38)                                                                                                                                   agtaAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   <  1:720477/59‑1 (MQ=255)                                                                                                                                   agtaAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                   >  2:720477/1‑59 (MQ=255)                                                                                                                                    gtaAAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAgg  <  1:467206/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        AATTATTTGTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGG  >  NC_000913/3984059‑3984328

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