Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,725,582	C→A	C44F (TGC→TTC)	ysaB ←	putative lipoprotein YsaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,725,582	0	C	A	100.0%	202.7 / NA	52	C44F (TGC→TTC)	ysaB	putative lipoprotein YsaB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (27/25);  total (27/25)

AAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATTATTTATATCGGTATAGTGGCTGATATCTTACGAATTTCCTTAAGCCTCTATAAGCTGTAGTGATAAAACTCAGATTATTC  >  NC_000913/3725349‑3725826                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     aaggaaggGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTTCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGGTTTGCCGCCTGATCCTTGAATAACGCatcggcctcg                                                                                                                                                                                                                                     <  2:487850‑M1/251‑11 (MQ=255)  aggaaggGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcgg                                                                                                                                                                                                                                    <  1:554087‑M1/251‑12 (MQ=255)             aGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgttt                                                                                                                                                                                                                          >  1:46907‑M1/1‑229 (MQ=255)                      gCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttg                                                                                                                                                                                                                >  2:223491‑M1/1‑220 (MQ=255)                                 gATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctg                                                                                                                                                                                                     >  1:93352‑M1/1‑209 (MQ=255)                                    ntacAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCACTGGAACTGTTCGAAGGCGGGGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcgccctcggattaacgtttatcagtgttgtttttgggctggcc                                                                                                                                                                                                  >  2:260146‑M1/2‑206 (MQ=255)                                            caaaacgGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGGATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGGGACGTCGATTTTTTGCTCGCCGGGGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggtctcggattcacgtttatcagtgttgtttttgggctggcagccaga                                                                                                                                                                                            >  1:59123‑M1/3‑198 (MQ=255)                                              aaaCGGGAAAAAGCAGGGCTTAACGCAGGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATGTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggatgcacgtttatcagtgttgtttttgggctggcagccagaaggg                                                                                                                                                                                        <  1:120283‑M1/251‑56 (MQ=255)                                                   ggAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggc                                                                                                                                                                                                   <  2:512307‑M1/235‑45 (MQ=255)                                                   ggAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggc                                                                                                                                                                                                   >  1:512307‑M1/1‑191 (MQ=255)                                                           cAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgata                                                                                                                                                                           <  2:127035‑M1/251‑69 (MQ=255)                                                            aGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatat                                                                                                                                                                          >  2:537721‑M1/1‑182 (MQ=255)                                                                       cATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaata                                                                                                                                                                >  2:297005‑M1/1‑171 (MQ=255)                                                                            aaaGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgg                                                                                                                                                          >  2:80171‑M1/1‑166 (MQ=255)                                                                            aaaGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgg                                                                                                                                                          >  1:152728‑M1/1‑166 (MQ=255)                                                                                aTGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttc                                                                                                                                                      >  1:117466‑M1/1‑162 (MQ=255)                                                                                  gCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgc                                                                                                                                                    >  2:556106‑M1/1‑160 (MQ=255)                                                                                      aaaaTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgatt                                                                                                                                                <  2:117466‑M1/251‑96 (MQ=255)                                                                                         attgccaTCCGCGTCAAAAGAACAGTCAAACTGCTCTTTACGGAAGGTATAACCTCGCATTTCATAGCTTCTCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagatgggagtcagtctgatattttgacaataatccgggttcgcgattcgc                                                                                                                                             <  2:380601‑M1/248‑99 (MQ=255)                                                                                         aTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctc                                                                                                                                             <  1:223491‑M1/251‑99 (MQ=255)                                                                                         aTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctc                                                                                                                                             <  1:537721‑M1/251‑99 (MQ=255)                                                                                             ccATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgcca                                                                                                                                         <  2:152728‑M1/251‑103 (MQ=255)                                                                                               aTCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccata                                                                                                                                       <  1:398572‑M1/251‑105 (MQ=255)                                                                                                 ccGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataac                                                                                                                                     >  2:370833‑M1/1‑145 (MQ=255)                                                                                                        aaaaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaa                                                                                                                                >  1:36214‑M1/1‑138 (MQ=255)                                                                                                          aaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaa                                                                                                                            <  2:265681‑M1/251‑116 (MQ=255)                                                                                                               cAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccacagaataatt                                                                                                                       >  2:528461‑M1/1‑131 (MQ=255)                                                                                                                        tGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggt                                                                                                              >  1:287281‑M1/1‑122 (MQ=255)                                                                                                                               tgcggAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgag                                                                                                        <  1:422581‑M1/248‑136 (MQ=255)                                                                                                                                  ggAAGGTATAACCGCTCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgc                                                                                                    <  1:536440‑M1/251‑140 (MQ=255)                                                                                                                                     aGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgctt                                                                                                  >  1:341875‑M1/1‑109 (MQ=255)                                                                                                                                              cgcgTATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacata                                                                                         >  1:247260‑M1/1‑100 (MQ=255)                                                                                                                                                       natAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccggggtcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtat                                                                               >  2:246802‑M1/2‑91 (MQ=255)                                                                                                                                                          aGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgac                                                                            <  1:546729‑M1/251‑164 (MQ=255)                                                                                                                                                              cccctGGAACTGTTCGAAGGCGGGGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcgcgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataattttaagaggtgatgagttgcttagttacataacgattgtatgacgaag                                                                        <  2:46907‑M1/250‑168 (MQ=255)                                                                                                                                                                  tGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcat                                                                    >  2:122949‑M1/1‑80 (MQ=255)                                                                                                                                                                  tGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcat                                                                    >  1:251228‑M1/1‑80 (MQ=255)                                                                                                                                                                   ggAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcat                                                                    <  1:333062‑M1/250‑172 (MQ=255)                                                                                                                                                                     aaCTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataac                                                                 <  1:122949‑M1/251‑175 (MQ=255)                                                                                                                                                                               aGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgt                                                          >  1:288146‑M1/1‑67 (MQ=255)                                                                                                                                                                                    gTGACGTCGATTGTTTGCTCGCCGGCGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcac                                                   >  1:390527‑M1/1‑62 (MQ=255)                                                                                                                                                                                       aCGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatca                                               <  1:481942‑M1/251‑193 (MQ=255)                                                                                                                                                                                        cGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatc                                                >  2:60886‑M1/1‑58 (MQ=255)                                                                                                                                                                                             aTTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagccaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagattaatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtga                                          >  2:552682‑M1/1‑53 (MQ=255)                                                                                                                                                                                                           ggTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtggacgccgtaagaaaat                           <  1:297005‑M1/251‑213 (MQ=255)                                                                                                                                                                                                                   aaCGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtg                   <  1:331734‑M1/251‑221 (MQ=255)                                                                                                                                                                                                                      gTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgatt                <  1:556106‑M1/251‑224 (MQ=255)                                                                                                                                                                                                                          tGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcag             <  2:287281‑M1/250‑227 (MQ=255)                                                                                                                                                                                                                               ccTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcac        >  1:63312‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                  gATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgattttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcacaa      >  1:481282‑M1/1‑16 (MQ=255)                                                                                                                                                                                                                                   aTCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcacaaaga   <  2:288146‑M1/251‑237 (MQ=255)                                                                                                                                                                                                                                     ccTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcacaaagat  <  1:449488‑M1/250‑238 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                     AAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATTATTTATATCGGTATAGTGGCTGATATCTTACGAATTTCCTTAAGCCTCTATAAGCTGTAGTGATAAAACTCAGATTATTC  >  NC_000913/3725349‑3725826

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