Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,240	G→A	L21L (TTG→TTA)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,240	0	G	A	98.0%	165.0 / NA	50	L21L (TTG→TTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/25);  minor base C (1/0);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAAC  >  NC_000913/4637994‑4638484                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     aaGCACCGGTCTGGGGTTGGAGATCGTCAGTGAAGTCGCCAGTTTGTGTAAAGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGGGGCGTGCTGGACTCGTTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTATTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTc                                                                                                                                                                                                                                                  >  2:289778/1‑251 (MQ=255)      gcGGTCTGGGGGTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGGTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCCTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATa                                                                                                                                                                                                                                             <  1:39897/251‑1 (MQ=255)       cGGTCCGGGGCTGGAGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCCCATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGCCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCATGATTGGTGCAGTATTGCTGTTGATAATTCCgcta                                                                                                                                                                                                                                             >  2:219067/1‑247 (MQ=255)         gTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATg                                                                                                                                                                                                                                          >  2:12120/1‑251 (MQ=255)             ggggTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGa                                                                                                                                                                                                                                      <  2:40709/251‑1 (MQ=255)                      gTTCGTCAGTGAGGTCGCCCGTTGGTTTAACGGCGAAGTCACGCTGCTCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCCCCGTCACTTCACATAGCGTCAAATTCTTCCCACTTATGCTTCGTTGCCGGCTGCCATTGCAAAGGAGAATACTATGTTGAAATCCCCCTGGTTCTGGAAAATGAGTAGCCTGGTTGGGGCAGTATTGCTGTTGTTAATTCCGATAATGTTGATTCGGCAgg                                                                                                                                                                                                                             <  2:508241/251‑1 (MQ=255)                                 aGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa                                                                                                                                                                                                                   >  1:399758/1‑250 (MQ=255)                                   gTCGCCTGTTTGTGTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCCCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGACTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTTAAATCCCCCCTGTTCTTGAAAATGACTAGTCCGTTTGGTGCAGTATTGCTGTTGTGAATGCCGAGAATGCTGATTCGGCAGGTGATTGTCGAAc                                                                                                                                                                                                                 <  2:232090/250‑1 (MQ=255)                                        ccGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTg                                                                                                                                                                                                           >  1:229624/1‑251 (MQ=255)                                            ttgtttATCGGCGAAGTCACGTTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGGTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTTTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGAtt                                                                                                                                                                                                        <  1:268245/250‑1 (MQ=255)                                                ttAACGGCGAAGTCACGCTGCGCAACTTGCAGGAAGGTGGCGTGCTTGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCg                                                                                                                                                                                                    <  2:300962/250‑1 (MQ=255)                                                             cACGCGGCGCAGCGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACGTCACATATCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGaa                                                                                                                                                                                       <  1:344735/250‑1 (MQ=255)                                                                gCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCCCCGTCACTTCACATAGCTTCAAATTCTGCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGTAGAAGACTATGTTGAAATCCCCCCGGTTCTGGAAAATGACTAGCCTGTTTGTTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTCCCGTAGCGATGTGTAAGATg                                                                                                                                                                                   <  2:434793/251‑1 (MQ=255)                                                                gCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATg                                                                                                                                                                                   >  1:153842/1‑251 (MQ=255)                                                                      cAACGTGCAGGACCGTGGCGTGCTGGCCTCGCTTCGCCTTAACCGTCACTTCACATAGCTTCAAATTCTTCCAACATAGTATTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTc                                                                                                                                                                             >  2:248213/1‑251 (MQ=255)                                                                                    gTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCCGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGc                                                                                                                                                               <  2:198766/251‑1 (MQ=255)                                                                                      ggCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCTGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCgg                                                                                                                                                             <  1:469059/251‑1 (MQ=255)                                                                                              ggCCTGGCTACGACTTCACAGGCAACTCCCATAGCTTCAACTTCTTCCCCCATAGTCTTCGTATACTGCTGCCATTGCAAAGGAGATGACTCTGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGCTGCAGTATTGCTGTTGTTCATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAACGATGTGGAAGATGCGATTCGACAAAGTACCAGCGGGCAGCaa                                                                                                                                                      >  2:474106/1‑250 (MQ=255)                                                                                                ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCaaaaa                                                                                                                                                   <  1:161265/251‑1 (MQ=255)                                                                                                 cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAAc                                                                                                                                                  >  1:225024/1‑251 (MQ=255)                                                                                                    gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCg                                                                                                                                               >  1:357667/1‑251 (MQ=255)                                                                                                            ttCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCg                                                                                                                                       >  1:139221/1‑251 (MQ=255)                                                                                                                        tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAATGAGAATACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGCTTCGGCAGGGGATTGTCTAACGTGCTGCTTAACGTAGCGATGTGGACGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCATTGGGCAGCTCATCGCTAtt                                                                                                                           >  2:379137/1‑251 (MQ=255)                                                                                                                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCt                                                                                                                        >  1:239747/1‑251 (MQ=255)                                                                                                                           caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGGTGAAAGCCCCCCGTTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCGCATCGCTATTCCt                                                                                                                        <  2:42866/251‑1 (MQ=255)                                                                                                                                     aaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGc                                                                                                              >  1:56341/1‑251 (MQ=255)                                                                                                                                         cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGCCTCTGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAAGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCCAAGTACCAGCGGGCCGCCACAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCttt                                                                                                           >  1:261790/1‑250 (MQ=255)                                                                                                                                           tcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTtat                                                                                                         <  1:88640/250‑1 (MQ=255)                                                                                                                                                    gTCTTCGTATCCTTCCTCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCGGTTCTGGAAAATGACTAGCCTGTTTGGTGCATTATTGCTGTTGTTAATTCCGATAATGCGGATTCGGCAGGTGATTGTCTAACGTGCTGATTACCGTAGCTAGGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCAGCGCGATTCCTGTGACCGAGCTTTATACGGTGCAgg                                                                                               <  2:287522/251‑1 (MQ=255)                                                                                                                                                               ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGCTTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATaa                                                                                     >  1:66127/1‑250 (MQ=255)                                                                                                                                                               ctgctgCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATaa                                                                                     >  1:199688/1‑250 (MQ=255)                                                                                                                                                                      cATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTgg                                                                             <  1:170671/251‑1 (MQ=255)                                                                                                                                                                          gCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGc                                                                          <  2:412310/250‑1 (MQ=255)                                                                                                                                                                                   aaggcTATGTTGAAATCCCCCCTGTTCTGGAAAAGGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTTATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAg                                                                <  1:309073/247‑1 (MQ=255)                                                                                                                                                                                   aagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGCTTGTCGACCGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGa                                                                  >  1:197169/1‑249 (MQ=255)                                                                                                                                                                                       ttATGGTGAAATCCCCCCTGTTCTGGAAAATTACGTTCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCTGAGAATGCTGATGCGGAAGGTGATTGTCGAACGTTCTGATTACCGGAGCGATGTGGAAGATGCGATGCGCCAAATTACCAGCGGGCCCCCAAAACTCGTTGGGCCGATCATCGTTATGCCTGTGACCGAGCTTTATACGGTGCAGGAATAGGATAAATCCGTGGAGCGGTAACGAAGtct                                                             <  2:281562/249‑3 (MQ=255)                                                                                                                                                                                         aTGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGAAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTAt                                                          >  1:236322/1‑251 (MQ=255)                                                                                                                                                                                             tGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTAGTTGGGCCGCTCATCGCTATTCCTGAGACCGAGCTTAATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGATTTATACAt                                                      >  1:347902/1‑251 (MQ=255)                                                                                                                                                                                              gAAATCCCCCCTGTTCTTGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGGTTTTAATTCCTATAATGCTTATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCTATGTGGAAGATGCGATTCGTCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGTTATTCCTGTGACCGAGCTTTATACGGTGCTGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCAt                                                      <  1:102086/250‑1 (MQ=255)                                                                                                                                                                                                    cccccTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTAGGCAGGGGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGt                                               >  2:469060/1‑251 (MQ=255)                                                                                                                                                                                                     ccccTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGtt                                              <  1:352027/251‑1 (MQ=255)                                                                                                                                                                                                          gTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTg                                         >  1:136327/1‑251 (MQ=255)                                                                                                                                                                                                            tCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGa                                        <  2:153575/250‑1 (MQ=255)                                                                                                                                                                                                               ggAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGGAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCAGCGCTATTCCTGTGAGCGAGCTGTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTc                                     <  1:342865/250‑1 (MQ=255)                                                                                                                                                                                                                 aaaaTGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGATCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAtt                                  <  2:503393/251‑1 (MQ=255)                                                                                                                                                                                                                           gCCTGTTTGGTGCAGTAGTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCAtggatggttgat                        <  2:213984/251‑1 (MQ=255)                                                                                                                                                                                                                             cTGTTTGGTGCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAACAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGCGCTTTATACGGTGCAGGGAGAGGATAAAACAGGGGAGCGGAAACGAAGTTTTATCCATATTTGGGTACCTGAGTCAttgatggttgctgg                      >  2:254522/1‑251 (MQ=255)                                                                                                                                                                                                                                      gCAGTATTGCTGTTGTTAATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGa              >  1:55972/1‑250 (MQ=255)                                                                                                                                                                                                                                       cAGTATGGCTGTTGTTAATTCCGATAAATCTGATTCGGCAGGTGATTTTCGAACGTGCTGATTACCGTAGCTATGGGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGCGCAGGAAGAGGATAAAACCGTCTAGCGGAAACGAAGTTTTATCCATTTTTGGTTCCCTGAGTCATTGATGGTTGATGGCAATCAGAAc            <  1:492081/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gttgttAATTCCGCTAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAAc  >  2:138565/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     AAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAAC  >  NC_000913/4637994‑4638484

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