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%	204.3 / NA	53	C44F (TGC→TTC)	ysaB	putative lipoprotein YsaB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (30/23);  total (30/23)

ACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATTATTTATATCGGTATAGTGGCTGATATCTTACGAATTTCCTTAAGCCTCTATAAGCTGTAGTGATAAAACTCAGATT  >  NC_000913/3725336‑3725822                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 aCGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATaa                                                                                                                                                                                                                                              <  2:335169/251‑1 (MQ=255)       ggATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatc                                                                                                                                                                                                                                        <  2:492170‑M1/251‑4 (MQ=255)                   aggGATCAGTATACCCGCTGGATGGAAGATATCCAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattc                                                                                                                                                                                                                            <  1:70851‑M1/251‑16 (MQ=255)                             ataCCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatc                                                                                                                                                                                                                   <  1:380847‑M1/250‑25 (MQ=255)                               aCCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagt                                                                                                                                                                                                                >  1:232316‑M1/1‑224 (MQ=255)                                 cagctGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGGCATCCGCGTCAAAAGAACAGACAAACTGGTCTTTACGGAAGGTATAACCGCGGATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGAATTTTTGCTCGCCGGTGTTATAACGTTTTGCTGCCTGATCCTTGAATAACGCatcggcctcggattcacgttcatcagagt                                                                                                                                                                                                              >  2:521304‑M1/3‑222 (MQ=255)                                        atggaAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgttttt                                                                                                                                                                                                       >  1:487766‑M1/1‑215 (MQ=255)                                               atatACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctgg                                                                                                                                                                                                <  2:312864‑M1/251‑44 (MQ=255)                                                  taCAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggc                                                                                                                                                                                               >  2:115292‑M1/1‑205 (MQ=255)                                                                ggAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtc                                                                                                                                                                                >  1:364001‑M1/1‑191 (MQ=255)                                                                  aaaaaGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcagg                                                                                                                                                                             >  1:118786‑M1/1‑189 (MQ=255)                                                                    aaaGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggct                                                                                                                                                                           >  2:180723‑M1/1‑187 (MQ=255)                                                                     aaGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctg                                                                                                                                                                          >  2:121180‑M1/1‑186 (MQ=255)                                                                     aaGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctg                                                                                                                                                                          >  2:310888‑M1/1‑186 (MQ=255)                                                                      tgcAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctg                                                                                                                                                                          >  2:45155‑M1/2‑185 (MQ=255)                                                                       gCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgat                                                                                                                                                                        <  2:25867‑M1/251‑68 (MQ=255)                                                                        cAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgat                                                                                                                                                                        >  1:432016‑M1/1‑183 (MQ=255)                                                                               tAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttg                                                                                                                                                                  <  1:488291‑M1/249‑74 (MQ=255)                                                                                        gAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccg                                                                                                                                                       <  1:470559‑M1/251‑85 (MQ=255)                                                                                        gAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccg                                                                                                                                                       <  1:220892‑M1/251‑85 (MQ=255)                                                                                                  aaaaaTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgat                                                                                                                                             >  1:466819‑M1/1‑157 (MQ=255)                                                                                                     aaTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattct                                                                                                                                          >  1:204278‑M1/1‑154 (MQ=255)                                                                                                            aTCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccat                                                                                                                                    >  1:176639‑M1/1‑147 (MQ=255)                                                                                                                 cgTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacacc                                                                                                                              >  2:60899‑M1/1‑142 (MQ=255)                                                                                                                    cAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaa                                                                                                                            <  2:487766‑M1/250‑112 (MQ=255)                                                                                                                     aaaaGAACCGACAAACTGCTCTTTACGGAAGGTCTAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagcaagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaa                                                                                                                            >  1:451566‑M1/1‑138 (MQ=255)                                                                                                                     aaaaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaag                                                                                                                          <  1:413399‑M1/251‑114 (MQ=255)                                                                                                                     aaaaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaa                                                                                                                            >  1:330695‑M1/1‑138 (MQ=255)                                                                                                                         gAACAGCCAAACTTCTCTTTACGGAAGGTATAACCTCGCATTTCATACCTTCCCTGGACCTGTTCGAAGGCAGTGACGTCGATTTTTTTCTCCCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATACCGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaaggtagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaata                                                                                                                      <  1:285614‑M1/251‑118 (MQ=255)                                                                                                                         gAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGAGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCagcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaata                                                                                                                      <  1:521304‑M1/251‑118 (MQ=255)                                                                                                                         gAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTGCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTAGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCTCCTGATCCTTGAATAACGCatcggcctcggattcacttttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaata                                                                                                                      <  1:285595‑M1/251‑118 (MQ=255)                                                                                                                          aaCAGACAAACTGCTCTTTACGGAAGGTATACCCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGAGGTCGATTTTTTGCTCGCCGGTTTTATAACGTGTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataa                                                                                                                     <  1:76332‑M1/251‑119 (MQ=255)                                                                                                                           aCAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataat                                                                                                                    <  2:432016‑M1/251‑120 (MQ=255)                                                                                                                            cAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatt                                                                                                                   >  2:14176‑M1/1‑131 (MQ=255)                                                                                                                                                naagGTATAACCGAGCATTTAATAGATTCCCTGGAATTGTTAGAAGGTGGTGACGTCGATTTTTTGATCGCCGGTGTTATAACGTTTTGCCGCCTGATACTTGAATAACGTatcggcctcggattaacgtttataagtgttgtttttggggtggcagccagaagggagtgaggatgatattttgataataatccgggttcgcgattctcgcgataacaccaaagaataatttttagaggtgatgagttgtt                                                                                               >  2:247907‑M1/2‑111 (MQ=255)                                                                                                                                                             cgcATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacg                                                                                  >  2:448970‑M1/1‑98 (MQ=255)                                                                                                                                                                     aTAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacggttgtatg                                                                          <  1:213583‑M1/251‑162 (MQ=255)                                                                                                                                                                     aTAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatg                                                                          >  2:475209‑M1/1‑90 (MQ=255)                                                                                                                                                                            ccatggaACTGTTCGAAGGCGGTGGCGTCGATTTTTGGCTCGCCGGTGTTATAACAGTTTGCCGCCTGATCCTTGAATAACTCatcggcctgggattcacgtttatcagtgttgtttttggtctggcagccagaagggagtaaggctgatattttgacaataatccggggtcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagtgacataacgattgtatgacgaag                                                                    <  2:66241‑M1/247‑168 (MQ=255)                                                                                                                                                                                 gAACTGTTCGAAGGCGGGGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcata                                                               >  2:121246‑M1/1‑78 (MQ=255)                                                                                                                                                                                           aaGGCGGTGAGGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtag                                                    <  2:376467‑M1/251‑184 (MQ=255)                                                                                                                                                                                               cGGTGACGTCGATTTTGTGCTCCCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcgtcctcggattctcgtttatcagtgttgtttttgggctggcagacagaagcgagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagagca                                                <  2:179295‑M1/251‑188 (MQ=255)                                                                                                                                                                                                agtGAGGTCGATTTTTTGCTCGCCGGTGTTCTTACGTTTTGACGCCTGCTCCTTGAATAACGCatcggcctcggattcacgtttatttgtgttgtttttgggctggcagccagaagggtttcaggctgatattttgacaataaacagggttagcgattctcgccataacaccaaagaataatttttagaagtgaggcgttgattagttacataacaattgtatgacgaaggcataacatgctgtagatcac                                               >  2:382599‑M1/2‑63 (MQ=255)                                                                                                                                                                                                      gTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcagg                                         >  2:440397‑M1/1‑57 (MQ=255)                                                                                                                                                                                                         gATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtg                                       >  1:549713‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                   tCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgta                              >  2:76511‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                   ccgCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaag                            <  1:121246‑M1/250‑208 (MQ=255)                                                                                                                                                                                                                            ttATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatct                   >  1:371816‑M1/1‑35 (MQ=255)                                                                                                                                                                                                                                aaCGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgt                <  1:310888‑M1/250‑220 (MQ=255)                                                                                                                                                                                                                                   gTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgatt            >  2:135034‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                                     tttGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattc           >  1:555876‑M1/1‑26 (MQ=255)                                                                                                                                                                                                                                           gccTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcac    <  1:115292‑M1/251‑232 (MQ=255)                                                                                                                                                                                                                                              tGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcacaa  >  1:407432‑M1/1‑17 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 ACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATTATTTATATCGGTATAGTGGCTGATATCTTACGAATTTCCTTAAGCCTCTATAAGCTGTAGTGATAAAACTCAGATT  >  NC_000913/3725336‑3725822

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