Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,851,841	C→T	Q328* (CAA→TAA)	ansA →	L‑asparaginase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,851,841	0	C	T	98.0%	163.0 / NA	50	Q328* (CAA→TAA)	ansA	L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/29);  minor base G (1/0);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTG  >  NC_000913/1851601‑1852085                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     ccACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTgcg                                                                                                                                                                                                                                            >  2:168758/1‑251 (MQ=255) ccACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTgcg                                                                                                                                                                                                                                            >  1:62524/1‑251 (MQ=255)        acaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGaa                                                                                                                                                                                                                                     <  1:278460/251‑1 (MQ=255)          aaaGCCTTCCTGCAGGAATTAAAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTACTTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCAATACCTGCTTATCCAGGAACTGGAAACTGAAAACATTCGACAGGCCTTGAGATAAAACCTGCGAGGCGAACt                                                                                                                                                                                                                                   >  2:410394/1‑251 (MQ=255)           aaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTg                                                                                                                                                                                                                                  <  1:362835/251‑1 (MQ=255)             gcctgcctGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTTGTGGTCAACCTGACACAATGTATGGCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATTACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGAc                                                                                                                                                                                                                                <  1:541445/251‑1 (MQ=255)                      cAGGAATTACAAGAAGCCAGCGATCGCGGGTTTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGGGGTTATGCCACCGGTAGTGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTTCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAGCCTGCGCGGCGAACTGACTCtttgctc                                                                                                                                                                                                                       <  2:329718‑M1/251‑8 (MQ=255)                          aaTTACAAGAAGCGATCGATCGCGGGATGGTGGTGGTCACCCTGACAGAAGGTATGTCCGGTAAAGTGAACATGGGTGGTTAGGCCGCCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGATCTGGATACTGAAGCCATCCGCAAGGCCATGAGCTAAATCCTGCGCGGCGACCTGACGCtttgctcacat                                                                                                                                                                                                                   <  2:486577‑M1/251‑12 (MQ=255)                             tACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatct                                                                                                                                                                                                                 <  2:139940‑M1/250‑14 (MQ=255)                                   aagCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcacttt                                                                                                                                                                                                           >  1:435558‑M1/1‑231 (MQ=255)                                        aGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgt                                                                                                                                                                                                     <  1:110602‑M1/251‑26 (MQ=255)                                                cgGTATGGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacat                                                                                                                                                                                             <  2:57520‑M1/251‑34 (MQ=255)                                                 ggTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTTCCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacatt                                                                                                                                                                                            <  1:557227‑M1/251‑35 (MQ=255)                                                  gTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTACCGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGGAGAAGCCACGCTAACCCACCTGCATTCCCTTCTTAGCCAGGAACTAGATACTGAAACCATTCGCAAGGCCATGCGCGAAAACCTGCGCGGCGAACTGACGCtttnctcccacctccctttaatcgtgctcccatt                                                                                                                                                                                            >  2:275061‑M1/1‑216 (MQ=255)                                                             gtCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgatt                                                                                                                                                                                <  1:97084‑M1/251‑47 (MQ=255)                                                             gtCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgatt                                                                                                                                                                                <  1:221602‑M1/251‑47 (MQ=255)                                                                               atgtCCGGTAAAGTGAACCTGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATCTGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATACGCAAGGACATGCGCTAAAACCTGCGCGGCGAACTGACGCtttccccacatctcactttaatagtgctcacattacgtgaatgattctaacaaaaccttaacac                                                                                                                                                              >  2:152281‑M1/1‑187 (MQ=255)                                                                                tgtCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtaactgattctaacaaacccttaacacc                                                                                                                                                             >  2:127086‑M1/1‑186 (MQ=255)                                                                                      gTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactgg                                                                                                                                                       <  1:33982‑M1/251‑72 (MQ=255)                                                                                             gAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaat                                                                                                                                                 <  1:191401‑M1/250‑78 (MQ=255)                                                                                                   gggTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtc                                                                                                                                           <  1:35952‑M1/250‑84 (MQ=255)                                                                                                    ggtggtTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGCTACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaacaactggcaaaattttgtcct                                                                                                                                         >  1:331278‑M1/1‑166 (MQ=255)                                                                                                    gggggttatGCCACCGGTAACGCCCCCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCCTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcct                                                                                                                                         <  2:517570‑M1/248‑86 (MQ=255)                                                                                                         ttATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACATGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaact                                                                                                                                    >  1:390472‑M1/1‑161 (MQ=255)                                                                                                          tATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgttactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaactt                                                                                                                                   >  1:153071‑M1/1‑160 (MQ=255)                                                                                                             gCCACCTGTAACGCCCTCGCCCATGCCGGCGTAATTTGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCGGCGCGGCGAACTGACGCtttgctcacatctccctttaatcgtgctcacattacgtgactgattctaacaaaacattatcaccaactggcaaaattttgtcctatacttgat                                                                                                                                <  1:23487‑M1/251‑95 (MQ=255)                                                                                                                               gCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgc                                                                                                               >  1:446945‑M1/1‑139 (MQ=255)                                                                                                                                         gCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgt                                                                                                    <  1:16715‑M1/251‑123 (MQ=255)                                                                                                                                          cGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTTAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacgtgatctcgacgaaatggctgcacctaaatcgtg                                                                                                   <  1:87308‑M1/251‑124 (MQ=255)                                                                                                                                          cGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtg                                                                                                   <  2:435558‑M1/251‑124 (MQ=255)                                                                                                                                                   cGGTGCAGATATGACGGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCATGTACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaactccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaaacgtgatgaaaatc                                                                                          <  1:363974‑M1/251‑133 (MQ=255)                                                                                                                                                       tcaGATATGACTGTAGAAGCCACGCTAACCAAACTGCCTTACCTGCTTAGCCAGGAACTGGATACTGACAACATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctaactttaatcgtgctcccattacgtgactgattctaacaaaacattcacaccaactggcacaattttgtcctaaacttgatctagccgaaatggctgcacctaaatcgtgatgaaaatcacat                                                                                      >  2:341249‑M1/2‑115 (MQ=255)                                                                                                                                                           atatGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacattttt                                                                                  >  1:212534‑M1/1‑111 (MQ=255)                                                                                                                                                               gACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcg                                                                              >  2:163324‑M1/1‑107 (MQ=255)                                                                                                                                                                   ggagaagacaCGCTAACCAACCTGCATTACCTGCTTAGCCAGGCACTGGATACTGAAACCATTCGCAATGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatagtgctcccattaagtgacagattctaacaaaacattaaaaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcaactanatcgtgttgaaaatcacatttttatcgtaat                                                                          >  2:296172‑M1/3‑103 (MQ=255)                                                                                                                                                                     aGAAGCCACGCTAACCAAACTGCATGTCCTNCTTAGCCAGGAACTGGATACTGAAACCATGCGCAAGGCCATGAGCTAAAACCTGCGCGGCGATCTCACGCtttgctcacatctcaccttcatcgtgctcacattacgtgactgattctaacaaaacattaacaccatctggcaaaattttgtcataaacttgacctctacgaaatggctgcacctaaatcgtgatgaaaatcacatctttatcgtaatt                                                                         <  2:284527‑M1/250‑150 (MQ=255)                                                                                                                                                                       aaGCCCCGCTAACCACACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatcgcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgcc                                                                      <  1:296172‑M1/251‑153 (MQ=255)                                                                                                                                                                                  aCCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattc                                                           >  1:575196‑M1/1‑88 (MQ=255)                                                                                                                                                                                       aCTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattcggggc                                                      >  2:175099‑M1/1‑83 (MQ=255)                                                                                                                                                                                          gCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattcggggcgcc                                                   <  1:302080‑M1/251‑172 (MQ=255)                                                                                                                                                                                           cATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattcggggcgccg                                                  <  1:571589‑M1/251‑173 (MQ=255)                                                                                                                                                                                                        gCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcacgttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaacaggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttagcgtaattgccctttaaaattcggggcgccgaccccatgtggtc                                     <  2:111382‑M1/251‑186 (MQ=255)                                                                                                                                                                                                         ccAGGAACTGGGTACTGAAACCATTCGCAAGGCCATGAGCTAAAACCGGCGCGGCGAACTGACGCtttgctcacatctcactttaatcttgctcacattacgtgactgattctaacaatacattaaccccaactggcaaaattttgtcctacacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattcggggcgccgaccccatgtggtct                                    <  2:14829‑M1/251‑187 (MQ=255)                                                                                                                                                                                                             gAACTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattcggggcgccgaccccatgtggtctca                                  >  1:110742‑M1/1‑61 (MQ=255)                                                                                                                                                                                                              aaCTGGATACTGAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattcggggcgccgaccccatgtggtctcaagc                               >  1:308518‑M1/1‑60 (MQ=255)                                                                                                                                                                                                                         gAAACCATTCGCAAGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatgtctgcacctaaatcgtgatgaaaatcgcatttttatcgtaagtgtcctttaaaattcggggcgccgaccccatgtggtctcaagcccaaaggaag                     <  1:341249‑M1/250‑202 (MQ=255)                                                                                                                                                                                                                                     aagtccatgaGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgcctgattctaacaaaacattaacaccaacaggcaaaattttgtcctaaacgtgatctcgacgatatggctgcacctaaaacgtgatgaaaatcacatttttatcgtaattgcccttgaaaattcggggcgccgactcaatgtggtctcaagcccaaaggaagagtgaggcgagtc        <  2:427479‑M1/247‑215 (MQ=255)                                                                                                                                                                                                                                      aGGCCATGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattcggggcgccgaccccatgtggtctcaagcccaaaggaagagtgaggcgagtc        >  1:545009‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                                           aTGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattcggggcgccgaccccatgtggtctcaagcccaaaggaagagtgaggcgagtcagtcgc  >  2:360506‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                                           aTGAGCTAAAACCTGCGCGGCGAACTGACGCtttgctcacatctcactttaatcgtgctcacattacgtgactgattctaacaaaacattaacaccaactggcaaaattttgtcctaaacttgatctcgacgaaatggctgcacctaaatcgtgatgaaaatcacatttttatcgtaattgccctttaaaattcggggcgccgaccccatgtggtctcaagcccaaaggaagagtgaggcgagtcagtcgc  <  1:168758‑M1/251‑221 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     CCACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACAATGTATGTCCGGTAAAGTGAACATGGGTGGTTATGCCACCGGTAACGCCCTCGCCCATGCCGGCGTAATTGGCGGTGCAGATATGACTGTAGAAGCCACGCTAACCAAACTGCATTACCTGCTTAGCCAGGAACTGGATACTGAAACCATTCGCAAGGCCATGAGCCAAAACCTGCGCGGCGAACTGACGCCGGATGATTAAGGAGACTGTAATGCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTG  >  NC_000913/1851601‑1852085

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