Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,583,800	C→T	D988N (GAC→AAC)	hsdR ←	type I restriction enzyme EcoKI endonuclease subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,583,800	0	C	T	95.2%	127.9 / NA	43	D988N (GAC→AAC)	hsdR	type I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (21/19);  minor base G (0/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAAC‑‑CACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGAT  >  NC_000913/4583579‑4584043                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    ccATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAAC‑‑CACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTc                                                                                                                                                                                                                         <  2:114595/250‑1 (MQ=255)   aTGCTTTGCGCAGGGAAGATTCCTCAAAGGGCTGGCGGTCAAAC‑‑CACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCAGCGTTGATATCGAGGAAGATCGGCGCATCATTCAgg                                                                                                                                                                                                                      <  2:67384/251‑1 (MQ=255)         tGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAAC‑‑CACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGgttgtt                                                                                                                                                                                                                 <  1:2978/250‑1 (MQ=255)                 gAAGATTCCTCAAAGTGCTGGCGGTCAAAC‑‑CACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGt                                                                                                                                                                                                        >  2:30548/1‑251 (MQ=255)                    gATTCCTCAAAGGGCTGGCGG‑‑AAACCACACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGCTCGCGCGGGCGATTAATAACTGCCTGCAATTCCGGTTGCTCGTTCGGGGAAAGTTGCCCCCACGCGGCAAAGGCTTCAAGGAAATCCTGCGGCGAGTCGTAATCCCCGGAAAGCGATTTTACACTCCCCACTTCATCGTTGATAACGAGGAAGAACGGGGcatccttacggtttttgtagtacgtttt                                                                                                                                                                                                     >  2:214883/1‑227 (MQ=255)                    gATTCCTCAAAGCGCTGGCGGACAAAA‑‑CACTCCTGTAGCTCGGCCAACCCTTTCCGGGTGAGATCGCGCGGGCGCTTAATAACTGCCTGACATGCCGGTTGCGCGTTCGGGGGAAGATGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATACCCGTACAGCGCTTTTACACTCCCCCCTTCATCGTTGATATCGAGGAAGACCGGAGCATCATTCCGGCTGTTGATGTCCGtttt                                                                                                                                                                                                     >  2:118951/1‑251 (MQ=255)                                          cAAAC‑‑CACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATa                                                                                                                                                                               <  1:4036/251‑1 (MQ=255)                                           aaaC‑‑CACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATa                                                                                                                                                                               >  1:217827/1‑250 (MQ=255)                                             aagacACTCCTGTAGCTTGACCAGCCCTTTACGGGGGAGATAGCGCGGGCGATTAATAACTGCCGGCAATGCCGGTTGAGCGATCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTGGAGATCGAGGAAGATCGGCGCATCATTCAGGATGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATaa                                                                                                                                                                              <  2:60831/247‑1 (MQ=255)                                                 cACTCCTGTAGCTCGACCAGCCCTTTACGGGGGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGGCCGTTTTCAGCGTTGCCAGACGGGCGATAAAGcc                                                                                                                                                                          <  2:164672/251‑1 (MQ=255)                                                     ccTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCAGCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAg                                                                                                                                                                        <  1:149659/249‑1 (MQ=255)                                                               gACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTaa                                                                                                                                                            <  2:54800/251‑1 (MQ=255)                                                                                  agaTCGAGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTg                                                                                                                                         >  2:248933/1‑251 (MQ=255)                                                                                             gCGATTAATTACTGCCGGCAATGCCAGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGAAAATCCCGCGGCGGGTCGTGATCACCGTACAGCGCTTTTACACCCCCCAATTCATCGTTGATATCGAGGAAGATAGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGATTTTCCTGTTGGGCGATACAGCCAGGTAATTTCTTAAAAACTTCGCCGCTCCAGTCCGTCCCttttt                                                                                                                              >  2:164067/1‑251 (MQ=255)                                                                                              cGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTc                                                                                                                             >  1:129592/1‑251 (MQ=255)                                                                                                              gcAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCg                                                                                                             >  1:110791/1‑251 (MQ=255)                                                                                                                   gCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGcc                                                                                                        <  1:30548/251‑1 (MQ=255)                                                                                                                         tGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCCCCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTaa                                                                                                   >  1:135880/1‑250 (MQ=255)                                                                                                                               ttCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCAc                                                                                            >  1:219486/1‑251 (MQ=255)                                                                                                                                        cGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAAACACCGTACAGCGATTTTACACACACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGCGGCGAAGCCGTTAAAGTTCAcgcccgcagc                                                                                   >  2:241768/1‑248 (MQ=255)                                                                                                                                         gTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCcgcg                                                                                  >  1:68279/1‑251 (MQ=255)                                                                                                                                             cACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCt                                                                              <  1:139695/251‑1 (MQ=255)                                                                                                                                              cccAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGGGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCCCGCCCGCCGCGTCCTg                                                                             <  2:209833/250‑1 (MQ=255)                                                                                                                                                cAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAAGCACCGTACAGCGATTGTACACGCACCGCTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGGCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGGCAGGCAGTTTGTTAAAGACTGCCTCGCTCCAGTGCGGCCCTTTTGCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGc                                                                           <  1:256038/251‑1 (MQ=255)                                                                                                                                                     aGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGTCTTc                                                                                                                                                     <  2:138677/172‑1 (MQ=255)                                                                                                                                                     aGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGGTTGTTAAAGACTTc                                                                                                                                                     >  1:138677/1‑172 (MQ=255)                                                                                                                                                           aGGCTCCGATGAAACCCTGCGGCTTGTCGTAATCCCCGTACAGCGATTTTACCCTCACCACTGCATCGTTGATATCGAGGAAGATCGGCGCTTCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGGTGGTTAAAGTCTTCGCCGCTCCAGTGCGGCCCTTTTTCCCGCCGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCTCGTCCTGGCATAGCTCATcc                                                                <  2:256804/251‑1 (MQ=255)                                                                                                                                                                 cGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTTAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAgacg                                                           >  1:187712/1‑250 (MQ=255)                                                                                                                                                                   aGGAAATCCTGCGGCGTGTCGTAATCACCGTACATCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGacgacg                                                        <  2:205053/251‑1 (MQ=255)                                                                                                                                                                    ggAAATCCTGCGGCGTGTCGTAATCCCCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAgacgaagc                                                       >  2:154335/1‑248 (MQ=255)                                                                                                                                                                          ccTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTGGAGATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGTCAGTATGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTg                                                  <  1:255970/250‑1 (MQ=255)                                                                                                                                                                               tgagtgtCGTAATCACCGTAGAGAGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCTGTTAACGTTCACGCCCGACGCGTCCTGGCACAGCTCATACAGACGACGCACCTGCTtatc                                            >  2:53292/4‑251 (MQ=255)                                                                                                                                                                                   tgtCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTtatctatc                                        >  1:31179/1‑251 (MQ=255)                                                                                                                                                                                     tCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGt                                      >  1:36590/1‑251 (MQ=255)                                                                                                                                                                                                gTACAGCGATTTTACACTCCCCCCTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGt                           <  2:217827/251‑1 (MQ=255)                                                                                                                                                                                                gTACAGCGATGTTACACCCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGCGCGGCCCTTTTCCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGt                           <  2:193710/251‑1 (MQ=255)                                                                                                                                                                                                        aTTTTACACTCACCACTTCAGCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTGTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGGTTCGCTGCGGTCACGGt                     <  2:219486/249‑1 (MQ=255)                                                                                                                                                                                                          tttACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGCCGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGtt                    >  2:144386/1‑248 (MQ=255)                                                                                                                                                                                                          tttACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGtt                    <  1:144386/248‑1 (MQ=255)                                                                                                                                                                                                           ttACACTCCCCTCGTCATCGTTGATAACGCGGAAGATCGGCTCATCATTCAGTTTGTGGATGTCCGTTTGCAGCTTTTCCCGACGGGCGATAAAGCCAGGCAGTTTTTTAAAGACTTCGGCGCCCCAGTGCGTCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAATTTCCCGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTAGCTATCGTTTCGCTGCGGTCACGTTTaaa                 <  2:214986/250‑1 (MQ=255)                                                                                                                                                                                                            tACACTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAAc                <  1:91162/250‑1 (MQ=255)                                                                                                                                                                                                             acacTCACCACTTCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGt              >  1:59627/1‑251 (MQ=255)                                                                                                                                                                                                                         tCATCGTTGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGCTAAGGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGACGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCCGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCgat  >  2:128880/1‑251 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    CCATGCTTTGCGCAGGGAAGATTCCTCAAAGTGCTGGCGGTCAAAC‑‑CACTCCTGTAGCTCGACCAGCCCTTTACGGGTGAGATCGCGCGGGCGATTAATAACTGCCTGCAATGCCGGTTGCGCGTTCGGGGAACGTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCCGTTTTCAGCTTTTCCAGACGGGCGATAAAGCCAGGCAGTTTGTTAAAGACTTCGGCGCTCCAGTGCGGCCCTTTTTCCCGCAGGCGCGAGGCGAAGCCGTTAAAGTTCACGCCCGCCGCGTCCTGGCATAGCTCATCCAGACGACGCACCTGTTTATCTATCGTTTCGCTGCGGTCACGGTTAAACGTGGCCAGACCGAT  >  NC_000913/4583579‑4584043

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