Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,715,661	C→T	D77N (GAT→AAT)	pdxY ←	pyridoxal kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,715,661	0	C	T	97.8%	141.4 / ‑3.6	45	D77N (GAT→AAT)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (22/22);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

CCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAAC  >  NC_000913/1715418‑1715906                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      ccGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTGGGCGAAATGATATCGCGGGCAGGAAACCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTGTCCGGATGACCCAGTACCTGATCGCAAAAATATGTCGCCTGCGGAGTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑AGTTTCACCCTGCTCCGCCGATCCCAGATATCCACTTAATACGGCATTACAGgtg                                                                                                                                                                                                                                                <  1:332327/251‑1 (MQ=255)    cATGCTCACAGAGAATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTAACGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑CTGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGgtgtgt                                                                                                                                                                                                                                             >  1:261770/1‑251 (MQ=255)       gCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTaa                                                                                                                                                                                                                                           <  1:59357/250‑1 (MQ=255)          cacaGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAAtt                                                                                                                                                                                                                                         <  1:354160/249‑1 (MQ=255)          cacaGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCAGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCTCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAAtt                                                                                                                                                                                                                                         <  1:435/249‑1 (MQ=255)           ccaGAGTATTTCCCGCTCAACCAGCTTTGGCGCAATGATATCGATGGCAGGCAAACCGTGCCGCAAATGAAACTCTGCGACACACGGTGCAAAGATACAGCCTTTTTCCGGATGACCCATTACCGTGTCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTTACCTGACGGACGATACCGAGGCT‑ATGCTCACCCTGCTCCGCCGATCCCAGGTAGCCACTTAATACGGCATTACAGGAGTGTAATTTAt                                                                                                                                                                                                                                      >  2:198453/2‑251 (MQ=255)                             aaCCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGGGCCGCACATGAAACTCTGCGACCCCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCtt                                                                                                                                                                                                                    <  2:334304/251‑1 (MQ=255)                              aCCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTg                                                                                                                                                                                                                   >  1:212768/1‑251 (MQ=255)                              aCCAGATTTGGCGCAATGAAATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACTGCATTACAGGTGTGTAATTTATCAATGGCGGCAATACCTTg                                                                                                                                                                                                                   >  2:387889/1‑251 (MQ=255)                                       ggCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGAGCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTc                                                                                                                                                                                                          <  2:225315/251‑1 (MQ=255)                                         cgcAATGCTATCGCTGGCAGGCAAACCGTGCCGCACATGACACTCTGCGACACCCGGTGCCACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCACCAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑CTGTTCACCCTGCTCCGCCGATCCCCGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCCATGGCGGCCATGCCTTGCACAATTTCgg                                                                                                                                                                                                        >  2:213373/1‑251 (MQ=255)                                              tGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCCCCTGACGGACGATACCTAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTaa                                                                                                                                                                                                    >  2:309049/1‑250 (MQ=255)                                                    cgcgggcaggcaAACCTCTCCGCACATGAAACTCTGCGACACCCTGTGCAACGATACACCCTTTTTCCGGATGCCCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTAGCCTGACGTACGATCCCGAGGAT‑AGGTTTCCCCTGTTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATTGCGGCAATGCCTTGCACAATTTCGGTTAAAt                                                                                                                                                                                                  <  1:289060/242‑1 (MQ=255)                                                       tGGCAGGCAAACCGTGGCGCACATGAAACTCTGCGACACCCGGTGCCACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCTCAAAACTCTTTCGTCTGCGGATTCGCGGCCTTCACCTGCCGGCCGTTACC‑AGGATGATGTTCAACCTGCTCCGCCAATCCAAGCAACCCACTTTAGACGGCATTGCCGATGTGTAATTTATCACAGTCGGCAATGTCGTGATCAATTTCGCTGAAATGGCTCggc                                                                                                                                                                                          >  2:319797/1‑251 (MQ=255)                                                        ggcaggcaACCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATCCCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggc                                                                                                                                                                                          <  2:542730/250‑1 (MQ=255)                                                           aggcaAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGCCGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGCTAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggcggcc                                                                                                                                                                                      >  1:580367/1‑250 (MQ=255)                                                               aaaCCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCAc                                                                                                                                                                                  >  1:322638/1‑251 (MQ=255)                                                                         gCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTACTTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCCTCACGCAGCCAGTc                                                                                                                                                                        >  2:212857/1‑251 (MQ=255)                                                                           acaTGCAACTCTGCCACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCCTTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGAAGGACGATACCGAGGGT‑ATGCTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCATTGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACTCAGCCAGTcc                                                                                                                                                                       >  2:244131/1‑250 (MQ=255)                                                                                     cggCGACACCCGGTGCAACGATCCAGCCTTTTTCCGGAGGACCCATTCCCGGATCGAAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTAt                                                                                                                                                            <  2:322638/249‑1 (MQ=255)                                                                                         gacaaacgccgcaccGATACAGCATTTTTCCGGCTGACCCAATCACGCATCGCAAAACTCTTTCACCTGCGGATTCGCGGCTTTCACCCGACGGACGATCCCGCGGAC‑ATGTTCACCATGCTCCGCCGATCCCAGATAGCCACTTCATACGGCCTTACAGGTGCGTAATTTATCAATGGCGGCAATGACTTGCACAATTTCGGTTAAATGGCTGGGCGGCCTCCCGCAGCCCGTCCCTTTGCCATATTggg                                                                                                                                                        >  2:346180/15‑251 (MQ=255)                                                                                                   gCAACGATACAGCCGTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTGCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAg                                                                                                                                               <  2:99005/250‑1 (MQ=255)                                                                                                       cGATACAGCCTTTTTCCGGATGACCCATTCCCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGGCGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGAGAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAt                                                                                                                                          <  1:63464/251‑1 (MQ=255)                                                                                                         aTACATCCTTTTCCCGGATGACCCTTTACCGTATCGCAAAAATATTTCGCCGGCGGATTCGCGGTTTTCACCTGCCGGACGATCCCGAGGAT‑ATGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTg                                                                                                                                        <  2:491492/251‑1 (MQ=255)                                                                                                                   tttCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTc                                                                                                                                                                                                          <  1:111431/175‑1 (MQ=255)                                                                                                                   tttCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTc                                                                                                                                                                                                          >  2:111431/1‑175 (MQ=255)                                                                                                                                    cccGGATCGCAAAAATATTTCGCCTGCGGTTTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGGCGTTcgc                                                                                                             <  2:538493/250‑1 (MQ=255)                                                                                                                                           cGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGCCGGACGATACCGCGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGAATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTAGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGCACGGTGTTCAGCGGCCAGACGTTCGCTCCCAgg                                                                                                      >  2:8491/1‑251 (MQ=255)                                                                                                                                            gCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAggc                                                                                                     >  2:478683/1‑251 (MQ=255)                                                                                                                                                              cGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTg                                                                                   >  1:242110/1‑251 (MQ=255)                                                                                                                                                                                      gacgATACCGAGGTT‑ATGTTCACCCTGCTCCGCCGATCCCAGCTAGCCACTTAATACGGCATTACAGGTGTGTAATTGATCAATGGCTGCAATGCCTTGCACAATGTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTTATTAGAAAATTGACCGGTGTTCAGCGGCCAGACGTTCGCGCCCATGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCAt                                                           <  1:577928/251‑1 (MQ=255)                                                                                                                                                                                       acgATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCAt                                                           >  2:461288/1‑250 (MQ=255)                                                                                                                                                                                            aCCGAGGAT‑ATGGTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAaacaa                                                     <  1:387889/251‑1 (MQ=255)                                                                                                                                                                                             ccggggAT‑AGGTTCCCCCTGTTCCACCGATCCCAGCTAGCCACGTAATAGGGCATTACAGGTGTGTAATGTAGCAATGGCGGCAATGCCTTTCACAATTTCGGTTAAAGGGCTGGGCGGCCTCACGCAGCCACTCCATTGGCCGTGTGGGGTGTGATTAGAAAATTGAACGTTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGACGCATCGTAAACTCTGCCGCGGTGTTACCCGCATGAGCATAaacaac                                                    <  2:243629/247‑1 (MQ=255)                                                                                                                                                                                              cGAGGAT‑ATGGTCCCCCTGCTCAGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGTGCGGCGTCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGTAACTTGAACGGTGTCCAGCGGCCAGTCGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGACCTGTTCCCCGCATGACCATAACCAACg                                                   >  2:402718/1‑251 (MQ=255)                                                                                                                                                                                               gAGGAT‑ATGTTCACCCTGCTCCGCCGCTCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGt                                                  >  2:245164/1‑251 (MQ=255)                                                                                                                                                                                               gAGGAT‑AGGTTCACCCGTCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAAGGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAGTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGt                                                  <  2:580367/251‑1 (MQ=255)                                                                                                                                                                                                              ccTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACGCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCg                                    <  1:245164/251‑1 (MQ=255)                                                                                                                                                                                                                         gATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAAAGTGAGACTGGATCGCGAGAATAttttt                         >  2:322394/1‑251 (MQ=255)                                                                                                                                                                                                                           tCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACACTTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGTGAATATTTTtc                        >  2:92199/1‑250 (MQ=255)                                                                                                                                                                                                                                      cACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCt            >  1:51029/1‑251 (MQ=255)                                                                                                                                                                                                                                       aCTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTg           <  1:443177/251‑1 (MQ=255)                                                                                                                                                                                                                                          tAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGt          <  1:383626/249‑1 (MQ=255)                                                                                                                                                                                                                                              aCGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATaaa     <  1:268086/250‑1 (MQ=255)                                                                                                                                                                                                                                                ggCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCGGCCGCACTGTTACCCGCATGACCATAAACAACGGGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAAc  <  1:213373/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGAT‑ATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAAC  >  NC_000913/1715418‑1715906

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