Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,285,995	G→A	A459A (GCC→GCT)	acs ←	acetyl‑CoA synthetase (AMP‑forming)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,285,995	0	G	A	94.3%	221.5 / NA	70	A459A (GCC→GCT)	acs	acetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (33/33);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGA  >  NC_000913/4285765‑4286226                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        aGCGCAGACGCAATCGCCGCCGTTCCCAGACGGTGACCGGAGACGTTCAGCACATCGTCCACTCGCCCGGTTATCCAGTAATAGCCATCGTCATCGCGACGCGCGCCGGCGCCGCTGAAATACATAGGTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTCCGCGCCTGACCCGGCCAGGAGTCGTTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTAcc                                                                                                                                                                                                                     <  2:416429/251‑1 (MQ=255)     ccGACTCAATCTCCGCCGTTCCCAGACGGTGACCGGAGACGTTCAGCACATCGTCCACACGCCCGGTTATCCAGTAATAGCCAGCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCg                                                                                                                                                                                                                 <  2:298778/251‑1 (MQ=255)           cAATCTCCGCCGTTCCCAGACGGTGACCGGAGACGTTCAGCACATCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATc                                                                                                                                                                                                            <  1:354121/250‑1 (MQ=255)            aaTCTCCGCCGTTCCCAGACGGTGACCGGAGACGTTCAGCACATCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCg                                                                                                                                                                                                           >  1:270316/1‑250 (MQ=255)               ccccgCCGTTCCCAGACGGTGACCGGAGACGTTCAGCACATCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCa                                                                                                                                                                                                       <  2:427509/247‑1 (MQ=255)                 ccgCCGTTCCCAGACGGTGACCGGAGACGTTCAGCACATCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGGCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAgc                                                                                                                                                                                                     <  1:461889/249‑1 (MQ=255)                         cccAGACGGTGACCGGAGACGTTCAGCACATCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGtt                                                                                                                                                                                              <  1:232850/250‑1 (MQ=255)                             gACGGTGACCGGAGACGTTCAGCACATCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACg                                                                                                                                                                                         >  2:424339/1‑251 (MQ=255)                                     ccGGAGACGTTCAGCACATCGTCCACACGCCCTGTTATCCAGTAATAGCCATCTTCAGCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCgaagaa                                                                                                                                                                                 <  2:323378/251‑1 (MQ=255)                                       ggAGACGTTCAGCACATCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACg                                                                                                                                                                               <  1:377339/251‑1 (MQ=255)                                                     catcgtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCgc                                                                                                                                                                  >  1:300596/3‑249 (MQ=255)                                                         gtcCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTTAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCgg                                                                                                                                                             <  2:408566/251‑1 (MQ=255)                                                             ccacGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCGGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGTCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCttt                                                                                                                                                         <  2:114674/250‑1 (MQ=255)                                                             acacGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGCGTCGGTGATCACCAGGCTGCCTTCGGAAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGGGTTGCCGAACCGGCttt                                                                                                                                                         >  2:362513/1‑251 (MQ=255)                                                              cacGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAGCAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCGCCAGGCTGCCTTCGGTATCTCCTTCCAGCGGGTTACCTTCGTTATCGACCATCGCGGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTc                                                                                                                                                        <  1:223472/251‑1 (MQ=255)                                                               acGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCa                                                                                                                                                       <  1:83563/251‑1 (MQ=255)                                                                      gTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGt                                                                                                                                                >  1:339757/1‑251 (MQ=255)                                                                           ccAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAgcg                                                                                                                                            <  1:449623/250‑1 (MQ=255)                                                                             aGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGcc                                                                                                                                          >  2:120528/1‑250 (MQ=255)                                                                                  tAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGCCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCATTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGc                                                                                                                                      >  2:393415/1‑249 (MQ=255)                                                                                   aGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGc                                                                                                                                   >  1:144347/1‑251 (MQ=255)                                                                                     ccATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCgg                                                                                                                                 <  1:243120/251‑1 (MQ=255)                                                                                       aTCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCggcg                                                                                                                               >  1:163612/1‑249 (MQ=255)                                                                                       aTCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGCTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCggcg                                                                                                                               >  2:341882/1‑249 (MQ=255)                                                                                           tCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGa                                                                                                                            <  1:16633/250‑1 (MQ=255)                                                                                                  gACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGaaa                                                                                                                    <  1:95067/251‑1 (MQ=255)                                                                                                  gACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGGAGCGCCAGGCAGCGGCGTGATCATGaa                                                                                                                     >  1:18577/1‑250 (MQ=255)                                                                                                   aCGCGCGCAGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGACGAACGGCCGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCACCAGGCAGCGGCGTGATCaagcaac                                                                                                                   >  2:355054/1‑245 (MQ=255)                                                                                                        cgcCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAAccgccc                                                                                                              >  1:14835/1‑250 (MQ=255)                                                                                                         gcCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAAccgcc                                                                                                               >  2:447619/1‑249 (MQ=255)                                                                                                              cgccgcTGAAATACATATTTTTGAAGGTGGAGAATTATGTCTGTTCAAAACG‑TTCGGGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCCCCAGTCTGCCTTCGGTGGCGCCTTCCAGCGGGTTACTTTCGTTATCGACCAGCGCCGGTTGCACGCTGAAGAACGGACGTGTTGCCGATCGGTCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCAGTTTc                                                                                                         <  1:214004/250‑1 (MQ=255)                                                                                                                     gAAATACATATTTTTGAAGGTGGAGAAGTAGTTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTCCCTTCGTTATCGACCAGCGCCGTTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTg                                                                                                   <  1:88456/249‑1 (MQ=255)                                                                                                                        aTACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGCTCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGACGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCAAGAAACCGCCAGTTTCGGTCTaccacc                                                                                              >  2:183341/1‑251 (MQ=255)                                                                                                                                tttagAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGGGATCGCCAAACAGCGTACGCGCCTGACCCGGCCGGGAGTCGGTGATCACCAGGCTGCTTTCGGTAGCGCCTGCCAGGGGTTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGGGATCATAAAACCGCCAGTTTCGGTCTGCCACCAGGTATc                                                                                       <  2:100815/246‑1 (MQ=255)                                                                                                                                   ggAAGGTGGAGAAGTAGGCCTGTTCAAACCT‑GTCGTGAGCGCCAAACAGCGTACGCGCCTGCCCCGGCCAGGAGTCGGGGATCACCAGGCTGCCTTCGGTAGCGCTTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGGGATCATTAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGac                                                                                    <  1:373110/249‑1 (MQ=255)                                                                                                                                    gAAGGTGGAGAAGTAGGTCTGT‑CAAAACGGTTCGTGACCGCCAAACGGCGTCCGCCCCTGCCCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTGCGGTAGCGCCTTCCGGCGGGTTCCCTCCGTTACCGACCGGCGCCGGTTGCACGCCGAAGAGCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGacca                                                                                  <  2:163612/251‑1 (MQ=255)                                                                                                                                     aaGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTc                                                                                                                                                   <  1:462037/185‑1 (MQ=255)                                                                                                                                     aaGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCATGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGCAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTc                                                                                                                                                   >  2:462037/1‑185 (MQ=255)                                                                                                                                          ggAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTCCGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGAc                                                                            >  1:192012/1‑251 (MQ=255)                                                                                                                                               aGTAGGTCTGTTAAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGAGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGGATCGACCACCgggaatttt                                                                       >  2:348821/1‑244 (MQ=255)                                                                                                                                                   gttCTGTTCAAACCG‑TTCGGGATCGCCAAACAGCGTTCGCGCCTGCCCCGGCCAGGAGTCGGTGATCCCCAGGCGGCCTTCGTTAGCGCCTTCCGGCGGTTTACCTTCGTTATCTACCAGCGCGGGTTGCACGCCGAAGAAGGGACGTGTTGCCAAACCGGTTTTCAGCTCGGTAGCGCCGGGCAGCGGCGTGATCTTGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGt                                                                   <  1:327803/249‑1 (MQ=255)                                                                                                                                                   ggTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCCGGGATCGACCCCCGGACATTTCTCGt                                                                   >  2:81619/1‑251 (MQ=255)                                                                                                                                                    gTCTGTTCAAAACG‑TTCGTGATCTCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGCTCATGAAACCGCCAGGTTCGGTCTGCCACCAGGTATCGACCACCGGCCATTTCTCGtt                                                                  >  2:119355/1‑251 (MQ=255)                                                                                                                                                        gTTCAAAACG‑TTCGTGTTCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGcc                                                               <  1:257088/250‑1 (MQ=255)                                                                                                                                                           cAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGCCAAGCGCCGGTTGCACGCCGAAGCACGGACGTGTTGCCGCACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGGTTCGGTCTGCCACCAGGAATAGCCCACCGGACATTTCTCGTTGCCGAtt                                                           >  2:367406/1‑251 (MQ=255)                                                                                                                                                           aaaaaCG‑TTCGGGACCGCCAAACAGCGTACGCCCCTGCCCCGGCCAGGAGTCGGGGATCACCCGGCTCCCTTGGGGAGCGCCTTCCCGCGGGTTACCGTCGTTTTCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTTTTGCCGATCCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGGGATAATGAACCCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAtt                                                           <  2:14835/250‑1 (MQ=255)                                                                                                                                                            aaaaCG‑TTCGGGACCGCCAAACAGCGGTCGCGCCTGCCCCGGCCGGGAGCCGGTGATCACCGGGCTGCCTCCGGTGGCGCCTTCGAGCGGGTTCCCTTCGTTATCGACCGGCGCCGGTGGCCCGCCGAAGAAGGGACGGGTGGCGGAACCGGTTTTAGGCTCGGTACCGCTAGGCAGCGGTGTGATCATGAAACCGCCAGTTTCGGTCTGTCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAttt                                                          <  1:185210/251‑1 (MQ=255)                                                                                                                                                             aaaCG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAtttt                                                         >  1:209398/1‑251 (MQ=255)                                                                                                                                                             aaaCG‑TTCGTGACCGCCAAACGGCGTACGCGCCGGCCCCGGCCAGGAGTGGGTGATCCCCGGTCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCCCGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAtttt                                                         <  2:144347/251‑1 (MQ=255)                                                                                                                                                                cG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGACCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCCGGTACCGACCACCGGACATTTCTCGTTGCCGAttttttt                                                      >  2:158356/1‑251 (MQ=255)                                                                                                                                                                 cattCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGTGACGGTGATCACCAGTCTGCCCGCGGTAGCGCCTTCCAGCGGGTTACCTTAGTTATCGACCAGCGCCGGTTACACGCCGAAGAACGCGCGTGTTGCCGATCCAGCTTGCCGCTCGGTAGCGCCCGGCAGCGGCGAGATCATGAAACCGGCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTATCGTTGCCGCttttttt                                                      >  2:372775/3‑251 (MQ=255)                                                                                                                                                                   ttCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTCCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTcc                                                    >  2:182955/1‑251 (MQ=255)                                                                                                                                                                     cGTGATCGCCAACAAGCGCACGCCCCTGCCCCGGCAGGGTGTCGGTGATCACCGGGCTGCCTGCGGTAGCGCCTTCCAGCGGGTTACTTTCGTTATCGACCAGCGCCGGTTGCAAGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATAAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAg                                                  <  1:338567/251‑1 (MQ=255)                                                                                                                                                                      gTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGCACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGt                                                 >  1:334375/1‑251 (MQ=255)                                                                                                                                                                      gTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGCGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGCTTTTTTTCCAGt                                                 >  2:347450/1‑251 (MQ=255)                                                                                                                                                                         aTCGCCAAACCGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGCACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGCAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTAcc                                              >  2:149880/1‑251 (MQ=255)                                                                                                                                                                         aTCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCCCCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGCCCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGTTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTAcc                                              <  1:348821/251‑1 (MQ=255)                                                                                                                                                                            cccaacacgcgtccccgcCTGCCCCGGACAGGAGGCGGGGATCCCAAGGTTGCCTTCGGTGGCGCTTTCCAGCGGGTTCCTTTCGTTACCGACCAGCTCCGGTTGCCCGCCGAAGAAGGGACGTGTCGCCGACCCGTCTTTCAGCTCGGTAGCGCCAGGCGGCGGGGTGATAAGTAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACt                                           <  2:325236/236‑1 (MQ=255)                                                                                                                                                                               aaaCAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTccc                                        >  1:418915/1‑251 (MQ=255)                                                                                                                                                                                     cGTACGCGCCTGACCCGGCCAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGCTTTTTTTCCAGTACCACTCCCCCGCtt                                  >  2:299675/1‑251 (MQ=255)                                                                                                                                                                                     cGTACGCGACTGACCCGGAAAGGAGTCGGTGATCACAAGGCTGCCTTCGGTAGCGACTTCCAGCGGGTTACCTTCGTTATAGACCAGCGCCGGTTGCACGCCAAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGACAGGCAGAGGCGTGATCATGAAAACGCCAGTTTCGGGCTGCCACAAGGTAGCGACCACCGGACAGTTCTCGTTGACGATTTTTTTACAGTACAACTCACACGCtt                                  >  2:384415/1‑251 (MQ=255)                                                                                                                                                                                            cccTGACCCGGCCAGGAGTCGGTGATCCCCAGGCGGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCGGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGtt                           <  1:447582/250‑1 (MQ=255)                                                                                                                                                                                                gacccggctgggctggggggttcaacAGGCTGCCTGCGGTGGCGCTCTCCAGCGGGTGACCGTCGTTAGCGACCCTCCCGGGTTGCGTGCGGAAGACGGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTGGCGCCAGGCAGCGGCGTGATCAGGAAGCCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCTATTTGTTTCCAGTACCACTCCCACGCTTCCGTGTTGAtt                       <  2:192012/226‑1 (MQ=255)                                                                                                                                                                                                     ggCCAGGAGTCGGTTATCCCCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTCCCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTc                  <  1:447619/251‑1 (MQ=255)                                                                                                                                                                                                     ggCCAGGAGTCGGTGATCCCCAGGCTGCCTTCGGTAGCGCCTTCCAGCGTGTTACCTTCGTTATCGCCCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTc                  <  1:358359/251‑1 (MQ=255)                                                                                                                                                                                                        cAGGAGTCGGTGATCACCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTGGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGcc               <  1:184367/251‑1 (MQ=255)                                                                                                                                                                                                            aGTCGGTGATCACCGGGCGGCCTTGGGTAGCGCCTTCCAGCGGGTTACCTTCTTTAGCGACCAGCGCCGGTTGCACGCCGAAGACCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACg           <  1:147195/251‑1 (MQ=255)                                                                                                                                                                                                                      ccccAGGCTGCCTTCGGTAGCGCCTCCCAGCGGTTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCgag   <  1:428566/247‑1 (MQ=255)                                                                                                                                                                                                                      caCCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCCCGGAACCgaga  >  2:209174/2‑250 (MQ=255)                                                                                                                                                                                                                      caCCAGGCTGCCTTCGGTAGCGCCTTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGCGTGATCATGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTGATTGGCTCGCCCACGGAACCgag   >  2:428566/2‑249 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        AGCGCCGACTCAATCTCTGCCGTCCCCAGACGGTGACCGGAGACGTTCAGCACGTCGTCCACACGCCCGGTTATCCAGTAATAGCCATCTTCATCGCGACGCGCGCCGTCGCCGCTGAAATACATATTTTTGAAGGTGGAGAAGTAGGTCTGTTCAAAACG‑TTCGTGATCGCCAAACAGCGTACGCGCCTGACCCGGCCAGGAGTCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGA  >  NC_000913/4285765‑4286226

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