Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,401	C→T	G78G (GGC→GGT)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,401	0	C	T	94.5%	214.3 / ‑6.6	73	G78G (GGC→GGT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (40/29);  minor base G (1/3);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACAC  >  NC_000913/14163‑14636                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             aaaaGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGggcggcggc                                                                                                                                                                                                                                  <  1:319034/251‑1 (MQ=255)   aaGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGcggcggcgg                                                                                                                                                                                                                                <  2:107992/251‑1 (MQ=255)    aGATGGCTAAGCAAGATGATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGt                                                                                                                                                                                                                               <  2:322042/251‑1 (MQ=255)         gCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGCGATCAAGGCAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGGGGTATGGGCGGCGGCGGtttttg                                                                                                                                                                                                                          >  2:668959/1‑249 (MQ=255)            aaGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTggcgg                                                                                                                                                                                                                       >  1:460353/1‑251 (MQ=255)               cAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAATCGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGAAGACAAAGAGGCCGAGGAGAAATTTAAAGAGATAAAGGACGCTTCTGCAGTTCTTACCGACTCGCAAAAACTGGCGGCTAACGCTCATAATGGTAATGCTGCGTTTGAGCAAGTAGGTATGGACGGCGGATGTTTTgtcggcgg                                                                                                                                                                                                                    >  2:594888/1‑251 (MQ=255)                aagcttattaCGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTggcggcggg                                                                                                                                                                                                                   >  2:501538/5‑250 (MQ=255)                             gaTTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCg                                                                                                                                                                                                      >  1:206841/1‑251 (MQ=255)                                   aGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGAGGTAttt                                                                                                                                                                                                >  2:146322/1‑251 (MQ=255)                                          cccAAAACAGCGAAAGAGCATGAAACCAAAAAGGCCGACAAACGTCTGGCCATGAAATACCTCCCGGACCGTAACCAGTGTGACAAAGAGGTCGAGTCGAAATTTAAAGCGATCAAGTAATCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGGTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCTTAGACTTCATCGATATTTTTGGTGa                                                                                                                                                                                         <  2:398963/250‑1 (MQ=255)                                             aaaaCAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGt                                                                                                                                                                                      >  1:616821/1‑251 (MQ=255)                                                      gAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCg                                                                                                                                                                              >  1:286495/1‑250 (MQ=255)                                                        agagCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGTCCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATTCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGTCGata                                                                                                                                                                           <  1:243486/251‑1 (MQ=255)                                                        agagCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGTAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGGTATTTTTGGTGACGTTTTCGGCgct                                                                                                                                                                            >  1:240593/1‑248 (MQ=255)                                                        agagCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGACATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGAGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGGTTTCGGCGatt                                                                                                                                                                           >  2:175303/1‑250 (MQ=255)                                                             atGAAATCAAAAAGGCCTCCAAACGCCTGGCCATGAAATACCACCCGGACCGTAATCAGGGTGACAAAGAGGCCGAGTCGAAATTTAAAGAGATCATGGAAGCTTATGAAGGTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTTGTGACGTTTTCGGCGATAttttt                                                                                                                                                                      <  2:600350/250‑1 (MQ=255)                                                                              tacatacgtttgtccATGAAATACCACCCGGACCGTTGCCAGTGTTACAAAGAGGCCGAGTCGAAATTAAAAGAGATCTAGGAAGCTTATGAAGTTCTGACCGCCTCGCAAAAACGTGCGGCATACTATCAGTATGGTCATGCTTCGTTTGTTCACGGGGGTATGTGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCATCGATATTTTGTGTGACGTTTTCGGCGATATTTTTGTCTGCGGACTTGgtc                                                                                                                                                      <  2:319104/237‑1 (MQ=255)                                                                                 aaccgacTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTTGCGGCGGAGCCGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGTCATGgtcgtc                                                                                                                                                   >  2:119651/7‑250 (MQ=255)                                                                                  aaCGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCaa                                                                                                                                                 <  2:206841/251‑1 (MQ=255)                                                                                    cGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACg                                                                                                                                               >  1:538526/1‑251 (MQ=255)                                                                                    cGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATAGTTTTGGCGGCGGACGTGGTCGTCAACg                                                                                                                                               <  1:162009/251‑1 (MQ=255)                                                                                               gAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGcgcgcg                                                                                                                                    >  1:14059/1‑251 (MQ=255)                                                                                                   tACCACCCGGACCGTAACCAGGGAGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGGCTTCAGCGGTATTTTTGGTGGCGGTTTCGGCGATAGTTTTGGCGGCGGCGGGGGTCGCCAACGATCGGCGCGCGGTgg                                                                                                                                >  2:120028/1‑250 (MQ=255)                                                                                                         cgggaacGTAACCAGGGTGACAAAGATACCGAGGCGAAATGTAAAGAGATCAATGATGCTTATGAAGTTCCTACCGACTCGCAAAAACGGACGGCATACGATCAGTATGGTTATGCTGCGTGTGAGCAGGGTGGGATGGGCGGCGGCGGTTGTGGCGGCGGTGCAGACTTCAGCGATATTTTTGGTGACGTGTTCGGCTATGATTTTGGCGGAGGACGTGGTCGTCAACCTGAGGCGGGCGGTGCTGAtga                                                                                                                          <  2:225056/245‑3 (MQ=255)                                                                                                               cGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCta                                                                                                                    <  1:373443/251‑1 (MQ=255)                                                                                                                gTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCtat                                                                                                                   <  1:97021/251‑1 (MQ=255)                                                                                                                  aaCCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTg                                                                                                                                                         <  2:625235/211‑1 (MQ=255)                                                                                                                  aaCCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTg                                                                                                                                                         >  1:625235/1‑211 (MQ=255)                                                                                                                  aaCCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATaa                                                                                                                 <  1:411133/251‑1 (MQ=255)                                                                                                                  aaCCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTT‑CGGCGATATTTTTGGCGGCGGACGTGGTAGTCAAAGAGCGGCGCGCGGTGCTGATTTAAGCTATAAc                                                                                                                >  1:8315/1‑251 (MQ=255)                                                                                                                   aCCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAAc                                                                                                                >  1:360239/1‑251 (MQ=255)                                                                                                                      aGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATg                                                                                                             >  1:393207/1‑251 (MQ=255)                                                                                                                        ggTGACAAAGAGGCCGAGGCGAAATTTAAAGAGAACAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATGTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGa                                                                                                           <  2:616821/251‑1 (MQ=255)                                                                                                                            aCAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGGTGATTTACGCTATAACATGGAGCTc                                                                                                       >  1:279058/1‑251 (MQ=255)                                                                                                                              aaaGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATTTTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGGGCGGTGCTGATTTCGGCTATAACCTgggtcgcac                                                                                                     >  2:311903/1‑244 (MQ=255)                                                                                                                                agagGCCGAGGCGAAATTTAAAGAGATCAAGGAAGTTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCAccc                                                                                                   <  2:66729/251‑1 (MQ=255)                                                                                                                                agagGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCAccc                                                                                                   >  1:539441/1‑251 (MQ=255)                                                                                                                                agagGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTGCTGACCGACTCGCAAAACCGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGGATGGGCGGCGGCGGTTTTGGCGGCGGCGCCGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATTGAGCTCAccc                                                                                                   <  1:630681/251‑1 (MQ=255)                                                                                                                                 gagGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGTTATTTTTGGCGGCGGGCGGGGTCGTCAACTTGCGGCGTGCGGGGCTGATTTACGCTATAACaaggggcccccccc                                                                                                  >  2:132665/1‑238 (MQ=255)                                                                                                                                 gagGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCt                                                                                                  <  1:71381/251‑1 (MQ=255)                                                                                                                                  agGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTc                                                                                                 <  1:269623/251‑1 (MQ=255)                                                                                                                                          gCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAAAGTGCGGCGTACGATCAGTCTGGTCATGCTGCGGTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGGCTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCCACGTGCGGAGGTCGGGGCTGATTTCCACGATAGCCTGGAGCTCAccccggatgaggg                                                                                         >  2:173501/1‑241 (MQ=255)                                                                                                                                                 ttAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCACGGTGGTATGGGCGGCGGCGGATTTGGCGGCGGCGCAGACTTCAGCGCTATTTTTGGTGACGTTTTCGGCGATATTTTTCGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATCACATGGAGGTCACCCTCGAAGAAGCTGTACGt                                                                                  >  1:38342/1‑251 (MQ=255)                                                                                                                                                 ttAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACgg                                                                                  >  1:678971/1‑250 (MQ=255)                                                                                                                                                    aaGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGAACAGTATGGTCAGGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAATAAGCTGTACGTgga                                                                               <  2:417677/251‑2 (MQ=255)                                                                                                                                                       agaTCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGGCTTCAGCGATACTTTTGGCGACGTTTTCTGCGATATTTTTGGCGGCCGACGCGGTCGTCAACGTGCAGCGCGCGGTGCTGATTTACGCTATAACACGGAGCTCACCCTCGAAGAAGCTGTACCTggaggg                                                                            >  1:143083/1‑247 (MQ=255)                                                                                                                                                       agaTCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATAAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGCCGTTTTCGGCGATATTTTTGGCGGCGGATGTGGTCGTAAAAGGGGGGCGGGCGGGGCCGATTTACGAGAAAACATAGAGGTCACCCACGAAGAAGCTGAACGTGGCggg                                                                            >  2:379560/1‑249 (MQ=255)                                                                                                                                                           cAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAc‑c                                                                        >  1:22600/1‑250 (MQ=255)                                                                                                                                                              ggAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAag                                                                    <  1:146322/251‑1 (MQ=255)                                                                                                                                                               gAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGGTTTACGCTTTAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGCC‑Caa                                                                      >  1:449685/1‑248 (MQ=255)                                                                                                                                                                    ttATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCTTGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGAGCAGAATTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGAAGTGGTCGGCAACTGGCGGCGTGCGGGGCTGGTTTTCGGTATAACATGGAGCTCACCCTCGAAGAAGCTGTACAGGGCGTGAC‑CAAagagttac                                                              >  2:350004/1‑247 (MQ=255)                                                                                                                                                                    ttATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGatat                                                                                                                                                                          >  2:93290/1‑144 (MQ=255)                                                                                                                                                                    ttATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGatat                                                                                                                                                                          <  1:93290/144‑1 (MQ=255)                                                                                                                                                                                gACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCACCGTGCGGCGCGCGGTGCTGATTTTCGCTATAAAATGGCGCTCACCCTCGACAAAGCCGTACGGGGCG‑CACGCAAAGCGATCCGCATTcccaccc                                                  >  2:637219/1‑246 (MQ=255)                                                                                                                                                                                  ccGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTg                                                <  1:501538/251‑1 (MQ=255)                                                                                                                                                                                            aaaaCGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGtgtg                                      <  1:350004/251‑1 (MQ=255)                                                                                                                                                                                                cGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGTGCGGCGGCGGTTTTGGCTGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGt                                  <  1:677028/251‑1 (MQ=255)                                                                                                                                                                                                    cGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAACGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGc                              >  1:420864/1‑251 (MQ=255)                                                                                                                                                                                                       cATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCAc                           <  1:175303/251‑1 (MQ=255)                                                                                                                                                                                                         tACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACgg                         >  1:648939/1‑251 (MQ=255)                                                                                                                                                                                                          aCGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACgg                         <  1:311903/250‑1 (MQ=255)                                                                                                                                                                                                             aTCAGTTTGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGTCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTGTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGc                     <  2:539441/251‑1 (MQ=255)                                                                                                                                                                                                                aGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAA‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGt                  >  1:663931/1‑251 (MQ=255)                                                                                                                                                                                                                 gTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTa                                                                                     >  1:575219/1‑184 (MQ=255)                                                                                                                                                                                                                 gTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGGTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCCCCCTCGAAGAAGCTGTa                                                                                     <  2:575219/184‑1 (MQ=255)                                                                                                                                                                                                                     ggTCATGCTGCGTTTGAGCAAGGTGGGATGGGCGGCGGCGGATTCGGCGGCGGCGCAGACTTCAGCGATCTTTTTGGCGACGTTTTCGGCGAAATTCTTGGAGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTTTAACCTGGAGCTCACCCTCGAAGAAGCAGTACGTGGAGTGAC‑CAAAGAGATCAGCATTCCGACTCTCGCAGAGTGTGACGTTTGCAACGGTCGCGGCGCaa              >  1:399004/1‑250 (MQ=255)                                                                                                                                                                                                                         aTGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCACAAcc          >  1:372882/1‑250 (MQ=255)                                                                                                                                                                                                                         aTGCTGCGGTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAAcc          >  1:243108/1‑250 (MQ=255)                                                                                                                                                                                                                           gctgcGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAgg       <  1:120028/251‑1 (MQ=255)                                                                                                                                                                                                                             tgcGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGAGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGGGGCGTGAC‑CAACGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGt      >  1:460719/1‑250 (MQ=255)                                                                                                                                                                                                                              gcGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTac    <  1:674134/251‑1 (MQ=255)                                                                                                                                                                                                                                tggggAGCAAGGGGGGATGGGCGGCGGCGGTTGTGGCGTAGGCTAAGATGTCAGCGTTATTGTTGGTGACGGTTTCTGCTATAGTTTTGGCGGCGGACGTGGTCGTCTACGTGCGGCGCGCGGTGCTGATTTACGCTATAACAAGGATCTAACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGTATTTCTATTCTGGAAGAGGGTGACGTTTGCCACGTTAGCGGTTCAACACCAGTTaca   <  2:280997/246‑1 (MQ=255)                                                                                                                                                                                                                                gTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTacac  >  1:663216/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             AAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGAC‑CAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACAC  >  NC_000913/14163‑14636

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