Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,930,364	G→A	W263* (TGG→TGA)	sdaB →	L‑serine deaminase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,930,364	0	G	A	92.0%	144.1 / ‑0.2	51	W263* (TGG→TGA)	sdaB	L‑serine deaminase II
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (20/26);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

AAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGGATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTG  >  NC_000913/2930123‑2930610                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         aaaCCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATg                                                                                                                                                                                                                                                 <  1:79315/251‑1 (MQ=255)   aCCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTAATGCGCGGCGGCATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGtt                                                                                                                                                                                                                                               >  1:158928/1‑251 (MQ=255)    ccGGGCTGTCACTCTCTGGCCTGATGATGACAATCGAGCTAGCGCTGCTCAGCAAAGAAGAGGTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGTCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCTCCACTGACCCGATGGCGGTTGTTTACTTAATCAACATGttt                                                                                                                                                                                                                                              <  2:98656/251‑1 (MQ=255)             cACTCTCTGGCCTGATGAGGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCa                                                                                                                                                                                                                                     <  1:236064/251‑1 (MQ=255)                          gatgatgaAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAagag                                                                                                                                                                                                                        <  1:370683/250‑1 (MQ=255)                          gatgatgaAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAgaga                                                                                                                                                                                                                       <  1:259719/251‑1 (MQ=255)                                gaAAAACGAGCTGGCGCTGCACAGCACAGAAGAGCTGGAACCGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGAGGTATTTCCACCGAAGGCGTGTTGCCTGGCACACTGCGCGTTCCACGCCGCGCTGCGGCCCTACGCCAGATGATTGTCAGCCCGGCTAAAACCCCCACTGACCCGCGGGCGGTTGTTGGATTAATCAACATGTTTGCACTAGCAGTG‑AAccgcgagacagcgg                                                                                                                                                                                                                 >  2:207752/1‑238 (MQ=255)                                    aaCGAGCTGGCGCTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCACGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGGATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTgg                                                                                                                                                                                                             <  1:51508/251‑1 (MQ=255)                                          cTGGCGCTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTcg                                                                                                                                                                                                       >  1:50533/1‑251 (MQ=255)                                                cTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCCCTGACCCGATGGCGGTTGTTGACTGGATACACATGTTTGCACAGGCAGGG‑AACGAAGAGAACGCTGCTGGCGGGcgcggggg                                                                                                                                                                                                 >  2:363199/1‑247 (MQ=255)                                                    gcaGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGGATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCTCTGACCCGATGGCGGGTGTTGACTGAATCAACATGTTTGCACGGGCAGTT‑AACTAAGAGAACTCTGCTGGCGGTCGCGTGGTGACt                                                                                                                                                                                             <  1:393719/250‑1 (MQ=255)                                                          aagaagAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGGGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGCTGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGACGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCg                                                                                                                                                                                       >  1:272754/1‑251 (MQ=255)                                                                   tGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTCACggg                                                                                                                                                                              >  1:333396/1‑250 (MQ=255)                                                                              cTGGCGAAAGTCTGGGAAGTAATGCGAGGCGGTATTGAGCGCGGTATGTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGATGGACAGCCAGGATAAAACCACCACTGACCCGCTGGCGGTTGGTGACTGCATCAACATTCTTGCACTGGCAGTGGAAC‑AAGAGAACGCTGCTGCCGGTCGCGTGGTGAATTCGCCGACTTACGgtgggcgagggat                                                                                                                                                                   >  2:310259/1‑241 (MQ=255)                                                                                 gCGAACGGCGGGGGAGTCAGGCGCGGCGGTATTGAGCGCGGGATTTCCACCGAAGGCGTGTGGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTAGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTAt                                                                                                                                                                <  2:17784/251‑1 (MQ=255)                                                                                    aaCGTCGGGGAAGTCATGCTCCGCGGTATGGATCGCGGTATTTCCACCGACGGCGTTTTCTCTGGCAAACTGCGCTTTCTACGCCGTGCTGCGGCTCTACGCCGTATGCTGGGCAGCCAGGATAAAACCACCACTGTCCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGGGGTGACTGCGCTGAATAACGGTGCGTGCGGGATTATcca                                                                                                                                                             <  2:348269/251‑2 (MQ=255)                                                                                       gTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTAAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGc                                                                                                                                                          <  1:316565/251‑1 (MQ=255)                                                                                         cTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACt                                                                                                                                                                                             <  1:254318/214‑1 (MQ=255)                                                                                               aGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTgg                                                                                                                                                  <  1:392858/251‑1 (MQ=255)                                                                                                  cATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCAGCGGCACTACGCCGGATGCTGGTCGGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGCCTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACGTCGCCGCATACAGGTGCGTGGGGGATTATccccgccgtatccgcgg                                                                                                                                               >  2:123188/1‑237 (MQ=255)                                                                                                     gcgcGGCGGTAGTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGGCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTTGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTTGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGtact                                                                                                                                            <  1:394711/251‑1 (MQ=255)                                                                                                                       cgGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGGATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGGGCGGGATTATCCCGGCAGTTCTGGCGTACTAAGACAAGTTTATCcgcg                                                                                                                          >  1:115855/1‑251 (MQ=255)                                                                                                                         gTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATAAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGa                                                                                                                         >  1:49945/1‑250 (MQ=255)                                                                                                                          tATTTCCACCGAGGGCGAGCTGCACAGCAAAATGCGACTTCCACGCCGCGCCGCGACACTCCGCCGGATGCTGGTAAGCCAGGATAACCCCACCACTGACCCGATGGCGGTTGATGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGGGGGGACTGCGCCGACTAACCGCGCGGGCGGGATTATCCCCGCCATTCAGGCCTACTACCAACAGTTTATcccggca                                                                                                                        >  2:452644/1‑245 (MQ=255)                                                                                                                             ttCCTCCGAGGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGTTGCTGGTCAGCTAGGGTAAAAGTACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGGGGGGGGTCGCGTGTGGGCTGCGCGGTCTAGCGGTGCGTGGGGGATTATCCCGGGAGTGCGGGCGtactgggacgattttattcttgaagtg                                                                                                                     <  2:184211/250‑24 (MQ=255)                                                                                                                              tCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCcgcggagggg                                                                                                                    >  1:9859/1‑244 (MQ=255)                                                                                                                               ccACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAAc                                                                                                                  <  1:104655/251‑1 (MQ=255)                                                                                                                                     aaGGCGCGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTTTCCGCGAAGGGAAAGCTAAc                                                                                                            >  1:118679/1‑251 (MQ=255)                                                                                                                                               gCCTGGCAAACGTCGCGTTCCACGCCGTGCTGCGGCACTCCGCCGGATGCGGGTCAGCCAGGATAAAACCACCACGGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTAGCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGcgc                                                                                                 <  2:287432/251‑3 (MQ=255)                                                                                                                                                      aaaCTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGAATGCTGGTCAGCCAGGATAAAACCACCACTGACCCCATGGCGGTTGTTGAATGAATCAACATGTTTGCACTGGCAGTG‑AACGACAAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCCACTAACGGGGCGGGGGGGATTATTCCGGCAGGTCTGGCGAAATACGACAAAGATATCCGCGACGTGAACGCTAACTC‑ACAGGcccggtaaac                                                                                          >  1:65864/1‑242 (MQ=255)                                                                                                                                                             gcgTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTa                                                                                   <  1:320085/251‑1 (MQ=255)                                                                                                                                                                   cACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAATCCACCACTGATCCGATGGCGGTTGTTGACTGAATCAACATGGTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCTGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAgc                                                                             <  1:412143/251‑1 (MQ=255)                                                                                                                                                                     cGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGGTCGGGATTATACCGGCAGTTCTGGCGGACTACGACAAGTTTATCCGCGAAAGGAACGCTAACTC‑ACTGGCTCGTTACCTTCTGGGAGCCAgcgc                                                                           >  1:13142/1‑251 (MQ=255)                                                                                                                                                                      gCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGcc                                                                          >  1:167540/1‑251 (MQ=255)                                                                                                                                                                        cGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCAt                                                                        <  1:24201/251‑1 (MQ=255)                                                                                                                                                                              gcGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGt                                                                    <  1:432752/249‑1 (MQ=255)                                                                                                                                                                               cGGCACTACGCCGGATGCTGGTAACCCAGGATAACACCACCACTGACCCGATGGAGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGCACGCTGCTGGCGGTCGCGTGGTGGCTGCGCCCACCAACGGGGCGGGCGGGATTATCCCGGCAGTTCTAGCGCACCAAGACAAGTTTATCCGCGGAGGGAACGCCAACTC‑ACCGGGTCGTTACatgccggcaaccaacccccttggttct                                                                 >  2:200806/1‑228 (MQ=255)                                                                                                                                                                                   aCTACGCCGGATGCTGGTCAGCCAGGCTAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGCAGAGAACGCTGCTGGCGGTCGCGTGGTGACTTCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGCTTATCCGCGAATGGACCGCTAACTC‑ACTGGCTCGTTAACTGCTTGTGGCCAGCGCCATTGGTTCTtttt                                                             >  2:118709/1‑251 (MQ=255)                                                                                                                                                                                    cTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGCAGAGAACGCTGCTGGCGGACGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGGTCTGGCGGACTACGACCAGTTTATCCGCGGAGTGGACGCTAACTCCAC‑GGGTCGTTGCCTGCTTGTAGCCCACGCCCTTTGTTCTCtttt                                                            >  2:159381/1‑250 (MQ=255)                                                                                                                                                                                                  tgtccgCCAGGATAAAACCTCCACTGTCCCGATGGCGGTTGTTGACTGGATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑TCTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATTAGATGATcgc                                                <  2:115855/244‑1 (MQ=255)                                                                                                                                                                                                  ggTCAGACAGGATAAAGCCACCACTGACCCGATGGCGGCTGTTGACTGAATCAACATGTTTGCACTGGCATTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCATTTCTGGCGTACTCCGACAAGTTTATCCGCGAAGTGAACTCTAACTC‑ACTTGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAAcgc                                                <  2:305080/249‑1 (MQ=255)                                                                                                                                                                                                        ccAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGGCTGCGCCGACTAACGGGGGGTGCGGGATTATCCCGGCAGTTCTGGCGGACTAAGACCAGTTTATCCGCGCAGGTAACGCTAACTC‑AATGGCCCGTTACCTGCTGGGCACCCGCGCCCTTGGTTTTCTTTTTAAGATGAAAGCGTCGGttt                                        >  2:455785/1‑251 (MQ=255)                                                                                                                                                                                                            gATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGCGTCGGTTTCTgg                                    >  1:199319/1‑251 (MQ=255)                                                                                                                                                                                                                aaatctccacTGTCCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCATTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑CCTTGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGCGTCGATTTCTGGTGca                                <  2:364227/245‑2 (MQ=255)                                                                                                                                                                                                                 aaCCACCACTGACCCGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTAGCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGCGTCGATTTCTGGTGCTg                               <  2:290361/251‑1 (MQ=255)                                                                                                                                                                                                                         cTGACCTGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGGAGCCAGCGCCATTAGTTCTCTTTTTAAGATGAACGCGTCGATTTCTGGGGCTGAAGTGGGt                       >  1:127632/1‑251 (MQ=255)                                                                                                                                                                                                                               cGATGGCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGGACTACGACAAGTTTATCCGCGAAGTGAACGCTAACGC‑ACTTGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGTGTCGATTTCTGGTGCTGAAGTGGGTTGCCAg                 <  1:452644/251‑1 (MQ=255)                                                                                                                                                                                                                                 aTGGCGGTTGTTGACTGAATTAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCTAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGCTTCGATTTCTGGTGCTGAAGTGGGTTGCCAggg               <  1:207752/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggCGGTTGTTGACTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGCGTCGATTTCTGGTGCTGGAGTGGGTTGCCAGGGTg             >  1:106082/1‑251 (MQ=255)                                                                                                                                                                                                                                    gCGGTTGTTGACTGAAGCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTTCGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGa            <  2:395116/251‑1 (MQ=255)                                                                                                                                                                                                                                              aCTGAATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTtggcgtg  <  1:233202/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         AAACCGGGCTGTCACTCTCTGGCCTGATGATGAAAAACGAGCTGGCGCTGCACAGCAAAGAAGAGCTGGAACAGCACCTGGCGAACGTCTGGGAAGTCATGCGCGGCGGTATTGAGCGCGGTATTTCCACCGAAGGCGTGTTGCCTGGCAAACTGCGCGTTCCACGCCGTGCTGCGGCACTACGCCGGATGCTGGTCAGCCAGGATAAAACCACCACTGACCCGATGGCGGTTGTTGACTGGATCAACATGTTTGCACTGGCAGTG‑AACGAAGAGAACGCTGCTGGCGGTCGCGTGGTGACTGCGCCGACTAACGGTGCGTGCGGGATTATCCCGGCAGTTCTGGCGTACTACGACAAGTTTATCCGCGAAGTGAACGCTAACTC‑ACTGGCTCGTTACCTGCTGGTAGCCAGCGCCATTGGTTCTCTTTATAAGATGAACGCGTCGATTTCTGGTGCTGAAGTGGGTTGCCAGGGTGAAGTTGGCGTG  >  NC_000913/2930123‑2930610

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