Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,025	A→G	E75E (GAA→GAG)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,025	0	A	G	98.0%	140.5 / NA	50	E75E (GAA→GAG)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (28/21);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑ACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAA  >  NC_000913/2796‑3268                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    cTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGGTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGTTGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCCACACCATCGGACGCGTTGACGATAAGCTGACGTCCGAACCCCGGGAAAATATCGTTTATCCGTGCTTGGAGCGGTTTTGCCAGGA‑GCTGGGAACGCAAATTCCAGTg                                                                                                                                                                                                                                >  2:326903/1‑251 (MQ=255)    aCATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCg                                                                                                                                                                                                                             <  1:138113/251‑1 (MQ=255)         gTTAAAGTTTATGCCCCGGCTTCCAGGGCCAATATGAGCGTCGGGGTTGATGTGCTCTAGGCGGCGGTGACACCTGTTGATGGGGCATTGCTCGGAGTTGTAGGCACGGTTGAGGCGGCAGGTACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATCAGCTGCCGTCAGAAACACGGGAAAATATCGTTTATCCGAGCCGGGAGCGGTTTTGCCAGGA‑GCTTGGCCAGCAAATTACTGTGGCggaggc                                                                                                                                                                                                                        >  2:414371/1‑246 (MQ=255)                      ccccGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAat                                                                                                                                                                                                           <  1:85710/251‑1 (MQ=255)                         cGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGGAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAAACTCGGACGCTTTTCCGGTAAGCTGCCGTCAGAACCCAGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCCGTGGCGCTGACTCTGGGAAAGAAAATg                                                                                                                                                                                                        >  2:230096/1‑251 (MQ=255)                               ccAGGGCCAATATGTGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCGCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATc                                                                                                                                                                                                  <  2:486221/251‑1 (MQ=255)                                cAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTTCCGATAAGCTGCCGTCAGAACCCCGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCCAATTCCAGTGGCGATGACTCTGGGAAAGAATATGCCAATCg                                                                                                                                                                                                 >  2:195378/1‑251 (MQ=255)                                 aGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCa                                                                                                                                                                                                                                   <  2:234159/216‑1 (MQ=255)                                 aGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCa                                                                                                                                                                                                                                   >  1:234159/1‑216 (MQ=255)                                      cAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCgg                                                                                                                                                                                           >  1:121006/1‑251 (MQ=255)                                       aaTATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGGCATTCAGTCTCAACAACCACGGACGCTTTGACGATAAGCTGACGTCAGAACCACGGGAAAATATCGGTTATCAGAGCTGGGAGGGAGTTTGCCCGGAGACTGG‑CAAGCACATTCCAGTGGGGGTGACTCTGGGAAAGAattgccaaatcgattcggg                                                                                                                                                                                          >  2:386380/1‑234 (MQ=255)                                       aaTATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGCGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCggg                                                                                                                                                                                          <  1:176078/251‑1 (MQ=255)                                          atGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCtt                                                                                                                                                                                       <  2:136238/251‑1 (MQ=255)                                             aGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGGGCGGTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAgg                                                                                                                                                                                    >  1:366256/1‑251 (MQ=255)                                                     gttggATGTGCTCGGGGCGGCGGTGACACCTGTTGCTGGTGCCTTGCTCGGAGATGGAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTcaagcg                                                                                                                                                                            >  2:217992/5‑246 (MQ=255)                                                     gattgctgtgCTCGGGGCGGCTGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGGCATTAAGGCGCAACAACATCGGACGCTTTGACGGTAAGCAGCAGGCAGAACCACGGGAGAAAATCTGTTATCAGTGCTGGGAGCGTTTTTGACAGGA‑GCTTGGCAAGAAACTACCAGTGGCGATGCCTCTGGAAAAGTACATTCCCATCGGCTTGGGCTTAGGGTCCaacg                                                                                                                                                                            >  2:362061/7‑248 (MQ=255)                                                            gtgCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTc                                                                                                                                                                     <  1:100320/251‑1 (MQ=255)                                                                       ggcggTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCgcgg                                                                                                                                                          >  1:144697/1‑251 (MQ=255)                                                                       ggcggTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCgcgg                                                                                                                                                          <  1:195378/251‑1 (MQ=255)                                                                           gTGACACCTGTTGATGGAGCATTGCTCGGAGATATAGTCACCGTTGAGGCGGCCGCGACTTTCAGCCCCAACAACCTCGGACGCTTTGCCGATACGCCGCCGTCAGAACCACGGGAAAATATAGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGAAAGCCAATTCCCGTGGCGATGACTCTGGCAAAGAATTTGCCAATAGGCTAGGGGTTAGGCTCCAACGCCCGTTCGATGGTCggggcgat                                                                                                                                                      >  2:207515/1‑249 (MQ=255)                                                                            tGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGAc                                                                                                                                                                                                                        <  1:479978/184‑1 (MQ=255)                                                                               cacCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGACCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCATTGGCGATGCCTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATg                                                                                                                                                  >  1:269674/1‑251 (MQ=255)                                                                                              cATTGCTCGGAGATGAAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTGTGCCGCTAAGCTGCCGTACGCACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGGGGCGACGACTTTGGAAAAGAATATGCCAATCGGATCGGGCGTAGGCACCAGCGCCTGTTTGGTGGTCGCGGCGCAGGTTGCGATGAATGTacac                                                                                                                                   >  2:313586/1‑251 (MQ=255)                                                                                                      ggAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTCTCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGGGGTCGCGGCGCTGTAGGCGATGAATGAacacggtggca                                                                                                                            >  2:85660/1‑243 (MQ=255)                                                                                                           tGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGc                                                                                                                      >  1:299963/1‑251 (MQ=255)                                                                                                                   cGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGac                                                                                                              >  1:429023/1‑251 (MQ=255)                                                                                                                    ggTTGCGGCGGCAGAGAAATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGCACCCCGGGAAAATATCGTTTATCAGTGCTGGGCGCGTTTTTGCCAGGC‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGTTGGCGGTGAATGCACACTGTGGCAAGCCGCTTAATGac                                                                                                              >  2:23863/1‑250 (MQ=255)                                                                                                                        gAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCg                                                                                                         >  2:424535/1‑251 (MQ=255)                                                                                                                            cggcAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTACTGACACTCGTTTg                                                                                                     >  2:136292/1‑251 (MQ=255)                                                                                                                               ccgagaCATTCAGGCTCAACAACCTCGGTCGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATTTCGTTTATCATTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACTCTGTGGCAAGCCGCTTAATGACACTCGTTTGCTg                                                                                                  <  2:475592/249‑1 (MQ=255)                                                                                                                                     cATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGATCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTg                                                                                            >  2:259270/1‑251 (MQ=255)                                                                                                                                      aTTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGTGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACTCTCGTTTGCTGGCTTTGa                                                                                           <  1:161723/251‑1 (MQ=255)                                                                                                                                                   aacCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCAAGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGGGGTCGCGGCGCTGATGGCGATGAATGACCACTGTGGCAAGCCGCTTAATGCCACTCGTATGCTGGCTTTGATGGGCGAGCTgg                                                                               >  1:128926/1‑250 (MQ=255)                                                                                                                                                     ccccggacgtctggCCGTTAAGCTGCCTGCAGAACCCCGGGAAACTATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCGTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAg                                                                            <  2:158633/238‑1 (MQ=255)                                                                                                                                                        cGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAggac                                                                         <  1:302608/251‑3 (MQ=255)                                                                                                                                                          gACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAggacg                                                                        >  1:348065/1‑247 (MQ=255)                                                                                                                                                                ttGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTGCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAAtct                                                                   <  1:215774/249‑1 (MQ=255)                                                                                                                                                                  gCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGCAGGACGAATCTCCgg                                                               >  1:147606/1‑251 (MQ=255)                                                                                                                                                                      aTAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGgcagc                                                           >  1:411007/1‑251 (MQ=255)                                                                                                                                                                        atgCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGTGATGATTCTGGAAAAGAATATGCCAATCGGTTCGGTCTTAGTCTCCAGCGCCCGTTCTGTGGTCGCGGCGCTGATGGCGTTGAATGAACGCTGTGGTAAGCCGCTTAATGACACTCGTTTGCTGGCTTGGGTTGGCGAGCTGGAAGGACGAATCTCCGgcagaat                                                         <  1:23326/249‑4 (MQ=255)                                                                                                                                                                          gCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGcattc                                                       >  2:406777/1‑251 (MQ=255)                                                                                                                                                                               cGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTAc                                                  <  2:272209/251‑1 (MQ=255)                                                                                                                                                                                    gAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGACAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGCTGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACa                                              >  2:46101/1‑250 (MQ=255)                                                                                                                                                                                                   atCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGttt                              <  1:424535/251‑1 (MQ=255)                                                                                                                                                                                                     cGTTTATCAGTGCTGGGGGCGTGTGTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACGCGGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGGTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGTCACCGTGtttt                             <  2:3633/250‑1 (MQ=255)                                                                                                                                                                                                           tCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTggtg                      >  1:195949/1‑251 (MQ=255)                                                                                                                                                                                                                          tttttGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGCCACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGAACCGTGTTTTCTTGGTGGTATTCAGTtgatg        >  2:358672/1‑250 (MQ=255)                                                                                                                                                                                                                             ttGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCAGTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATc     <  1:246690/250‑1 (MQ=255)                                                                                                                                                                                                                             ttGCCAGGA‑GCTGTGCAAGCAAATTCCAGTGGCGATGACTCGGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATc     <  2:429023/250‑1 (MQ=255)                                                                                                                                                                                                                               gCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCgaa  <  1:217992/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    CTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑ACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAA  >  NC_000913/2796‑3268

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