Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,030	C→T	G236G (GGC→GGT)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,030	0	C	T	96.4%	155.5 / ‑4.0	58	G236G (GGC→GGT)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (22/32);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

CTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGT  >  NC_000913/4617793‑4618270                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 ctcCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCGGCctgc                                                                                                                                                                                                                                     >  1:346753/1‑251 (MQ=255)   ccATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCct                                                                                                                                                                                                                                       <  1:204722/247‑1 (MQ=255)    cATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGCTATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTgg                                                                                                                                                                                                                                  <  2:170057/251‑1 (MQ=255)       cAAAGCGGGTGCGGACTGCATCAAATCCTCTACCGGTAAAGTGGCTGTGAGCGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATTTGGGCGTAGAAAAAACCGGTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCaa                                                                                                                                                                                                                               <  2:301789/251‑1 (MQ=255)       aaaaGCGGGTGCGGACGTCATCAAAACCGTAACCGGTAAAGTGGCTGTGAAAGCGACGCCGGAAAGCGCGAGCATCATGATGGAAGGGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGGACTGCGGAAGAGGCGCAGAAATATCTCGCCATCGCAGAGGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCAAGACTGCTGGCaa                                                                                                                                                                                                                               <  2:229911/250‑1 (MQ=255)            cGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGGGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAAATGTTCGGTGCTGACTGGGCAGTTGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCctg                                                                                                                                                                                                                          >  2:123714/1‑251 (MQ=255)                        tcaccaAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTGTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCt                                                                                                                                                                                                               <  2:432392/246‑1 (MQ=255)                          atcaAAACCTCTACCGGGAAAGGGGCTTTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTATCAAACCAGCGGGCGGCGTGCGTACGGCGGAAGAGGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGTGCAGATGCGCGGCACTACCGCTGTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTggg                                                                                                                                                                                                            <  2:207439/251‑1 (MQ=255)                             aaaaCCTCTGCCGGTAAAGGGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGAGGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCAGCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTgggcca                                                                                                                                                                                                         <  2:471507/251‑4 (MQ=255)                                            aaaGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTaaacg                                                                                                                                                                                          >  1:156207/1‑248 (MQ=255)                                               gTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAgcgc                                                                                                                                                                                       >  1:242092/1‑251 (MQ=255)                                                   cTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGATACAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCCGATGAACTGTTCGGTGCTGACTGGGCAGCTGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGGTGAAAGCGCTGGGTCACGGTGACGGTAAGAgcacaagc                                                                                                                                                                                   >  1:211589/1‑245 (MQ=255)                                                     gtgAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCcagcag                                                                                                                                                                                 <  1:322958/251‑1 (MQ=255)                                                            cgACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGGGACGGTAAGAGCGCCAGCAGctacta                                                                                                                                                                           >  2:292293/1‑250 (MQ=255)                                                                            gggcagcatGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATAGCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGc                                                                                                                                                          <  2:83253/245‑1 (MQ=255)                                                                                  caggatgGAAGTGATCCGTGATATGGGAGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtg                                                                                                                                                    <  2:156207/248‑1 (MQ=255)                                                                                    tgatgGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcg                                                                                                                                                  >  2:75450/1‑251 (MQ=255)                                                                                         gAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTgg                                                                                                                                                                                                             >  1:55154/1‑187 (MQ=255)                                                                                         gAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTgg                                                                                                                                                                                                             <  2:55154/187‑1 (MQ=255)                                                                                                        atGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTAc                                                                                                                               <  1:318213/250‑1 (MQ=255)                                                                                                        atGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACGGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTAc                                                                                                                               <  1:415351/250‑1 (MQ=255)                                                                                                                aGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTAt                                                                                                                      <  1:403661/251‑1 (MQ=255)                                                                                                                     aaaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTAcc                                                                                                                 >  2:369578/1‑251 (MQ=255)                                                                                                                     aaaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTAcc                                                                                                                 >  2:147167/1‑251 (MQ=255)                                                                                                                       aCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCg                                                                                                               >  1:202576/1‑251 (MQ=255)                                                                                                                        ccGTTGGTTTCAAACCCGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGt                                                                                                              >  1:147925/1‑251 (MQ=255)                                                                                                                              gTTTCAAAACAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCACGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGAAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACACGTAGGCCg                                                                                                        >  2:418127/1‑251 (MQ=255)                                                                                                                                    aaCCAGCGGGCGGCGTGCGTACTGCGGTAGATGCGCAGATATATCTCGCCATCGCTGATGGACTGTTCGGTGCTGCCTGGGCTGTTGCGCGTCTCTTCCTCTTTGGTGCTTCCCGCCTGCTGGCGAGCCTGTTGCATGCGCTGGGCCACGGTGACGGTAAGAGCGCCTGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATa                                                                                                    <  1:369578/249‑1 (MQ=255)                                                                                                                                           gggcggCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTGCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGa                                                                                           <  2:242092/251‑1 (MQ=255)                                                                                                                                           gggcggCGTGCGGACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGGTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGa                                                                                           <  2:103533/251‑1 (MQ=255)                                                                                                                                           gggcggCGGGCGGACTGCGGAAGAATCGCAGAAATATCGCCCCAACGCAGTTGATCTGTTCGGGGCTGCCTGGTCAGATCAGCTTCAGTACCGCTGTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGGAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACTCCTGATGCGCGGCGCTTATCAGGCCTACGAGACGTATCTAACCGTAGGCCGGATAAGGCGTAGa                                                                                           <  2:277692/251‑1 (MQ=255)                                                                                                                                              cggcgTACGTACTGCGGAAGATGCGCAGAAATATCGCGGCCTCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCCAGCCTGTTGAAAGCGCTGGGCCACGGGGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCGGCGCTTATCAGGCCTACGAGCCGTATCTACCCGTAGGCCGGATAAGGCggaaacgc                                                                                        >  2:416897/1‑244 (MQ=255)                                                                                                                                               ggcgTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGc                                                                                        >  2:47897/1‑250 (MQ=255)                                                                                                                                                 cgtgcgtACTGCGGAAGATGCTCAGAAATATCTCGCCAACGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATc                                                                                     <  2:129367/251‑1 (MQ=255)                                                                                                                                                        aCTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCaa                                                                               >  1:190817/1‑250 (MQ=255)                                                                                                                                                        aCTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCaa                                                                               >  1:123238/1‑250 (MQ=255)                                                                                                                                                           gCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAg                                                                            <  1:75450/250‑1 (MQ=255)                                                                                                                                                                          aaaTATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCttt                                                            >  1:193440/1‑251 (MQ=255)                                                                                                                                                                               tctcGCCACTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCt                                                        <  2:470521/250‑1 (MQ=255)                                                                                                                                                                                  cGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTc                                                       <  1:15521/248‑1 (MQ=255)                                                                                                                                                                                  cGCCATCGCAGATGAACTGTTCGGTGCGGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCg                                                      <  2:193440/249‑1 (MQ=255)                                                                                                                                                                                    ccATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATa                                                                  >  1:309887/1‑235 (MQ=255)                                                                                                                                                                                    ccATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAg                                                  <  1:47897/251‑1 (MQ=255)                                                                                                                                                                                    ccATCGCAGATGAACTGTTCGGTGCTGACGGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCTCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATa                                                                  <  2:309887/235‑1 (MQ=255)                                                                                                                                                                                        agCAGATGAACTGTTCGGGGCGGACTGGGAAGAGGGGCGGCACTACCGCTTTGGGGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGGAGACGCATACGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGtt                                               <  2:90221/249‑1 (MQ=255)                                                                                                                                                                                          cAGATGAACTGGGCGGGGCTGACTGGGAAGATGCGCGTCACTACCGCGTTGGGGCTTCCAGCCGGATGGCAAGCCTGGTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGAAGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTAcc                                            <  2:252660/251‑1 (MQ=255)                                                                                                                                                                                          cAGATGAACCGTTCGGTGCTGACTGGGCAGATTCTCGTCTCTCCCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCAACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGAGGCTTTACGCCTGAGGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGACGCCTCATACGCTTTTCCTCGGGAGGTTAcc                                            <  2:236254/251‑1 (MQ=255)                                                                                                                                                                                            gATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTCCCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCAAGGTGACGGTAAGAGCGCCAGCAGCTACTCAGTACGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACATATCTACCCGTAGGCCGGATAATGCGTAGACGCATCCGGCACAAGACGCCTCATACTCTTTTCCTCGGGAGGATACCtt                                          >  2:260854/1‑251 (MQ=255)                                                                                                                                                                                             aTGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTg                                         <  1:223738/251‑1 (MQ=255)                                                                                                                                                                                                   tGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTtctc                                   >  2:477041/1‑251 (MQ=255)                                                                                                                                                                                                    gtgcGGTGCTGCATGGGCAGAGGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGATGTTACCTTGTTTCTCg                                  <  2:341694/248‑1 (MQ=255)                                                                                                                                                                                                    gTTCGGTGCTGACTGGGCACATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCg                                  <  2:437563/251‑1 (MQ=255)                                                                                                                                                                                                    gTTCGGGGCTGACGGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACAAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTtct                                    <  2:93269/249‑1 (MQ=255)                                                                                                                                                                                                            cTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGCTAAGGCGTAAACGCATCCGGCAAAAGACGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCCCAAGa                           >  2:376472/1‑250 (MQ=255)                                                                                                                                                                                                                      aGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGt                 >  2:301046/1‑250 (MQ=255)                                                                                                                                                                                                                            gcgTCACTCCCTGTTTGGTGTTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCGGGGCCACGGGGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTTTTTCTCGCACAAGAAATTATTCGTAAAAAAc          <  2:297862/251‑1 (MQ=255)                                                                                                                                                                                                                            gcgTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAAc          <  1:97541/251‑1 (MQ=255)                                                                                                                                                                                                                                    aCCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGt  <  1:477041/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 CTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGT  >  NC_000913/4617793‑4618270

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