Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,142	C→T	C114C (TGC→TGT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,142	0	C	T	98.3%	175.9 / ‑4.3	60	C114C (TGC→TGT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (34/25);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

TGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCG  >  NC_000913/2902‑3387                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       tGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGc                                                                                                                                                                                                                                              >  1:299963/1‑251 (MQ=255)         cGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGac                                                                                                                                                                                                                                      >  1:429023/1‑251 (MQ=255)          ggTTGCGGCGGCAGAGAAATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGCACCCCGGGAAAATATCGTTTATCAGTGCTGGGCGCGTTTTTGCCAGGCGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGTTGGCGGTGAATGCACACTGTGGCAAGCCGCTTAATGac                                                                                                                                                                                                                                      >  2:23863/1‑250 (MQ=255)              gAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCg                                                                                                                                                                                                                                 >  2:424535/1‑251 (MQ=255)                  cggcAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTACTGACACTCGTTTg                                                                                                                                                                                                                             >  2:136292/1‑251 (MQ=255)                     ccgagaCATTCAGGCTCAACAACCTCGGTCGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATTTCGTTTATCATTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACTCTGTGGCAAGCCGCTTAATGACACTCGTTTGCTg                                                                                                                                                                                                                          <  2:475592/249‑1 (MQ=255)                           cATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGATCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTg                                                                                                                                                                                                                    >  2:259270/1‑251 (MQ=255)                            aTTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGTGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACTCTCGTTTGCTGGCTTTGa                                                                                                                                                                                                                   <  1:161723/251‑1 (MQ=255)                                         aacCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCAAGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGGGGTCGCGGCGCTGATGGCGATGAATGACCACTGTGGCAAGCCGCTTAATGCCACTCGTATGCTGGCTTTGATGGGCGAGCTgg                                                                                                                                                                                                       >  1:128926/1‑250 (MQ=255)                                           ccccggacgtctggCCGTTAAGCTGCCTGCAGAACCCCGGGAAACTATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCGTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAg                                                                                                                                                                                                    <  2:158633/238‑1 (MQ=255)                                              cGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAggac                                                                                                                                                                                                 <  1:302608/251‑3 (MQ=255)                                                gACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAggacg                                                                                                                                                                                                >  1:348065/1‑247 (MQ=255)                                                      ttGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTGCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAAtct                                                                                                                                                                                           <  1:215774/249‑1 (MQ=255)                                                        gCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGCAGGACGAATCTCCgg                                                                                                                                                                                       >  1:147606/1‑251 (MQ=255)                                                            aTAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGgcagc                                                                                                                                                                                   >  1:411007/1‑251 (MQ=255)                                                              atgCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGTGATGATTCTGGAAAAGAATATGCCAATCGGTTCGGTCTTAGTCTCCAGCGCCCGTTCTGTGGTCGCGGCGCTGATGGCGTTGAATGAACGCTGTGGTAAGCCGCTTAATGACACTCGTTTGCTGGCTTGGGTTGGCGAGCTGGAAGGACGAATCTCCGgcagaat                                                                                                                                                                                 <  1:23326/249‑4 (MQ=255)                                                                gCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGcattc                                                                                                                                                                               >  2:406777/1‑251 (MQ=255)                                                                     cGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTAc                                                                                                                                                                          <  2:272209/251‑1 (MQ=255)                                                                          gAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGACAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGCTGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACa                                                                                                                                                                      >  2:46101/1‑250 (MQ=255)                                                                                         atCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGttt                                                                                                                                                      <  1:424535/251‑1 (MQ=255)                                                                                           cGTTTATCAGTGCTGGGGGCGTGTGTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACGCGGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGGTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGTCACCGTGtttt                                                                                                                                                     <  2:3633/250‑1 (MQ=255)                                                                                                 tCAGTGCTGGGAGCGTTTTTGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTggtg                                                                                                                                              >  1:195949/1‑251 (MQ=255)                                                                                                                tttttGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGCCACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGAACCGTGTTTTCTTGGTGGTATTCAGTtgatg                                                                                                                                >  2:358672/1‑250 (MQ=255)                                                                                                                   ttgtccggagcttgggaaGCAAACTCCCGTGGCGATGACTCTTGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGACTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACTCTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTTCGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATc                                                                                                                             <  2:119142/234‑1 (MQ=255)                                                                                                                   ttGCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCAGTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATc                                                                                                                             <  1:246690/250‑1 (MQ=255)                                                                                                                   ttGCCAGGAGCTGTGCAAGCAAATTCCAGTGGCGATGACTCGGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATc                                                                                                                             <  2:429023/250‑1 (MQ=255)                                                                                                                     gCCAGGAGCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCgaa                                                                                                                          <  1:217992/251‑1 (MQ=255)                                                                                                                        aggagcttggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACCCTGTGGCAAGCCGCTTACTGACAATCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACGACCTGGCCACGTGTTTTATTGGTGGTATGCAGTTGATGATCgaaga                                                                                                                        >  2:192633/12‑250 (MQ=255)                                                                                                                        aggagcttggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCCTTCATTACGACAACGTGGCCCCGTGTTTTCTTGGTGGTATGCAGTTGATGATCgaaga                                                                                                                        >  2:243706/12‑250 (MQ=255)                                                                                                                                ggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACTCTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAgg                                                                                                                                                                                                   <  1:394805/164‑1 (MQ=255)                                                                                                                                ggcaaGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTTATGGGCGAGCTGGAAgg                                                                                                                                                                                                   >  2:394805/4‑167 (MQ=255)                                                                                                                                ggcaaGCAAATGCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCTGTGGTCGAGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAATGACGAATCTCCGGCAGCATTCATTACGACAACGCGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAACACGCcat                                                                                                               <  2:154563/248‑1 (MQ=255)                                                                                                                                            ccAGTGGCGATGACTCTGGAAAAGAATATGCCAAGCGGTTCGGGCTTAGGCTCCAGCGCCTGTGCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATcagccagc                                                                                                    <  2:121006/250‑1 (MQ=255)                                                                                                                                               gTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAgg                                                                                                >  1:81131/1‑250 (MQ=255)                                                                                                                                                   cgttGACTCTGGAAAAGAATATGCCCATCGGTTGGGGCATAGGCACCAGCACCTGTTCGACGGTCGCGGAGCTGACGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGGGCTGGATGGACGATTCTCCGGCAGCATTCATTACGACAACGAGGCACCGTGATTTCTTGGTGGTATGCAATCGCTGATCGAAGAAAACCACATCATCAGCCAGCAAGTggcg                                                                                            >  2:174705/4‑248 (MQ=255)                                                                                                                                                       gactcTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCATTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAggg                                                                                         <  1:143692/246‑1 (MQ=255)                                                                                                                                                            tGGAAAAGAATATGGAAAGCGGTTAGGGATTAGGCTAAAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGACTGAACACTGTGGCAAGCCGCTTAATGAAACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTtgat                                                                                   >  1:358998/1‑251 (MQ=255)                                                                                                                                                             ggAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGGGCCAGGGTTtgatg                                                                                  >  2:338302/1‑251 (MQ=255)                                                                                                                                                                aaaGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGt                                                                               >  1:310392/1‑251 (MQ=255)                                                                                                                                                                aaaGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGt                                                                               <  1:414371/251‑1 (MQ=255)                                                                                                                                                                        tgccactCGGCTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGTTGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGCGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAGGAAAACGACATCATCAGCCAGCCGGTGCCGGGGTTTGATGAGTGGCTGTgg                                                                       >  2:213645/7‑251 (MQ=255)                                                                                                                                                                                   tCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGt                                                            >  2:216472/1‑250 (MQ=255)                                                                                                                                                                                     gggCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGttt                                                                                                                                                      >  2:359051/1‑159 (MQ=255)                                                                                                                                                                                     gggCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGttt                                                                                                                                                      <  1:359051/159‑1 (MQ=255)                                                                                                                                                                                      ggCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATcc                                                        <  1:230096/251‑1 (MQ=255)                                                                                                                                                                                       gCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCg                                                       >  1:347416/1‑251 (MQ=255)                                                                                                                                                                                       gCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCg                                                       >  1:228234/1‑251 (MQ=255)                                                                                                                                                                                        cTTAGGCTCCGGCGCATGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACGGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCGGGAAGGACGAATCTCCGGCAGCGTTCATTACGACAACGTGGCACCGTGTTTTCTTGGAGGTATGCAGTTGATGATCGAAGAAACCGACATCATCAGCAAGCCAGTGCCAGGGTTTGATGAGTGGCTGTGGGTG‑TGAGCGTATCCgg                                                      >  2:129041/1‑251 (MQ=255)                                                                                                                                                                                                  agcgCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGt                                            >  1:1796/4‑251 (MQ=255)                                                                                                                                                                                                  agcgCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGATGTGGGTGCTG‑GCGTATCCGGGGGTTCAAGt                                            >  2:91435/4‑251 (MQ=255)                                                                                                                                                                                                     gCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGtctc                                         >  1:118650/1‑251 (MQ=255)                                                                                                                                                                                                      ccTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAATGTCTCg                                        >  1:266239/1‑251 (MQ=255)                                                                                                                                                                                                         gtgcGGTGGTCGCGGAGCTGAGGGCGTTGAATGAACACAGTGGCAAGCCGCGTAATGACTCTTGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAAGCTCCGGCAGGATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGTCTGTGGGTGCTG‑GCGTTTCCGGGGATTAAAGTCTCGAc                                      <  2:358998/247‑1 (MQ=255)                                                                                                                                                                                                                           tGATGGCGATGAATGAACACTGTGGCACGCCTCTTACTGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTAATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAACCGACATCATCAGCCAGCACGTGCCAGGGATTGATGAGTGGCTGTAGGTGCTG‑GCGAATCCGGGGATTAAAGTCTCGACGGCAGACACCAGGGCTAtt                   >  2:78141/1‑251 (MQ=255)                                                                                                                                                                                                                            gATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTAtt                   <  1:255399/250‑1 (MQ=255)                                                                                                                                                                                                                             aTGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGGTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTAtttt                 <  1:23863/251‑1 (MQ=255)                                                                                                                                                                                                                              tGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTCtttt                 >  1:483215/1‑250 (MQ=255)                                                                                                                                                                                                                                 cGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCg             <  1:156832/251‑1 (MQ=255)                                                                                                                                                                                                                                     gaatgaaCACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGgcg          <  1:456506/250‑1 (MQ=255)                                                                                                                                                                                                                                            cacTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATcg  <  1:406777/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       TGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGAACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCAGGGTTTGATGAGTGGCTGTGGGTGCTG‑GCGTATCCGGGGATTAAAGTCTCGACGGCAGAAGCCAGGGCTATTTTACCGGCGCAGTATCG  >  NC_000913/2902‑3387

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