Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	535,337	T→C	A474A (GCT→GCC)	gcl →	glyoxylate carboligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	535,337	0	T	C	98.0%	181.8 / NA	51	A474A (GCT→GCC)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/26);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TCACAAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATT  >  NC_000913/535101‑535577                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 tCACAAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGAt                                                                                                                                                                                                                                    >  1:380509/1‑251 (MQ=255)  cccaatgggCTGCGGCAAAAGTGCGGCATGTTTTAAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGAtt                                                                                                                                                                                                                                   <  2:363183/243‑1 (MQ=255)   acaAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTc                                                                                                                                                                                                                                  >  1:332861/1‑251 (MQ=255)             gCGGCACAAATGCTTCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTcaca                                                                                                                                                                                                                        <  1:77872/251‑1 (MQ=255)              gggCACAAAGGCTGCATGTCTTAAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACaa                                                                                                                                                                                                                       <  2:259374/250‑1 (MQ=255)              cGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTcaca                                                                                                                                                                                                                        >  1:277337/1‑250 (MQ=255)                 cacaAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAAcgc                                                                                                                                                                                                                    >  2:79037/1‑251 (MQ=255)                        gctgcATTTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTg                                                                                                                                                                                                             <  1:182346/251‑1 (MQ=255)                                  tttAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGact                                                                                                                                                                                                    >  1:332928/1‑250 (MQ=255)                                      aaGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGc                                                                                                                                                                                               >  1:267889/1‑251 (MQ=255)                                       aGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCg                                                                                                                                                                                              <  1:79037/251‑1 (MQ=255)                                       aGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCg                                                                                                                                                                                              >  1:100596/1‑251 (MQ=255)                                        gacggccacgtgATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCGGCGCTAGGGGGTTGTGCCGCTGATCCGAAACACAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTGCCTGATTGTAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGGACATCCATGTGCTGGTGAACAACGCCTATCTGGGGCTGATGCGTCAGTCACAGCGCGCTTTTGACATGGACCACTGCGt                                                                                                                                                                                             <  2:281578/240‑1 (MQ=255)                                                     tCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTc                                                                                                                                                                                >  1:136859/1‑251 (MQ=255)                                                           gtgGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAat                                                                                                                                                                          <  1:32636/251‑1 (MQ=255)                                                           gtgGTCAGGCTGGTCCGTTAGGCTAGACGATGCCGGCTGCGCTAGGGGTTTGTGCCGCTGTGCCTAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATGGAAGAGTTAGCGGTTGGCGCGCAGTTCAACATTCCGTACAGCCATGTGCTGGTCAACAACGCTTAGCTGGGGCTGATTCGTCTCTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAat                                                                                                                                                                          <  2:346168/251‑1 (MQ=255)                                                              gTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATc                                                                                                                                                                       <  1:289121/251‑1 (MQ=255)                                                                              aGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCagtgaagtg                                                                                                                                                        <  2:271652/250‑1 (MQ=255)                                                                                 cTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATg                                                                                                                                                    >  1:132733/1‑251 (MQ=255)                                                                                        aTTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACgg                                                                                                                                             >  2:193207/1‑251 (MQ=255)                                                                                          tCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGgtg                                                                                                                                           >  2:127458/1‑251 (MQ=255)                                                                                                tGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGAcc                                                                                                                                     >  2:223837/1‑251 (MQ=255)                                                                                                tGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAAAGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGAGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGAcc                                                                                                                                     >  2:360952/1‑251 (MQ=255)                                                                                                      aGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGt                                                                                                                                 <  1:228923/249‑1 (MQ=255)                                                                                                                      gATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGttt                                                                                                               >  1:317599/1‑251 (MQ=255)                                                                                                                         ccGAAACGCAATGTGGTGTCGATTTCTGGCGACTGTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGGTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTGTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAGGTGAATGGCTACGGTGTTGACCGCGTAAAAGTAGCGGAAGGTTTAgg                                                                                                            <  2:71044/251‑1 (MQ=255)                                                                                                                           gAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGtt                                                                                                          <  2:161542/251‑1 (MQ=255)                                                                                                                             aaCGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGt                                                                                                        <  1:75811/251‑1 (MQ=255)                                                                                                                                 cAATGTGGTGGCGATTTATGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCAAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAg                                                                                                    <  2:264386/251‑1 (MQ=255)                                                                                                                                  aaTGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGc                                                                                                   >  2:169358/1‑251 (MQ=255)                                                                                                                                          ggCGATTACTGGCGACTTTGACTTCTAGTTCCTGATTAAAGAGTTAGCTGTTGGCGCGCAGTTCCACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGGTTAGGTTGTAAAGCTATTCggg                                                                                           >  2:126927/1‑251 (MQ=255)                                                                                                                                                    ggCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCaaa                                                                                   <  2:277337/249‑1 (MQ=255)                                                                                                                                                                 tCCAGTTCCTGATTGAAGAGTTAGCGGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTGGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTgcg                                                                    >  1:2700/1‑251 (MQ=255)                                                                                                                                                                     gTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTAACTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGCAGCTATTGCGcc                                                                  >  2:132080/1‑249 (MQ=255)                                                                                                                                                                        ccTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGcc                                                              <  1:338682/250‑1 (MQ=255)                                                                                                                                                                        ccTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCt                                                             <  1:127458/251‑1 (MQ=255)                                                                                                                                                                         cTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCtt                                                            >  2:129639/1‑251 (MQ=255)                                                                                                                                                                         cTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCtt                                                            >  2:9759/1‑251 (MQ=255)                                                                                                                                                                           gATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCtt                                                            <  1:9759/249‑1 (MQ=255)                                                                                                                                                                              tGAAGATTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGTTTGTAAAGCTATTCGGGTCTTCAAACCGTAAGATATTGCGCCAGCCTTTGaa                                                        <  2:345031/250‑1 (MQ=255)                                                                                                                                                                               gAAGAGTTAGCTGTTGGCGCGCAGTTCAGCATTCCGTTCAGCCATGTGCTGGTCAACAACGCCTAGCTGGGGCTGATTCGTCAGGCACAACGCGCTTTTGACATGGACTACTGCGTGGAACTCGGTTTTGAGAATATCATCTCCAGTGAAGTGAATGGCGACGGTGTTGACGACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGTCTTTGAACa                                                      <  1:94118/251‑1 (MQ=255)                                                                                                                                                                               gAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACa                                                      <  2:99605/251‑1 (MQ=255)                                                                                                                                                                                aaGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAg                                                     <  1:139101/251‑1 (MQ=255)                                                                                                                                                                                             tGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTa                                         >  1:251291/1‑250 (MQ=255)                                                                                                                                                                                                cgcgcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATgg                                     <  2:380509/251‑1 (MQ=255)                                                                                                                                                                                                 gcgcAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGgc                                    >  2:80343/1‑251 (MQ=255)                                                                                                                                                                                                             atgccGTACATACATGTGCTGGACAACAACGCAGATCTGGGGCTGATTCATCAGTCACAACGAGCTTTTGACATGGACTACTGCGTGCAACTCGCTTGCGAGAATATCAACGCCAGTGAAGTGAATGGCTACGGTGTTGAACACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTAAAACCGGAATATAGTGCGACAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATAGGGt                        <  2:216104/248‑1 (MQ=255)                                                                                                                                                                                                                ccGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTAcc                     <  2:100596/251‑1 (MQ=255)                                                                                                                                                                                                                        ccATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTc              <  1:351103/250‑1 (MQ=255)                                                                                                                                                                                                                                 ggTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGtt     >  1:29424/1‑250 (MQ=255)                                                                                                                                                                                                                                    caacaaCGCCTAGCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATAGCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAGCAGGCGAAAGCCTTAATGGCGCAAGATCGGGTACCGGTAGTCGTGGAAGTTAtt  <  1:107820/250‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 TCACAAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATT  >  NC_000913/535101‑535577

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