Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,406,476	C→T	A187V (GCG→GTG)	accC →	biotin carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,406,476	0	C	T	98.0%	145.8 / NA	49	A187V (GCG→GTG)	accC	biotin carboxylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/19);  minor base A (1/0);  total (30/19)
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.33e-01

CAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGCGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGC  >  NC_000913/3406236‑3406720                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     cAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGAGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGaaa                                                                                                                                                                                                                                            <  1:74200/251‑1 (MQ=255)  aGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGGTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGaa                                                                                                                                                                                                                                             >  1:63865/1‑249 (MQ=255)  aGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGCAGCGaa                                                                                                                                                                                                                                             >  1:370681/1‑249 (MQ=255)              ccTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGc                                                                                                                                                                                                                                >  1:166044/1‑250 (MQ=255)                  aTGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCaccg                                                                                                                                                                                                                            >  2:281837/1‑247 (MQ=255)                    gggCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGat                                                                                                                                                                                                                          >  1:414552/1‑250 (MQ=255)                     ggCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGAGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGatat                                                                                                                                                                                                                        >  1:116786/1‑251 (MQ=255)                       cGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATgg                                                                                                                                                                                                                      <  1:411655/251‑1 (MQ=255)                       cGACAAAGTATCCGAGATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTGTCCGGTGATTATCAAAGCCGCCGGCGGCGCCGGCGGTCGCGGTATGCGCGTAGGGCGCGGCGACGCTGAACTGCCGCAATCCATCACCATGACCCTTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATgg                                                                                                                                                                                                                      <  1:367785/251‑1 (MQ=255)                                    gCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGaa                                                                                                                                                                                                         <  2:276053/251‑1 (MQ=255)                                      aaTCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAAt                                                                                                                                                                                                       <  1:101013/251‑1 (MQ=255)                                           cggcgATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTg                                                                                                                                                                                                  >  1:222737/1‑251 (MQ=255)                                            ggcgATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTATGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAAAGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGAGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATAACTgg                                                                                                                                                                                                 >  1:63462/1‑251 (MQ=255)                                                         gCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATCCCTGGAAAATCCTCGCCa                                                                                                                                                                                    <  2:346499/251‑1 (MQ=255)                                                                 cccTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCAGTGCTAAACGCATTGGTTATCCGGTGAGTATCAAAGCCTCCGGCGGCGGCGGCGGGCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCAGCGCCATGTCCCGGGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTGGcctgggcgaga                                                                                                                                                                            <  2:408477/251‑10 (MQ=255)                                                                                   tGACGGCCCGCTGGACGCCGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGGGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACg                                                                                                                                                          >  2:262644/1‑251 (MQ=255)                                                                                    gccGGCCCGCTGTGCGACGAGTTGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACg                                                                                                                                                          <  2:388025/248‑1 (MQ=255)                                                                                       ggCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAACAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTc                                                                                                                                                       >  1:213125/1‑250 (MQ=255)                                                                                         cccNCTGGGCGACGATATGGATAAAAACCGTGCCATGGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCTAGATTCAGGTACTGGCTGACGGTCAgg                                                                                                                                                    <  2:370681/251‑1 (MQ=255)                                                                                          ccGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAggg                                                                                                                                                   >  2:258663/1‑251 (MQ=255)                                                                                                  cgacgaTATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGct                                                                                                                                            >  1:108636/1‑250 (MQ=255)                                                                                                   gacgaTATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGctat                                                                                                                                          >  1:111595/1‑251 (MQ=255)                                                                                                    accatatGGAAAAAAACCGTGCCCTTGCTAAACACATGGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCCCGTAGTGCGCGGCGACCCTGAACTGGCACAATCCATCTCCATGACCCGTTTTGAAGCGAAAGCTGCTTTCAGCAACGATATGGGTTACATGGAGAAATACCCGGAAAATCCTCGCCACGTCGATATTCAGGTACTGGCTGACGGTCAGGGCAACGctatc                                                                                                                                         <  2:63462/248‑1 (MQ=255)                                                                                                       atatGGATAAAAACCGTGCCATTGCTAAACGCATTGATTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGTGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATAGGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGctatctat                                                                                                                                      <  2:308413/251‑1 (MQ=255)                                                                                                              aaaaaaCCGTGCCATGGCTAAACGAATTGGGTATCCGGTGATTATCAAAGCCTCCTGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGATCTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCGCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTAGCTATCTGGCga                                                                                                                               <  2:255236/250‑2 (MQ=255)                                                                                                               aaaaaCCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCACCACGTCGAGATTCAGGCACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCgg                                                                                                                               >  1:326048/1‑250 (MQ=255)                                                                                                                 aaaCCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAAACCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGATGGTCAGGGCAACGCAATCTATCTGGCGGAAc                                                                                                                            >  2:136234/1‑251 (MQ=255)                                                                                                                  aaCCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACg                                                                                                                           >  2:269672/1‑251 (MQ=255)                                                                                                                           aTTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGCACGTGACTGCTc                                                                                                                  >  1:47244/1‑251 (MQ=255)                                                                                                                                 aaaCGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGc                                                                                                             >  2:95292/1‑250 (MQ=255)                                                                                                                                  aaCGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCaa                                                                                                           <  1:341366/251‑1 (MQ=255)                                                                                                                                   aCGCATTGGTTATCCGGTGATTATCAAAGCCGCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAAc                                                                                                          <  2:47244/251‑1 (MQ=255)                                                                                                                                     gCATTGGTTATCAGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACGGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGCACTGGCTGACGGTCAGGGAAACGCTATCTATCTTGCGGAACGTGACTGCTCCATGCAAcgc                                                                                                        >  1:231131/1‑251 (MQ=255)                                                                                                                                      cATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAAcgcc                                                                                                       >  1:243111/1‑251 (MQ=255)                                                                                                                                          ggTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACgccgcc                                                                                                    >  2:270925/1‑250 (MQ=255)                                                                                                                                            ttATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGccacc                                                                                                 >  1:11893/1‑251 (MQ=255)                                                                                                                                                ccGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGACGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGa                                                                                              >  1:154599/1‑250 (MQ=255)                                                                                                                                                            aaaGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCga                                                                                    >  1:221299/1‑248 (MQ=255)                                                                                                                                                              aGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATCCAACGCCGCCACCAGAAAGTGGACgaaga                                                                                 >  2:344358/1‑249 (MQ=255)                                                                                                                                                                                              gcgTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGt                                               >  2:241847/1‑251 (MQ=255)                                                                                                                                                                                                  aGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGc                                            <  1:14648/250‑1 (MQ=255)                                                                                                                                                                                                                         cacaATCCATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTaaa                    <  2:255606/251‑1 (MQ=255)                                                                                                                                                                                                                                cATCTCCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTgg                                                                                                                                  <  1:404058/134‑1 (MQ=255)                                                                                                                                                                                                                                cATCTCCATGACCCGTGAGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTCCATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTgg                                                                                                                                  >  2:404058/1‑134 (MQ=255)                                                                                                                                                                                                                                   ctcCATGACCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTg          >  1:363706/1‑251 (MQ=255)                                                                                                                                                                                                                                   caccATGACCCGTTTGGAAGGGAAAGCTGTTTTCAGCAACGATAGGGTTTACATGTAGAAATACCTGTAAAAGCCGCGCCACGTCGAGATTCAGGTAGTTGCTGACGGTCAGGGCAACGCTAGCTATCTGTCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCGAGAACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACTTTTCGCTAAAGCGTGTGTTg          <  2:389292/249‑1 (MQ=255)                                                                                                                                                                                                                                        ttaCCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGAAGTCGAAGAAGCGCCAGCAGCGGGCATTACCCCGGCACTGCGTCGCTACATCGGCGCACGTTGAGCTAAAGCGAGTTTTGATATc     >  2:208316/3‑251 (MQ=255)                                                                                                                                                                                                                                          aCCCGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCgg   >  2:181434/1‑251 (MQ=255)                                                                                                                                                                                                                                           cccGTGTGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGc  <  2:111595/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     CAGAAACCATTCGCCTGATGGGCGACAAAGTATCCGCAATCGCGGCGATGAAAAAAGCGGGCGTCCCTTGCGTACCGGGTTCTGACGGCCCGCTGGGCGACGATATGGATAAAAACCGTGCCATTGCTAAACGCATTGGTTATCCGGTGATTATCAAAGCCTCCGGCGGCGGCGGCGGTCGCGGTATGCGCGTAGTGCGCGGCGACGCTGAACTGGCACAATCCATCTCCATGACCCGTGCGGAAGCGAAAGCTGCTTTCAGCAACGATATGGTTTACATGGAGAAATACCTGGAAAATCCTCGCCACGTCGAGATTCAGGTACTGGCTGACGGTCAGGGCAACGCTATCTATCTGGCGGAACGTGACTGCTCCATGCAACGCCGCCACCAGAAAGTGGTCGAAGAAGCGCCAGCACCGGGCATTACCCCGGAACTGCGTCGCTACATCGGCGAACGTTGCGCTAAAGCGTGTGTTGATATCGGC  >  NC_000913/3406236‑3406720

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