Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	789,580	C→A	E134* (GAA→TAA)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,580	0	C	A	94.9%	103.5 / 0.2	40	E134* (GAA→TAA)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (18/19);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

ACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGT  >  NC_000913/789344‑789826                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       aCTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGgaa                                                                                                                                                                                                                                           <  1:15490/250‑1 (MQ=255)                     ccACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGc                                                                                                                                                                                                                      >  1:268884/1‑251 (MQ=255)                      cacaCCTTTGGGAATGGAAACTGCTTTGGACCCCAGTGAGCGGCAATCGGTCAGCAAGGCATGTTCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAACCTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGAGTTACAGTGAAGAGGTAGACCTTAACCCGGCACCCTGCg                                                                                                                                                                                                                     >  2:162333/1‑251 (MQ=255)                         cccTTTGGGCATGGAAACTTCTTTGGTCCCCAGTTAGCGGAAATCGATCAGCAAGGCATGATCTTTCTGGCCGAGCGCGGAAATTAGTTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGGCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGca                                                                                                                                                                                                                  <  2:381008/250‑1 (MQ=255)                                 gCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCCGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAAcacgaggtaggaagagcgtagtgtt                                                                                                                                                                                                          >  2:69962/1‑227 (MQ=255)                                 gCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCg                                                                                                                                                                                                                               <  1:69962/230‑1 (MQ=255)                                  cATGGAAACTGCTTTGGTCCCCAGCGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTACCGAGCGCGGAAATTAGCTGATACATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCATCGCCAGATACTAAAGCTGCTGCAATCCGGTTCCGACCGCGACTTACAGTGAAGCGGAAGCACTTAACCCGGCACCCGGCGGCACAtatccgc                                                                                                                                                                                                         >  2:447717/1‑245 (MQ=255)                                      gaaaatggttgggacaccAGTGAGCGGCAATCGATAAGCAAGGCATGATCTTTCTTGACGCGCGCGGAAATTCGCTGATCCATGATCCCGCAGTTACAGCCTACACACTGGTTTTCTGCGTCCTCACCGTTAAGCGCGATTTGTGCGCCGTCCAGCCGCAGATGATAAAGAAGCTGCAATCTGGTTCGGACCGCGCCTCTCAGTGAAACGGAAGAAATTAACCCGGCACGCTGAGTACCATTGCCGCTgcg                                                                                                                                                                                                     >  2:430126/17‑249 (MQ=255)                                              tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGt                                                                                                                                                                                             >  1:369155/1‑251 (MQ=255)                                                    ccccAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAACTTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCAcgc                                                                                                                                                                                       >  1:448546/1‑251 (MQ=255)                                                             gCGGCAATCGATCAGCAAGGCATGTTCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGc                                                                                                                                                                              <  1:162333/251‑1 (MQ=255)                                                                     cGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCTCCATGTCCACGCCGCCGAAGCTgttgtt                                                                                                                                                                       <  1:402945/250‑1 (MQ=255)                                                                             aaGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTg                                                                                                                                                              >  1:55595/1‑251 (MQ=255)                                                                                gCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGc                                                                                                                                                             >  1:422486/1‑249 (MQ=255)                                                                                  aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGAt                                                                                                                                                         <  1:134077/251‑1 (MQ=255)                                                                                         ttcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTca                                                                                                                                                   <  1:269446/250‑1 (MQ=255)                                                                                             tgacGAGCGCGGAAATTAGCTGATCCAGGTACCCGCAGTTACAGCCTTCAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTcacca                                                                                                                                                <  2:210025/246‑1 (MQ=255)                                                                                                        gAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTa                                                                                                                                   <  1:137610/251‑1 (MQ=255)                                                                                                               gCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATCCGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAg                                                                                                                                                           >  1:53679/1‑220 (MQ=255)                                                                                                                  gATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATTTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGGACGTAGTTAGCCCAt                                                                                                                         <  1:410134/251‑1 (MQ=255)                                                                                                                     ccATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGa                                                                                                                      <  1:35742/251‑1 (MQ=255)                                                                                                                                   ttACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGc                                                                                                        >  1:450199/1‑251 (MQ=255)                                                                                                                                                   tggtttCTGCTTCCTGACTTTTAAGCTCGTTTTGGTCTCCGTTCTGCGTCAGATTAGAAAGTTGCTGCAATACGGTCCCGTCCGCGTCTTACAGTGAAGCGGAAGAACTTAACCTGGCACCCTGCGGCACATTGCCGAGGATCACCATGTCCACGCCGCCGAAGTTGTTGTTACGCAGTTGCAGATGTTTCACCACGTCACGAACGTAGTTATCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATc                                                                                         <  2:263562/247‑1 (MQ=255)                                                                                                                                                     tttacTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCCGATGATCAAGATGCTGCAATCCGGTTCCGACCGCGACTTACAGTGAAGTGGAAGAACTTAACACAGCACCCTGCAGCACATTTCCGCTGAGCACAATGGCCACGCCAACGAGGCTGTAGTTAAGCAGTTGCACATGATTCAACACACCCCACACGTAGTTAGCACATTAATAGTTTTCATTCGCGACacagggagactcacg                                                                                      >  2:267459/5‑237 (MQ=255)                                                                                                                                                      tttCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAgg                                                                                     >  1:174450/1‑251 (MQ=255)                                                                                                                                                      ttacTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAg                                                                                      <  1:441092/247‑1 (MQ=255)                                                                                                                                                        tCTGCTGCCTGAGCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACTGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTTCGGCACATTGCCGCTGATCTCCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAggg                                                                                    <  1:375524/250‑1 (MQ=255)                                                                                                                                                                cTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACtcgt                                                                           >  2:195401/1‑251 (MQ=255)                                                                                                                                                                  gACCGTTAAGCGCGATTTGGGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCAGCGAGGGAAAACtcgtcg                                                                         <  1:134368/251‑1 (MQ=255)                                                                                                                                                                  gACCCTTAAGCGCGATTTGGTCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTCCCGACCGCGACTTACAGTGAAGCGGAAGATCTTATCCCGTCACCCTGCGTCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCGGTTGTTACGCAGTTGCAGATGTTGCCCCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGTGCATCGAGGGAAAACtcgtcg                                                                         <  2:395182/251‑1 (MQ=255)                                                                                                                                                                        ttaGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATggtcg                                                                                             <  2:312927/223‑4 (MQ=255)                                                                                                                                                                        tAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATggtcg                                                                                             >  1:312927/1‑222 (MQ=255)                                                                                                                                                                                    gtgCGCCGTCCAGCGGCAGAGGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGATCTTACCCCGGCACCCTGCGGCACTTTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGCTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCg                                                       <  1:330220/251‑1 (MQ=255)                                                                                                                                                                                       cgcCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGACCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGAAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGcg                                                    >  1:117017/1‑250 (MQ=255)                                                                                                                                                                                                        tgagaaaGCTGCTGCAATACGCTTCCGACCGCGACTTACAGTGAAGCTGAAGAACTTAACCCGGCAACCTGCGGCCCATGTCCGCTGATCTCCATGTACACGCCGCCGAAGCTGTGGTTACGCAGTTGCAGATGTTGCTCCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACtt                                    <  1:165040/246‑1 (MQ=255)                                                                                                                                                                                                               gctgctGCAATACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACgccgcc                                                                                                                                                                                   >  1:40294/1‑100 (MQ=255)                                                                                                                                                                                                               gctgcgGCAATACGGTTCCGACCGCGACTGACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGTACAGTGCCGCTGATCACCATGTCCACgccgcc                                                                                                                                                                                   <  2:40294/100‑1 (MQ=255)                                                                                                                                                                                                                       aaTACGGTTCCGACCGCGACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCgtgg                    >  1:20487/1‑251 (MQ=255)                                                                                                                                                                                                                                     cgcgACTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCa      <  1:399369/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCTTACAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATccaggt  >  1:67995/1‑246 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       ACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGCACATTGCCGCTGATCACCATGTCCACGCCGCCGAAGCTGTTGTTACGCAGTTGCAGATGTTTCACCACGCCACGAACGTAGTTAGCCCATTGATAGTTTTCATGTGCGACAATGGGCGCATCGAGGGAAAACTCGTCGAGCTGATTTTCATAATCGGCTGCCATCACGCGAACTTTACGGTCATCGCGTGGTGCACAACTGATCACGGT  >  NC_000913/789344‑789826

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