Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,055	C→T	T85T (ACC→ACT)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,055	0	C	T	98.1%	148.8 / NA	53	T85T (ACC→ACT)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/21);  minor base G (0/1);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑ACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCA  >  NC_000913/2817‑3298                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    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)                 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:195378/251‑1 (MQ=255)                                                  ggcggTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCgcgg                                                                                                                                                                                        >  1:144697/1‑251 (MQ=255)                                                      gTGACACCTGTTGATGGAGCATTGCTCGGAGATATAGTCACCGTTGAGGCGGCCGCGACTTTCAGCCCCAACAACCTCGGACGCTTTGCCGATACGCCGCCGTCAGAACCACGGGAAAATATAGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑GCTTGGAAAGCCAATTCCCGTGGCGATGACTCTGGCAAAGAATTTGCCAATAGGCTAGGGGTTAGGCTCCAACGCCCGTTCGATGGTCggggcgat                                                                                                                                                                                    >  2:207515/1‑249 (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)                                                                                                                                                                                                         ttgtccggagcttgggaaGCAAACTCCCGTGGCGATGACTCTTGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGACTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACTCTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTTCGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATc                                   <  2:119142/234‑1 (MQ=255)                                                                                                                                                                                                          gCCAGGA‑GCTTGGCAAGCAAATTCCAGTGGCGATGACTCTGGAAAAGAATATGCCAATCGGTTCGGGCTTAGGCTCCAGCGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGTGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGACGAATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTTGGTGGTATGCAGTTGATGATCgaa                                <  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)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    CCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTATCAGTGCTGGGAGCGTTTTTGCCAGGA‑ACTGGGTAAGCAAATTCCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGAAGGCCGTATCTCCGGCAGCATTCATTACGACAACGTGGCACCGTGTTTTCTCGGTGGTATGCAGTTGATGATCGAAGAAAACGACATCATCAGCCAGCAAGTGCCA  >  NC_000913/2817‑3298

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