Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,625,248	C→T	D118D (GAC→GAT)	serB →	phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,625,248	0	C	T	96.0%	145.6 / ‑1.9	50	D118D (GAC→GAT)	serB	phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (28/20);  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 = 5.48e-01

GGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTG  >  NC_000913/4625009‑4625487                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                ggctggctGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCAACTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCTCCAt                                                                                                                                                                                                                                      >  1:258279/1‑251 (MQ=255)       ctgctgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATcc                                                                                                                                                                                                                                    >  2:324848/1‑247 (MQ=255)                ggTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCGAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCCACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCCCGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCCTCCAGATTGAATGTAt                                                                                                                                                                                                                       >  2:151433/1‑251 (MQ=255)                   cGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATtg                                                                                                                                                                                                                     <  1:324848/250‑1 (MQ=255)                          tGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAAtt                                                                                                                                                                                                             >  2:154791/1‑251 (MQ=255)                            gATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCGCATTAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGGTTGTTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGc                                                                                                                                                                                                           <  1:461333/251‑1 (MQ=255)                                   agcGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTg                                                                                                                                                                                                    >  1:157111/3‑251 (MQ=255)                                   agcGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGCGATTGTTGCCGCCTGGCGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGCTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGAGCAACTGGATGCCGCACCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGCGATGGACATGGATTCAACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTg                                                                                                                                                                                                    >  1:347024/3‑251 (MQ=255)                                       ccTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCCCCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCg                                                                                                                                                                                                <  1:423581/250‑1 (MQ=255)                                               aaTACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGc                                                                                                                                                                                        >  1:113615/1‑251 (MQ=255)                                                     agagCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAtg                                                                                                                                                                                  <  2:162532/251‑1 (MQ=255)                                                         cAAACTGGGTGCGGCGATGGTGATTGTGGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAtggtgg                                                                                                                                                                              <  1:421558/251‑1 (MQ=255)                                                          aaaCTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGc                                                                                                                                                                             <  1:114399/251‑1 (MQ=255)                                                           aaCTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCg                                                                                                                                                                            <  1:175268/251‑1 (MQ=255)                                                                 gTGCGTCGATTGTGATTGTTGCCGCCGGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTa                                                                                                                                                                      <  1:500240/251‑1 (MQ=255)                                                                  tGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCCCCGCACGGGCTCCGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAgga                                                                                                                                                                      >  2:69579/1‑248 (MQ=255)                                                                   gcggcgATGGTGATTGTTGCCGCCTGGGGCGTGGAAGATGATCAGGTGATTCGGCTGGCAGGTGCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTaa                                                                                                                                                                     <  2:234971/250‑1 (MQ=255)                                                                        gATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAAc                                                                                                                                                               >  1:33711/1‑251 (MQ=255)                                                                             ggATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGct                                                                                                                                                          <  1:380313/250‑2 (MQ=255)                                                                              gATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGCACGGGCgg                                                                                                                                                         >  2:161219/1‑250 (MQ=255)                                                                                          cTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAAc                                                                                                                                             >  1:215092/1‑251 (MQ=255)                                                                                                           aTCAGGGGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTAcagcacgc                                                                                                                            >  2:362486/1‑244 (MQ=255)                                                                                                            tCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGGGATGGTGGCGGAAGTAACCGAACGGGGGATGCGCGGCGAACTCGCTTTTACCGCCAGcc                                                                                                                           >  2:185881/1‑251 (MQ=255)                                                                                                              aGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTg                                                                                                                         <  1:259566/251‑1 (MQ=255)                                                                                                                gTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTgcg                                                                                                                       >  1:108666/1‑251 (MQ=255)                                                                                                                 tGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCGTGATGCGCATCTGGCTGTCGCCCCGCTGGGGCAAATCCCGCACCTGCGCACGCAGGGTTTGATAGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGAGCGAGCTGGTGGCGGAAGTAACCGAACGTGCGATGCGCGGCGAACTAGATGTTACCGCCAGCCTgcgc                                                                                                                      >  1:452102/1‑251 (MQ=255)                                                                                                                      cGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCGGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGcc                                                                                                                  <  1:422864/250‑1 (MQ=255)                                                                                                                        tCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGACAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACCCGGTTTTACCGCCAGCCTgccccgccgcg                                                                                                               >  2:449862/1‑242 (MQ=255)                                                                                                                         cTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCgcgt                                                                                                              >  1:60887/1‑248 (MQ=255)                                                                                                                             ccggTTCACTCACCGCACGGGCTACGCGCCTGGCCCGAGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAAGCCCGCCCCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCg                                                                                                          <  1:208754/249‑1 (MQ=255)                                                                                                                                     ctcacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGACTGGGCAGCCGCGAGGCGACGcggac                                                                                                  >  2:272965/1‑247 (MQ=255)                                                                                                                                            cACGGGCTACGCGCCTGGCCCATGACGAGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGCGCTGGCCATGGCTTCCACCGCCATCCAGATTGCATGAATTGATGAAATTGCCAAACTGGTCGGGACGGGAGCGATGGTGGCGGGAGTAACCGAACGGGCGAAGCGCGGCGAACTCGATTTTCCCCCCAACCTGCGCAGccacgggacgccgctgaaaggggat                                                                                           >  2:279469/1‑228 (MQ=255)                                                                                                                                               gggCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGcc                                                                                                                                                                                                          <  2:140238/137‑1 (MQ=255)                                                                                                                                               gggCTACACGCCTGGACCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGcc                                                                                                                                                                                                          >  1:140238/1‑137 (MQ=255)                                                                                                                                                           gggCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAACCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCCCCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTg                                                                            <  1:161219/250‑1 (MQ=255)                                                                                                                                                            ggCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCCATATTCTGc                                                                           >  2:481666/1‑251 (MQ=255)                                                                                                                                                             gcccatgaggcgcacctggagtgcGCCCCTCGGGGGAAAACCCCGCACCTGCGCACGCGGGGTTGGCTGGTGATGGACATGGTTTCCACCGCCATCCAGATTGAAGGTATTGATGAAATTGCCAACCTGGCCGGAACGGGCGAGAGGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCa                                                                          <  2:148632/228‑1 (MQ=255)                                                                                                                                                                catgaagcgcaatgggaggTCGCCCCGCTGGGGAAAATCCCGCACCGGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCCCCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACGCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACa                                                                       <  2:215092/233‑1 (MQ=255)                                                                                                                                                                        gcaacTGGAGGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATGGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTg                                                               <  1:457862/247‑1 (MQ=255)                                                                                                                                                                                       cccGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGCTGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTg                                                 >  1:284756/1‑250 (MQ=255)                                                                                                                                                                                               ggAAAATACCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTCTTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGc                                        >  1:329291/1‑251 (MQ=255)                                                                                                                                                                                                 aaaaTCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGGGCGTGAAAATCTGCCGCTGATGCCAGGCtt                                      >  2:132375/1‑251 (MQ=255)                                                                                                                                                                                                     ccccGCACCAGCGCACGGGGGGTTTGCTGTTGATGGACTTGGATTCCACCGCCCTCCAGATTGAAGGTATTGATGAAATGACCAAACTGGCCGGAACGGGCGAGATGGTGGGGGAAGTAACCGAACGGGCGATGCGCGGCGAACCCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGTTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTCACg                                  <  2:93015/250‑1 (MQ=255)                                                                                                                                                                                                         gCACCTGCGCACGCCGGGTTTGCTGGGGATGGACATGGATTCCACCGCCATCCAGATTGCATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGAGGCGGAAGTAACCGAACGGGCGATGCGCGCCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGACGGGCGCTGACGCCCATATTCTGCAACAGGGGCGTGAAAATCTTCCGCGGATGCCCGGCTTAACGCAAc                              >  2:69211/1‑251 (MQ=255)                                                                                                                                                                                                          ccccTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACt                             <  1:490294/249‑1 (MQ=255)                                                                                                                                                                                                           aCCTGCGAACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGCTTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAGGTAACCGGACGGGCGATGCGAGGCGAACTCGATTTTACCGCGGGCCTGCGCAGCCGCGTGGCGACGCTGAAAGACGCTGACGCCAATATTCTGCCAAAGGTACCTGAAAATCTGCCGCAGATGcccgcagtaaggcaacgg                            >  2:15678/1‑235 (MQ=255)                                                                                                                                                                                                                ggcACGCCGGGTTGGCTGGTGATGGAGATGGATTCCACCGCCATCCAGATTTAATGTATTGATGAAATGGCAAACCGGGCGGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCt                       <  2:200911/250‑1 (MQ=255)                                                                                                                                                                                                                  cACGACGGGTGTGCTGGTGATGGACCTGGATTCCACCGCCATCCAGATTGAATGAATTGATGAAATTGCCAAACTGGCCGGCACGGGCGAGGTTGGGGCGGAAGTAACCGAACGAGCGATGCGCAGCGAACTCGATTTTACCGCCCGACTGCGCAGCCGCGTGGCGACGCTGAAAGGAGCTGACGCCAATATTCTGCAACAGGTGCGTGACAATCTGCAGCTGATGCCAGGcatacctcgactggtacac                      >  2:75997/1‑232 (MQ=255)                                                                                                                                                                                                                              gCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgc         >  2:412162/1‑251 (MQ=255)                                                                                                                                                                                                                                     aTGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACgctgggctg  >  1:380973/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                GGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATTCTGCAACAGGTGCGTGAAAATCTGCCGCTGATGCCAGGCTTAACGCAACTGGTGCTCAAGCTGGAAACGCTGGGCTG  >  NC_000913/4625009‑4625487

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