Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,891,864	G→A	P12L (CCG→CTG)	cysJ ←	sulfite reductase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,891,864	0	G	A	100.0%	331.5 / NA	94	P12L (CCG→CTG)	cysJ	sulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/56);  total (38/56)

GGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAG  >  NC_000913/2891729‑2892000                                                                                                                                        |                                                                                                                                          ggcgctggcgTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAg                                                                                                                                       <  1:1193154/140‑1 (MQ=255)     ctggcgTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAAc                                                                                                                                   <  2:1054122/140‑1 (MQ=255)        gcgTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAAcgcg                                                                                                                                <  1:10013/140‑1 (MQ=255)         ggttgcggCAAGCGCAGCAGGCTGGTGATTGGGTACGCCCCAGCAATAGCCAGAAACCCAGGCAATCTGGGTGGGAGTTAAATCGCTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCgg                                                                                                                               <  2:522251/133‑1 (MQ=255)         cgTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCgg                                                                                                                               <  2:1276491/140‑1 (MQ=255)         cgTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCgg                                                                                                                               >  1:65061/1‑140 (MQ=255)         cgTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCgg                                                                                                                               <  1:1485653/140‑1 (MQ=255)             gCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑gg                                                                                                                          <  2:1348504/140‑1 (MQ=255)              cTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGt                                                                                                                         <  2:338283/140‑1 (MQ=255)                gCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTgg                                                                                                                       <  2:1257159/140‑1 (MQ=255)                  aaGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGa                                                                                                                     >  2:1185833/1‑140 (MQ=255)                  aaGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGa                                                                                                                     >  1:1423800/1‑140 (MQ=255)                   aGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGAc                                                                                                                    <  2:1306777/140‑1 (MQ=255)                   aGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGAc                                                                                                                    <  2:62626/140‑1 (MQ=255)                   aGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGAc                                                                                                                    <  1:1289269/140‑1 (MQ=255)                     cgcAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCt                                                                                                                  <  2:331461/140‑1 (MQ=255)                        agcagGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCtgtg                                                                                                               <  1:1340958/140‑1 (MQ=255)                        agcagGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCtgtg                                                                                                               >  2:845598/1‑140 (MQ=255)                        agcagGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCtgtg                                                                                                               <  2:1511677/140‑1 (MQ=255)                         gcagGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAAGGGTGGGACCtgtg                                                                                                               >  2:936105/1‑140 (MQ=255)                            ggCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTgtcgt                                                                                                           >  1:67295/1‑140 (MQ=255)                               tgctgATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCAt                                                                                                        >  1:472163/1‑140 (MQ=255)                                gctgATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATg                                                                                                       >  2:1068103/1‑140 (MQ=255)                                gctgATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATg                                                                                                       <  2:667802/140‑1 (MQ=255)                                  tgATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGcg                                                                                                     <  1:704494/140‑1 (MQ=255)                                  tgATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGcg                                                                                                     <  1:937412/140‑1 (MQ=255)                                  tgATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGcg                                                                                                     <  1:1005235/140‑1 (MQ=255)                                  tgATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGcg                                                                                                     <  1:281447/140‑1 (MQ=255)                                       gAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAGGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCGCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCACGCGTCGtt                                                                                                <  2:412464/140‑1 (MQ=255)                                       gAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGtt                                                                                                <  1:1071139/140‑1 (MQ=255)                                          tACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATg                                                                                             <  2:410731/140‑1 (MQ=255)                                          tACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATg                                                                                             <  2:722917/140‑1 (MQ=255)                                           aCGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGt                                                                                            >  1:876292/1‑140 (MQ=255)                                            cGCCCCAGAAATAGCCAGAGACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCTGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGtt                                                                                           <  2:397431/140‑1 (MQ=255)                                             gCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTc                                                                                          >  2:1109869/1‑140 (MQ=255)                                             gCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTc                                                                                          >  2:816173/1‑140 (MQ=255)                                             gCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTc                                                                                          >  2:517847/1‑140 (MQ=255)                                               cccAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCa                                                                                        >  1:619613/1‑140 (MQ=255)                                               cccAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCa                                                                                        <  2:1165401/140‑1 (MQ=255)                                                    aaaTAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAg                                                                                   >  2:1477748/1‑140 (MQ=255)                                                     aaTAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGc                                                                                  <  1:10426/140‑1 (MQ=255)                                                       tAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAGTCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCaa                                                                                <  1:253421/140‑1 (MQ=255)                                                         gCCAGAGACCCAAGCAAGCAGCGCGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAg                                                                              >  2:722947/1‑140 (MQ=255)                                                         gCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAg                                                                              >  1:383671/1‑140 (MQ=255)                                                         gCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAg                                                                              >  1:25097/1‑140 (MQ=255)                                                         gCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAg                                                                              >  1:1282132/1‑140 (MQ=255)                                                            aGAAACCCAGGCAAGCTGGGCGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTg                                                                           >  2:408977/1‑140 (MQ=255)                                                              aaaCCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGtt                                                                         >  1:633496/1‑140 (MQ=255)                                                              aaaCCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGtt                                                                         <  1:193842/140‑1 (MQ=255)                                                                aCCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTc                                                                       <  2:1499760/140‑1 (MQ=255)                                                                  ccAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTg                                                                     <  2:893999/140‑1 (MQ=255)                                                                     ggCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTgcgc                                                                  <  1:308205/140‑1 (MQ=255)                                                                      gCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGcc                                                                 <  2:67295/140‑1 (MQ=255)                                                                       cAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGccc                                                                <  2:810520/140‑1 (MQ=255)                                                                            tgggtgggAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCa                                                           >  2:466223/1‑140 (MQ=255)                                                                              ggtgggAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAgc                                                         <  2:50852/140‑1 (MQ=255)                                                                              ggtgggAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAgc                                                         >  1:1175990/1‑140 (MQ=255)                                                                              ggtgggAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAgc                                                         <  1:764390/140‑1 (MQ=255)                                                                              ggtgggAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAgc                                                         <  1:863399/140‑1 (MQ=255)                                                                                tgggAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAgcgc                                                       >  2:1165246/1‑140 (MQ=255)                                                                                tgggAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAgcgc                                                       >  2:1028189/1‑140 (MQ=255)                                                                                   gAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCa                                                    <  2:592722/140‑1 (MQ=255)                                                                                     gTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATa                                                  >  2:1333400/1‑140 (MQ=255)                                                                                        aaaTCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATaaaa                                               >  2:799187/1‑140 (MQ=255)                                                                                        aaaTCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATaaaa                                               >  2:316374/1‑140 (MQ=255)                                                                                          aTCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACa                                             >  2:150818/1‑140 (MQ=255)                                                                                          aTCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACa                                             >  2:1266803/1‑140 (MQ=255)                                                                                              gTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACagaag                                         >  2:704497/1‑140 (MQ=255)                                                                                              gTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACagaag                                         >  1:460611/1‑140 (MQ=255)                                                                                              gTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACagaag                                         >  2:1080982/1‑140 (MQ=255)                                                                                               tCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACagaaga                                        <  2:619613/140‑1 (MQ=255)                                                                                                cGTGGCCGCCTGAAGGTGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACagaagag                                       <  1:1338433/140‑1 (MQ=255)                                                                                                cGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACagaagag                                       <  1:1477748/140‑1 (MQ=255)                                                                                                     ccgccTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGaaggt                                  <  2:472163/140‑1 (MQ=255)                                                                                                         cTGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGaaggtaagg                              >  2:986222/1‑140 (MQ=255)                                                                                                          tGAAGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGaaggtaaggt                             <  1:891083/140‑1 (MQ=255)                                                                                                             aGGCGTGCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTaa                          <  1:1524574/140‑1 (MQ=255)                                                                                                                   gCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGATAAGGTAAGGTTAACGGGGc                    >  1:1142078/1‑140 (MQ=255)                                                                                                                   gCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGc                    >  2:1071761/1‑140 (MQ=255)                                                                                                                   gCCAGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGc                    <  2:1473882/140‑1 (MQ=255)                                                                                                                      aGTTGCTCCGGGTTCAACAGAAGCAACGCGGAAGGGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCaa                  <  1:936105/140‑1 (MQ=255)                                                                                                                      aGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCaaa                 <  2:201026/140‑1 (MQ=255)                                                                                                                      aGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCaaa                 <  2:1282132/140‑1 (MQ=255)                                                                                                                      aGTTGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCaaa                 <  2:1151546/140‑1 (MQ=255)                                                                                                                        ttGCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACg               >  2:186612/1‑140 (MQ=255)                                                                                                                          gCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGgt             <  2:326170/140‑1 (MQ=255)                                                                                                                          gCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGgt             <  1:1362897/140‑1 (MQ=255)                                                                                                                          gCTCCGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGgt             <  2:1187934/140‑1 (MQ=255)                                                                                                                             ccGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGtgtg          >  2:239209/1‑140 (MQ=255)                                                                                                                              cGGGTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTgg         >  2:875911/1‑140 (MQ=255)                                                                                                                               gggTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGa        <  1:150818/140‑1 (MQ=255)                                                                                                                               gggTTCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGa        <  1:1266803/140‑1 (MQ=255)                                                                                                                                   tCAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTaa    <  1:626759/140‑1 (MQ=255)                                                                                                                                    cAACAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTaaa   >  2:533636/1‑140 (MQ=255)                                                                                                                                     aaCAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAg  <  1:517847/140‑1 (MQ=255)                                                                                                                                     aaCAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAg  <  1:816173/140‑1 (MQ=255)                                                                                                                                     aaCAGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAg  <  1:1109869/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                          GGCGCTGGCGTCGCTGCAAGCGCAGCAGGCTGCTGATTGAGTACGCCCCAGAAATAGCCAGAAACCCAGGCAAGCTGGGTGGGAGTTAAATCGGTCGTGGCCGCCTGAAGGCGTGCCAGTTGCTCCGGGTTCAACGGAAGCAACGCGGAA‑GGTGGGACCTGTGTCGTCATGCGTCGTTATGTTCCAGTAAGCAAAGCTGTTTCTGCGCCCTGTCAGCGCCCATAAAACAGAAGAGAAGGTAAGGTTAACGGGGCAAACGGTGTGGATTAAAG  >  NC_000913/2891729‑2892000

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