Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,420	T→A	A229A (GCT→GCA)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,420	0	T	A	98.2%	174.5 / NA	58	A229A (GCT→GCA)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (27/29);  minor base C (1/0);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGA‑GACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG  >  NC_000913/4177‑4665                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        aCCGAATGTGAAAGTGGTTATTCTTTATTCACGAGGCAAAATCAGTACACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAATTGTCGCCATCGACGGCGATTTCGATGCCTGTCAGGCGATGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTGGCTACTACTTTGAAGCAGTTgcg                                                                                                                                                                                                                                                   <  2:277547/250‑1 (MQ=255) aCCGAAGGTGAAAGTGGTTATCCTCTATCATCGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTgcgc                                                                                                                                                                                                                                                  <  2:8070/251‑1 (MQ=255)  ccGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCa                                                                                                                                                                                                                                                 <  1:456669/251‑1 (MQ=255)      aTGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTg                                                                                                                                                                                                                                             >  1:23299/1‑251 (MQ=255)      aTGTGAAAGTGGTTATCATATATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTCCATTGGGCGGCAATATCGAAACTGTTGCCATCGCCGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAGCATCAGACGTTTGCTGGCGCCGATTTTCTACTACTTTGAAGAAGTTGCGCAGcgg                                                                                                                                                                                                                                             >  2:520112/1‑249 (MQ=255)       tgtgAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGAAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGCAGCAGTTGCGCAGCTGc                                                                                                                                                                                                                                            >  2:412990/1‑251 (MQ=255)            aaGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAg                                                                                                                                                                                                                                       <  1:414830/251‑1 (MQ=255)             aGTGGTTAACCTCTATCCACGAGGCCAAACCAGTCGAATGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAACCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGTTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGGGGCGATGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAGGCCGATGCGCAGCTACCGAAgg                                                                                                                                                                                                                                      >  2:379601/1‑251 (MQ=255)                    aTCCTCTATCCACGAGGCTAAATAAGTTCACTTCAAGAAAAACTTTACTGTCCATTGGTAGTCAATATCGAAACTGTTGCCATCGACGGCGATTTCTATGCCTGTCAGTCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAATCTCGGCTAACTCGATTAACATCAGCGGGATGCTGGCTCAGATTTCCTACTACTTTGAAGCTGTTTCGCAGCTGCCGCAGGT‑GActcg                                                                                                                                                                                                                              <  2:409211/251‑4 (MQ=255)                         ctATCCACGAGGCAAAATCATTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCTATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTGGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCaa                                                                                                                                                                                                                           <  2:585389/249‑1 (MQ=255)                         ctATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAAcc                                                                                                                                                                                                                         >  2:213540/1‑251 (MQ=255)                                    gCAAAATCAGTCCACTGCAAGAAAAACTTTTCTGTACATTGGGCTGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGGAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAATCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGtc                                                                                                                                                                                                              <  2:357277/251‑1 (MQ=255)                                    gCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACGTCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGtc                                                                                                                                                                                                              >  1:32305/1‑251 (MQ=255)                                       aaagcAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCg                                                                                                                                                                                                           <  1:604965/247‑1 (MQ=255)                                        aaTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGACACTGATGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCgg                                                                                                                                                                                                          >  2:149560/1‑251 (MQ=255)                                         aTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCCGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCggg                                                                                                                                                                                                         >  2:431329/1‑250 (MQ=255)                                                cACTGCAAGAAAAACTGTTCTGTTCAGTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGAGGCCTGTCAGGCGCTGTTGAAGCAGGCGTTTGATGATGTAGAGCTGAAAGTGGCGCTGGTGTTAAACGCAGCTAACTCGATGAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAATCAGTTGCCCAGCTGTCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGc                                                                                                                                                                                                  <  2:277319/251‑1 (MQ=255)                                                    gCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCATCCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                                                                                              <  2:526402/251‑1 (MQ=255)                                                     cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                                                                                              >  1:419795/1‑250 (MQ=255)                                                        gAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACtt                                                                                                                                                                                          >  2:330529/1‑251 (MQ=255)                                                           aaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAATCGGAATCTTCgg                                                                                                                                                                                       <  1:499889/251‑1 (MQ=255)                                                           aaaCTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCgg                                                                                                                                                                                       >  2:105044/1‑251 (MQ=255)                                                                 ttCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGAttt                                                                                                                                                                                 >  1:315340/1‑251 (MQ=255)                                                                   cTGTCCATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTg                                                                                                                                                                                >  2:129887/1‑250 (MQ=255)                                                                           tGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCgggg                                                                                                                                                                       >  2:511917/1‑250 (MQ=255)                                                                                   aaTATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTggg                                                                                                                                                               >  1:15404/1‑250 (MQ=255)                                                                                         gAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCggaggc                                                                                                                                                         >  2:174817/1‑246 (MQ=255)                                                                                               gtggcctacGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTGGATGATGAAGAACTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGAGTTGCTGCTACTGTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACGTCGGCGAGTTGACGGCGGGTCTGCTGGCGATGTCACTc                                                                                                                                                     <  2:547897/241‑1 (MQ=255)                                                                                               gTTGCCATCGACGGCTATTTAGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACGCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACGACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGTCACTCg                                                                                                                                                    <  2:451502/250‑1 (MQ=255)                                                                                                  gtcaacGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGGTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCACCCAGCTGGTTGTCTCGGTGCCAAGCGGAACCTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCt                                                                                                                                                <  2:290957/246‑1 (MQ=255)                                                                                                  gCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCgggggtgccggcggagtcccccgggct                                                                                                                                                >  2:231244/1‑227 (MQ=255)                                                                                                        gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGt                                                                                                                                          <  1:171426/251‑1 (MQ=255)                                                                                                              gATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACCCGGTCTGCCGGGGAAACg                                                                                                                                    >  1:172013/1‑251 (MQ=255)                                                                                                                  tCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGCGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGCAGTCACTCCGTCTGCCGGTGaaaaggttt                                                                                                                                >  2:14570/1‑245 (MQ=255)                                                                                                                   cGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGtttt                                                                                                                                <  1:529948/250‑1 (MQ=255)                                                                                                                                     ggTGAAGCCGGCGTTTGATGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTTTTTCTGCGGCCCACGTgg                                                                                                             >  1:200014/1‑250 (MQ=255)                                                                                                                                        gAAGCAGGTGGTTGAAGATGAAGAGCTGAAAGTGGCGCTGGGGTTAAATGCAGCTAACTCGATTAACCTCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGatcg                                                                                                          <  2:233944/251‑4 (MQ=255)                                                                                                                                                      atgatgaAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCCGGCGAAAGCACACGGTCGGCCGGTGGAACGTTTTATTGCTGCGACCCAAGGGGACGCAACCGTGCCACGt                                                                                            >  1:404046/1‑251 (MQ=255)                                                                                                                                                          tgaAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGCAGTCACCCGGTCTGCCGGGGTAACGTGTTTTTTCTGCGCCCAACGTGAACGAGAACGTGCCACggtttc                                                                                        >  2:506223/1‑246 (MQ=255)                                                                                                                                                            aagagCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAATTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTg                                                                                      <  1:563956/251‑1 (MQ=255)                                                                                                                                                               aaCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACGCGATTAACATCAGCCATTTGATGGCGCAGATTTGCTACTATTTTTAAGCAGTTGCGCAGCTGCCGCAGGAAGA‑GCGCAACCAGCTTGTTGTCTCGGTGCCAAGCGGAAACTTCGGCAATTGGACGGCGGGTCTGCTGTCGACGGCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATGCCTTGCCGCGTTTCCTGcgc                                                                                   <  2:7764/251‑3 (MQ=255)                                                                                                                                                                gcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCacg                                                                                  <  1:14570/250‑1 (MQ=255)                                                                                                                                                                   gAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGGTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCacgac                                                                                <  1:104435/250‑1 (MQ=255)                                                                                                                                                                         ggCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAg                                                                          <  1:129887/250‑1 (MQ=255)                                                                                                                                                                                 ggTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGGCACTCGGTCTGCCGGGGGAACGGTTTTTTGCTGCGACCAACGCGAAAGATACCGTGCCACGCTTTCTGCCAGACGGTCAGGGGGCA‑ccc                                                                >  2:430617/1‑251 (MQ=255)                                                                                                                                                                                          cagCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACt                                                       <  1:105044/249‑1 (MQ=255)                                                                                                                                                                                            gCTAACTCGATTAACATCAGCCGTTTGCGGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTTTCGGTGCCAAGCTGAAACTTCGGCTATTTGACGGCGGGTCTTCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATTCCGTGGCACGTTTCCTGCACGACGGTCAGTGTTCT‑CCCAAAGCGACTCa                                                     <  1:538589/251‑1 (MQ=255)                                                                                                                                                                                                aCTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACTCAGGCg                                                 <  1:337115/251‑1 (MQ=255)                                                                                                                                                                                                   cGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTTGGGGAGTCACTCGGTCTGCCGGGGGAACGTTTTTTTTCTGCGACCCACGTGAACGATACCGTGCCACGTTTTCTGCACGACGGGCCGTGGTCA‑CCCCAAGCGACCCAGGCGACg                                              >  1:2605/1‑251 (MQ=255)                                                                                                                                                                                                        aaCATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTTCTGGCGAAGTACCTCGGTCTGCCGGGGCAACGTGTTATTGGTGCGACCAACGGGAACGAAACCGTGCCCCCTTTCCTGCACGGCGGGCCGTGGTCAACCCAAA‑CGACTCAGGccgcggtatc                                         >  1:364483/1‑242 (MQ=255)                                                                                                                                                                                                        aaCATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACTCAGGCGACGTTAt                                          <  1:290014/250‑1 (MQ=255)                                                                                                                                                                                                         aCATCAGCCGTTTGCTGGCGCAGATTTGCTACTACGTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCGGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACTCAGGCGACGTTATcc                                        <  2:295962/251‑1 (MQ=255)                                                                                                                                                                                                                 cGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACTCGGGCTGCCGGTGGAACGTTTTATTGCTGCGACCAACGGGAACGATACCCTGCCACGCGTCCTGTCCGACGGGCCGTGGTCA‑CCCCAAACGACTCAGGCCACGGTATCCAAagcg                                  >  2:313235/1‑249 (MQ=255)                                                                                                                                                                                                                         ggCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGTTCA‑CCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGt                          <  1:520112/249‑1 (MQ=255)                                                                                                                                                                                                                                    gCTACTAATTTGAAGTAGTTGCGCAGCTGCCGCAGTA‑AGCGCGCAACCAGCGGGTTGTCGCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAATCGACTCATGCGCCGTTATCCAGCGCGATGTACGTTAGTCAGCCGaac             <  2:623721/251‑1 (MQ=255)                                                                                                                                                                                                                                             ttGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGcc     <  2:627517/250‑1 (MQ=255)                                                                                                                                                                                                                                              tGAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGACTGCTGGCGAAGGCACTCGGTCTGCCGGTGAAACGGTTTATTGCTGCGACCCACGTGAACGATACCGTGCCACGGTTCCTGCCCGCCGGTCAGTGGTCC‑CCCAAAACGACTCAGGAGGCGTTATCCCACGCGATGGACGCGAGGCAGACGAACAACTGGCcgc   >  1:257674/1‑251 (MQ=255)                                                                                                                                                                                                                                               gAAGCAGTTGCGCAGCTGCCGCAGGA‑AGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACTCGGTCTGCCGGGGAAACGTTTTATTGCTGCGACCAACCGGAACGATAACGTGCCACGTTTCCTGCACGACGGGCCGTGGTCA‑CCCCAAGCGACTCAAGCGACGGTATCCAACGCGGTGGACGGGGGTCCACCGAACAACTTGCcgcc  >  1:502626/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                        ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGA‑GACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCA‑CCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG  >  NC_000913/4177‑4665

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