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	98.5%	235.5 / NA	66	C44F (TGC→TTC)	ysaB	putative lipoprotein YsaB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (33/31);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATTATTTATATCGGTATAGTGGCTGATATCTTACGAATTTCCTTAAGCCTCTATAAGCTGTAGTGATAAAACTCAGATTATTCC  >  NC_000913/3725344‑3725827                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      aTTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAATACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGAAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTCTAAACGCGCAATTCATAGCTTCCCTGGAACTGTTCGCAGGCGGTGACGTCGATTTTTTGCTCGACGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACCCatcgg                                                                                                                                                                                                                                           >  2:608361‑M1/1‑246 (MQ=255)   ttAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGCCGTCGATTTTTTGCTCGCCGGTGTTATACCGTTTTGCCGCCTGATCCTTGAATAACGCatcggcc                                                                                                                                                                                                                                         >  1:470430‑M1/1‑244 (MQ=255)     aaaGGCAGGGATAAGTATACCCGCTGGATGGAAGATCTACAGTACAAAACGGGCAAAAGAAGGGCTTAACGCATGGCAAGATGCACAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGCATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATCACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctc                                                                                                                                                                                                                                       >  2:308057‑M1/1‑242 (MQ=255)     aaaGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctc                                                                                                                                                                                                                                       <  1:364393‑M1/251‑10 (MQ=255)      aaggaaggGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGTAACTGTTCGAAGGCGGTGACGTCTATTTTTTGCTCGCCGGTGTTATAACGTTGTGCCGCCTGATCCTTGAATAACGCatcggcctcg                                                                                                                                                                                                                                      <  2:520536‑M1/251‑11 (MQ=255)          aaggGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatctgcctcggatt                                                                                                                                                                                                                                  <  2:568667‑M1/251‑15 (MQ=255)          aaggGATCAGTATACCCGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggatt                                                                                                                                                                                                                                  >  1:277420‑M1/1‑237 (MQ=255)          aaggGATCAGTATACCAGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATACGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggatt                                                                                                                                                                                                                                  >  2:384597‑M1/1‑237 (MQ=255)                         ccGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgt                                                                                                                                                                                                                   <  2:173764‑M1/251‑30 (MQ=255)                         ccGCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgt                                                                                                                                                                                                                   <  2:43868‑M1/251‑30 (MQ=255)                           gCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACCGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCCTAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttg                                                                                                                                                                                                                 >  1:290553‑M1/1‑220 (MQ=255)                           gCTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttg                                                                                                                                                                                                                 >  2:45244‑M1/1‑220 (MQ=255)                            cTGGATGGAAGATATACAGTACAAAACGGGAAAAAGCCGGGATTAACGCCCGGAAAGATGCAAAACATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGAATAACCGCGCATCTCATAGATTCCCTGGAACTGTTCGAAGGCGGTGACGTTGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttctcagtgttgt                                                                                                                                                                                                                >  1:279873‑M1/1‑219 (MQ=255)                              ggatggaAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgttt                                                                                                                                                                                                              >  2:450613‑M1/1‑217 (MQ=255)                               gatggaAGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttt                                                                                                                                                                                                             >  1:57793‑M1/1‑216 (MQ=255)                                    aaGATATACAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggc                                                                                                                                                                                                        <  1:24808‑M1/251‑41 (MQ=255)                                            cAGTACAAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCAGCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGTAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctgtcagcc                                                                                                                                                                                                <  2:241810‑M1/251‑49 (MQ=255)                                                        ggAAAAAGCCGGGCTTAACGCATGGACAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTAAGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcc                                                                                                                                                                                    >  2:571219‑M1/1‑191 (MQ=255)                                                           aaaaGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggc                                                                                                                                                                                 >  2:607110‑M1/1‑188 (MQ=255)                                                                cAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatt                                                                                                                                                                            >  2:149819‑M1/1‑183 (MQ=255)                                                                 aGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatat                                                                                                                                                                           <  2:462554‑M1/251‑70 (MQ=255)                                                                 aGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatat                                                                                                                                                                           <  2:42798‑M1/251‑70 (MQ=255)                                                                          cGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaat                                                                                                                                                                  <  1:182590‑M1/251‑79 (MQ=255)                                                                           gCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATaa                                                                                                                                                                                                                                                   <  2:199848/169‑1 (MQ=255)                                                                           gCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATaa                                                                                                                                                                                                                                                   >  1:199848/1‑169 (MQ=255)                                                                           gCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAACTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGGTCCTTGCATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagcccgaagggagtcaggctgatattttgacaat                                                                                                                                                                  >  1:84452‑M1/1‑172 (MQ=255)                                                                             aTGGAAAGATGCAAAAAATGGCCCTCCGCGTCAACAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatccgtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaatatt                                                                                                                                                               >  2:107902‑M1/1‑170 (MQ=255)                                                                               ggAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatcc                                                                                                                                                             <  1:550132‑M1/251‑84 (MQ=255)                                                                                gAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccg                                                                                                                                                            <  1:234120‑M1/251‑85 (MQ=255)                                                                                  aaGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccggg                                                                                                                                                          >  1:168253‑M1/1‑165 (MQ=255)                                                                                            aaaTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattc                                                                                                                                                <  1:149819‑M1/251‑97 (MQ=255)                                                                                            aaaTGGCCATCCGCGTCAAAAGAACAGACAAAATGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattc                                                                                                                                                >  2:577207‑M1/1‑155 (MQ=255)                                                                                              aTGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctc                                                                                                                                              <  2:144055‑M1/251‑99 (MQ=255)                                                                                                ggCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgc                                                                                                                                            >  1:608815‑M1/1‑151 (MQ=255)                                                                                                     tCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACAGCGCATTTCATAGCTTCCATGGAACTGTTCGAAGGCGGTGACGCCGATTTTTTGCTCGACGGTGTTATCACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccata                                                                                                                                        >  1:552461‑M1/1‑146 (MQ=255)                                                                                                            cAAAAGAACAGACACACTGGTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGACGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgccaataatccgggttcgcgattctcgccataacacaaa                                                                                                                                 >  2:413757‑M1/1‑139 (MQ=255)                                                                                                            cAAAAGAACAGACAAACTTCTCTTTAAGGAAGTTATAACCGCGCATTCCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTACTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTGGAATAACGCatcggcctcggatgcacgtttatcagtgttgttttttggctggcagccagaagggagtcaggctgatattttaacaataatccgggttcgcgattctcgccataacaccaaa                                                                                                                                <  2:306356‑M1/251‑113 (MQ=255)                                                                                                            cAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTGCCCTGGAACTGTTCGAATGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTGGCCGCCTGTTCCTGGAATAACGCat                                                                                                                                                                                                                                              <  2:303255‑M1/141‑3 (MQ=255)                                                                                                            cAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGACTTTCATAGCTTCCCTGGAACTGTTCGAAAGCGGTGACGACGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCat                                                                                                                                                                                                                                              >  1:303255‑M1/1‑139 (MQ=255)                                                                                                             aaaaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgccaataatccgggttcgcgattctcgccataacaccaa                                                                                                                                 >  2:296377‑M1/1‑138 (MQ=255)                                                                                                               aaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaaga                                                                                                                              >  1:536551‑M1/1‑136 (MQ=255)                                                                                                                aGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGTGCCTTACATAGCTTCCCTGGAACTGTTCGCAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagcgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaat                                                                                                                            >  1:359060‑M1/1‑135 (MQ=255)                                                                                                                    cAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttagcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatt                                                                                                                        <  2:309119‑M1/251‑121 (MQ=255)                                                                                                                          aaCTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagg                                                                                                                  >  2:231188‑M1/1‑125 (MQ=255)                                                                                                                           aCTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgttgatcagtgttgtttttgggctggcagccagaagggagtcaggcttatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagag                                                                                                                 <  2:454583‑M1/251‑128 (MQ=255)                                                                                                                                 ctTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGACGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCctcggcctcggcttcacgtttatcagtgttgtttttgggctggccgccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatg                                                                                                           >  1:418339‑M1/1‑118 (MQ=255)                                                                                                                                 ctTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTAAAACCGTTAGCAGGCGGTGAAGTCGATTCTTTGCTCGCCGCTGTTATCACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcgcattcacgtttatcagcgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggatcgcgattctcgccataccaccaaagcataatttttagaggtgatg                                                                                                           >  1:390207‑M1/1‑118 (MQ=255)                                                                                                                                   ttACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgat                                                                                                            <  1:474274‑M1/248‑133 (MQ=255)                                                                                                                                   ttACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGGGCCGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgat                                                                                                            <  2:437968‑M1/248‑133 (MQ=255)                                                                                                                                      gggaAGGTATAACCGCGCATTTCATAGCTTCCCTGGATCTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgtcgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttg                                                                                                      <  2:10154‑M1/250‑139 (MQ=255)                                                                                                                                         aaGGTATACCAGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataattttagaggtgatgagttgctt                                                                                                    >  1:494467‑M1/1‑110 (MQ=255)                                                                                                                                            gTATCACCTCGCATTACATAACTTCCCTGGAACTGTTCGAAGTCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACTTTTTGCCGCCGGATGCTTGAATAACGCagcggcctcggattcacgtttatcagttgtgtttttgggctgccagccagaagggagtcaggcggatatttggaaaataatcctggttcgcgatgctcgccataacaccaaagaataatttttagaggtgatgagttgcttagt                                                                                                <  2:188082‑M1/251‑145 (MQ=255)                                                                                                                                              ataACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagtta                                                                                              <  1:413019‑M1/251‑147 (MQ=255)                                                                                                                                                 aCCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacat                                                                                           <  1:439229‑M1/251‑150 (MQ=255)                                                                                                                                                     cgcATTTCATAGCGTCACTGGATCTGTTCGAATGCGGTGACGTCGATTTTTTGCTCGCCGGTTTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaatattcctggttcgcgatgctcgccataacaccaaagaataagtttgagaggtgatgagttgcttagttacttaacg                                                                                       <  2:289159‑M1/251‑154 (MQ=255)                                                                                                                                                       cATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGGCGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATTACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgat                                                                                     <  2:513092‑M1/251‑156 (MQ=255)                                                                                                                                                         tttCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattg                                                                                   >  1:31527‑M1/1‑94 (MQ=255)                                                                                                                                                                                    aGGCGGTGACGTCGATTTTTTGATCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgt                                                           >  1:186626‑M1/1‑67 (MQ=255)                                                                                                                                                                                                 gATTTTGTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggagtcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatatttggacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtg                                            <  2:168253‑M1/250‑197 (MQ=255)                                                                                                                                                                                                        tGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccg                                     <  1:296377‑M1/250‑204 (MQ=255)                                                                                                                                                                                                         gCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgt                                    >  1:162213‑M1/1‑46 (MQ=255)                                                                                                                                                                                                             gCCGGTGTTATCACGTTTTGCCGCCTGATCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaaga                                >  2:401969‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                 gtgtTATAACGTTTTGCCGCCTGAGCCTTGAATAACGCatcggcctcggattaacgtttatcagggttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtaggacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaata                           <  2:269095‑M1/251‑214 (MQ=255)                                                                                                                                                                                                                                   gcttgatCCTTGTAGAACGCagcggcctcggattcacgttgaccagtgttgtttttgggctggcagccagaagggagtcaggctgatatttggaccataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacattttgtagatcacatcaggtcaacgccgtaagaaaatatcttgtgattcagagctc         <  2:424126‑M1/248‑232 (MQ=255)                                                                                                                                                                                                                                         tCCTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcacaa       >  1:125493‑M1/1‑14 (MQ=255)                                                                                                                                                                                                                                          ccTTGAATAACGCatcggcctcggattcacgtttatcagtgttgtttttgggctggcagccagaagggagtcaggctgatattttgacaataatccgggttcgcgattctcgccataacaccaaagaataatttttagaggtgatgagttgcttagttacataacgattgtatgacgaaggcataacatgctgtagatcacatcaggtgaacgccgtaagaaaatatcttgtgattcagatcacaaagatt  <  1:45244‑M1/251‑239 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      ATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGCTTAACGCATGGAAAGATGCAAAAAATGGCCATCCGCGTCAAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATTATTTATATCGGTATAGTGGCTGATATCTTACGAATTTCCTTAAGCCTCTATAAGCTGTAGTGATAAAACTCAGATTATTCC  >  NC_000913/3725344‑3725827

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