Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,282	T→C	A35A (GCT→GCC)	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,282	0	T	C	98.0%	149.2 / NA	50	A35A (GCT→GCC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/27);  minor base A (0/1);  total (22/28)
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

GTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGT  >  NC_000913/4638039‑4638522                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  gtttAACGGCGAAGTCTCGCGGCGCAACGTGCAGGAATGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTTCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATGCGGCAAGGGATTGTCGAACGTGCCTATTAcc                                                                                                                                                                                                                                            <  2:56344/251‑1 (MQ=255) gtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTAc                                                                                                                                                                                                                                             <  1:508561/250‑1 (MQ=255)        ggCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCg                                                                                                                                                                                                                                      >  1:103262/1‑250 (MQ=255)            aaGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGCCTTCACCGTCACGTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTATTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGAtgtg                                                                                                                                                                                                                                 <  1:586035/251‑1 (MQ=255)            aaGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGAtgtg                                                                                                                                                                                                                                 <  2:186844/251‑1 (MQ=255)                         cAACGTGCATGAAGGTGTCTTGCTGTCCTCGCTGTGACTTCACCTTCACTTCACATAGCTGCAAATTCTTACCACCTAGCCTTCGTATTCTGTTGCCATTGCAAAGGAGAAGACTATGTTGAAAGCCCCCCTGTTCTGGAAAATGACCACCTTGGTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTc                                                                                                                                                                                                                    <  2:587465/251‑1 (MQ=255)                                            gTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGcc                                                                                                                                                                                                 >  2:315728/1‑251 (MQ=255)                                             tGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCg                                                                                                                                                                                                >  1:294568/1‑251 (MQ=255)                                                       gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCg                                                                                                                                                                                      <  1:587711/251‑1 (MQ=255)                                                        cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGCTTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGCAGATGCGATTCGCCAAAGTACCAGCGGGCCGAAAAAACTCGt                                                                                                                                                                                     >  2:631852/1‑251 (MQ=255)                                                                 cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGTCCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCt                                                                                                                                                                            <  1:558172/251‑1 (MQ=255)                                                                                 aGCTTCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCtgtg                                                                                                                                                            <  1:598509/251‑1 (MQ=255)                                                                                     tCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTAc                                                                                                                                                                                                          <  2:220687/201‑1 (MQ=255)                                                                                     tCAAATTCTTCCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTAc                                                                                                                                                                                                          >  1:220687/1‑201 (MQ=255)                                                                                       aaaTTCTTCCCACATATCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATTTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAATCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAg                                                                                                                                                      <  2:418704/251‑1 (MQ=255)                                                                                           tcttcCCACATAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCtt                                                                                                                                                   <  1:504562/250‑1 (MQ=255)                                                                                                    aTAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGc                                                                                                                                         >  1:164429/1‑251 (MQ=255)                                                                                                    aTAGCCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTGGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGTAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGc                                                                                                                                         <  1:84071/251‑1 (MQ=255)                                                                                                                 ccTGCTGACATTGCAAAGGAGAAGTCTATGTTGAAATCCTCCCTGTTCTGGTACATGTCCCCCTTGTTTGGTGCAGTATTGCTGTTGTTTATTCCGATAAAGCTGACTCGGCAAGGGATTGTCGAACGTGGCGATTATCTTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATTCGGTGTAGGAAGAGGATa                                                                                                                             <  2:524426/250‑1 (MQ=255)                                                                                                                       gcCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAAcc                                                                                                                        <  1:315728/249‑1 (MQ=255)                                                                                                                          atagCAAAGGAGAATATTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTTGTGCAGTATTGCTGTTGGTGATTCCGATAATGTTGATTCGGCAAGTGATTGGCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAATCCGTGGa                                                                                                                   <  2:123975/248‑1 (MQ=255)                                                                                                                                  ggAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAAc                                                                                                           >  1:420196/1‑251 (MQ=255)                                                                                                                                  ggAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAAc                                                                                                           <  1:45422/251‑1 (MQ=255)                                                                                                                                    agaagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAACACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAACCCGTGGAg                                                                                                                  >  2:58300/1‑242 (MQ=255)                                                                                                                                    agaagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAg                                                                                                                  <  1:58300/242‑1 (MQ=255)                                                                                                                                      aagaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑a                                                                                                        >  1:232798/1‑249 (MQ=255)                                                                                                                                       agaCTATGTTGAAATCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGt                                                                                                     <  1:220653/251‑1 (MQ=255)                                                                                                                                                    aTCCCCCCTGTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCAttttt                                                                                        <  1:527712/251‑1 (MQ=255)                                                                                                                                                         cccgtttcgggaaaatgaccccattGTTTGGTGTAGTATTGATGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGATCGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGTGGGCCGCAAAAACTCGTTGGTCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATTCATTTTTGGTTa                                                                                   <  2:535481/227‑1 (MQ=255)                                                                                                                                                            tGTTCTGGAAAATGAACACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCt                                                                                >  1:352183/1‑251 (MQ=255)                                                                                                                                                             gTTCTGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTg                                                                               >  1:536789/1‑251 (MQ=255)                                                                                                                                                                 tGGAAAATGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGcc                                                                                                                                                                                                 >  1:445684/1‑134 (MQ=255)                                                                                                                                                                 tGGAAAATGACCACCGTGTGTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATGCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCATCGGGcc                                                                                                                                                                                                 <  2:445684/134‑1 (MQ=255)                                                                                                                                                                   gAAAATGACCACCTTGTTTGGAGCAGAATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCAAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACG‑AGGTTTTATCCATTATTGGTTACCTggggcat                                                                         >  2:555033/1‑245 (MQ=255)                                                                                                                                                                    aaaaTGACCACCTTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑ACGTTTTATCCATTTTTGGTTACCTGAGTCAtt                                                                        >  1:673996/1‑251 (MQ=255)                                                                                                                                                                                       gTGCAGTTTTGCTGTCGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGGGCCGATGACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACGCGTTGGGACGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGGTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAg                                                     <  2:181873/251‑1 (MQ=255)                                                                                                                                                                                       gTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCACGTGATAGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGAGTCGCCAAAGTACCATCGGGCCGCAAAAACTCGATGTGCCACTTCTCGCTATTACTGTGACCGAGCCTTATACGGTGCAGGAAGAGGATAAAACCGAGGCGCGGAC‑CGTAAGTTTTAACCATTCTTGGTTACCTGAGTCATTGATGGTTGCTGGCCATCa                                                      >  2:191940/1‑250 (MQ=255)                                                                                                                                                                                                  cTGTTGTTGATTCCGATAATGATGATTCGGCAAGTGATTGTCGAACGTGCCGCTTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTAAGTCATTGACGGTTGATAGAAATCAGAACGTGga                                             >  2:453262/1‑248 (MQ=255)                                                                                                                                                                                                        ttgATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTTCCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGaa                                                                                                             <  2:34448/179‑1 (MQ=255)                                                                                                                                                                                                        ttgATTCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGaa                                                                                                             >  1:34448/1‑179 (MQ=255)                                                                                                                                                                                                            ttCCGATAATGCTGATTCGGCCAGTGATTGTCGAACGTGCCGATCACCGTAGCGATGTGGAAGATGCGCTTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGATTATCGATATTCCTGTGACCGCGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AACTTTTATCCCTTTTTGGTTACCTGCGACATTGGTGGTTGATGGCATTCAGTACGTGGAAGAACGCAAGAt                                 >  2:38523/1‑250 (MQ=255)                                                                                                                                                                                                             tCCGATAATGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAACGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGCGCTTTATACGGTGCAGGAAGAGGATAAAAGCGTGGAGCGGAAACG‑AAGTTTTATCCATTTGTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGACGAACGCaggaa                                 >  2:563015/1‑245 (MQ=255)                                                                                                                                                                                                                   aaTGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGACGAACGCAAGATAGGGAttt                         >  1:102153/1‑251 (MQ=255)                                                                                                                                                                                                                   aaTGCTGATTCGGAAAGTGATTGTCGAACGTGCCGAGTACCGTAGAGATGGGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGAttt                         <  2:547763/251‑1 (MQ=255)                                                                                                                                                                                                                    atgctgagtgtgcaagTGATTGTCGACTGTGCCGAGTTCCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACGCGTTGGTCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTGTGGTTACCTGAGTCATTGATGGTTGATGGCAAGCAGATCGTGGAAGAACGCAAGATAGGGATTta                        <  2:386117/236‑1 (MQ=255)                                                                                                                                                                                                                     tGCTGATTCGGCAAGTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTtat                       >  1:201852/1‑251 (MQ=255)                                                                                                                                                                                                                          aTTCGGCAAGTGATTGTCGAACGTGCAGATTACCGTAGCGATGTGGAAGATGCGATTCTCCAAAGTACCAGGGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGCTTTATTCCgg                  <  2:562475/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTGATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑CAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCCAGATAGGGTTTTATACCGGTCAGGTCTg         >  2:262900/1‑251 (MQ=255)                                                                                                                                                                                                                                     gATTGTCGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTgg        <  1:62712/250‑1 (MQ=255)                                                                                                                                                                                                                                           cGAACGTGCCGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGt  <  1:644496/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  GTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCCCCCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACG‑AAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGT  >  NC_000913/4638039‑4638522

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: ‑