Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,516,393	C→T	G272G (GGC→GGT)	patD →	gamma‑aminobutyraldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,516,393	0	C	T	100.0%	223.1 / NA	61	G272G (GGC→GGT)	patD	gamma‑aminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/27);  total (34/27)

CCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAC  >  NC_000913/1516109‑1516685                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                      ccGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGa                                                                                                                                                                                                                                                                                        >  1:695741/1‑300 (MQ=255)       aCCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGAt                                                                                                                                                                                                                                                                                  <  1:754589/300‑1 (MQ=255)          gcgTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGt                                                                                                                                                                                                                                                                               >  1:433680/1‑300 (MQ=255)            tttGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGAGTAACAGACTGTGTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGTGCTCNATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGTTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTAc                                                                                                                                                                                                                                                                             <  2:196457/299‑1 (MQ=255)            gtttaaattGGCAGAACTGGCGAAAGATATCTTCCCTGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTAc                                                                                                                                                                                                                                                                             <  2:391832/293‑1 (MQ=255)            gTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTAc                                                                                                                                                                                                                                                                             <  1:454069/300‑1 (MQ=255)            gTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTAc                                                                                                                                                                                                                                                                             <  2:473120/300‑1 (MQ=255)                 aGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTTATTAACATACTGTNTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCgg                                                                                                                                                                                                                                                                        <  2:620476/300‑1 (MQ=255)                            cTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGTGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATc                                                                                                                                                                                                                                                             <  2:438055/300‑1 (MQ=255)                                    agctctCTTCCCGGCAGGCGTGATTAACATACTGTTCAGCCGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGTAGTGGTCGAAGGTGTACGTATATTTGGTTCTTCCAATGCTGGACAGGATTGTACTGCGGCTTGTCGGGTCTAcgtgaa                                                                                                                                                                                                                                                     >  2:170546/6‑296 (MQ=255)                                          cTTCCCGGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCNTCAGCCATACCGCGTCGTCCATTAAGNGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAgg                                                                                                                                                                                                                                               >  2:577077/1‑300 (MQ=255)                                                         gATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTCCTGCGGCTTGTCGGATCTACGCGCAAAACGGCCTTTCCTATACGCt                                                                                                                                                                                                                                >  2:342201/1‑300 (MQ=255)                                                         gATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCt                                                                                                                                                                                                                                >  2:495633/1‑300 (MQ=255)                                                         gATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCAAAAGTGCGGATGGTGTAGCTGACGGGCTCTATCGCCACCAGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGCTTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCTAAGGTGTACGTACATTTGGTTATCACAATGCTGGACAGGATTGTACTGCGGCTTGTAGGATCTACGCGCAAAGAGGCATTTACGATActct                                                                                                                                                                                                                                >  2:109421/1‑297 (MQ=255)                                                           ttAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCtgg                                                                                                                                                                                                                              >  2:360851/1‑300 (MQ=255)                                                                   cTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTANGCGTACTCATATGGAACTTGGTGGCANAGCGCCAGTGATTGTTTTTGATGATGTGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑aaaaa                                                                                                                                                                                                                     >  2:579089/1‑300 (MQ=255)                                                                   cTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCAACGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCACAAAGGCATTTACGATACGCTGGTGG‑aaaaa                                                                                                                                                                                                                     >  2:71678/1‑300 (MQ=255)                                                                                     aCGGTGGGTGATCCGCTGACCTGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAATCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCa                                                                                                                                                                                                   <  2:97925/300‑1 (MQ=255)                                                                                            gtgATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTNGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTaa                                                                                                                                                                                            >  1:620124/1‑300 (MQ=255)                                                                                               aTCCGCTGGTCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCGCTATCGCCACCTGCGAGTACATCATCAGCCATACCGCGTCGTCCATTAAGCGTATTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTATTTGTGATGATGCAGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTGGGTTACTACAATGCGGGACAGGATTGTACTGCGCCTTGTCGGATCTACGCACAAAAAGGAATTTACGATACGCTGGTGG‑AAAAACTGGGTGCCGCGGTGGCAACGTTAAAAt                                                                                                                                                                                         <  2:395819/300‑1 (MQ=255)                                                                                                                          ggATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGNGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACgg                                                                                                                                                              <  2:733910/300‑1 (MQ=255)                                                                                                                               gtgtCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCtt                                                                                                                                                         >  1:693985/1‑300 (MQ=255)                                                                                                                               gtgtCGCTGACGGGCTCTATCGCCAACGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGTGCACGTAATTGTTTTTGATGATGTGGATATTGAAGCAGCGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATCCGCTGGTGGCAAAAAC‑GGGTGCTACGGTGGGAACATTAAAATCTGGTGCGACAGATGACGAGTccacaaagcct                                                                                                                                                         >  2:307005/1‑290 (MQ=255)                                                                                                                                           ggCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCTTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGc                                                                                                                                             <  1:488212/300‑1 (MQ=255)                                                                                                                                             ctctATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAAAGTTAAAATCTGGTGCGCCAGGTGACGCGTCTACGGAGCTTGGACCTTTAAGCTc                                                                                                                                           >  2:122435/1‑300 (MQ=255)                                                                                                                                                    gCCACCTGCGAACACATCATCAGCCATACCGCGGCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCTAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGgcg                                                                                                                                    <  1:250798/300‑1 (MQ=255)                                                                                                                                                               gCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGTTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAAcg                                                                                                                         <  2:825584/300‑1 (MQ=255)                                                                                                                                                                  catcatCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGt                                                                                                                      >  2:395649/1‑300 (MQ=255)                                                                                                                                                                        caGCCAGACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCaa                                                                                                                <  2:530735/300‑1 (MQ=255)                                                                                                                                                                          gCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAgg                                                                                                              >  1:99790/1‑300 (MQ=255)                                                                                                                                                                          gCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAgg                                                                                                              >  2:822464/1‑300 (MQ=255)                                                                                                                                                                             aTACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAg                                                                                                           >  2:711851/1‑300 (MQ=255)                                                                                                                                                                             aTACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAg                                                                                                           >  1:167111/1‑300 (MQ=255)                                                                                                                                                                             aTACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGAAGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAg                                                                                                           >  2:3907/1‑300 (MQ=255)                                                                                                                                                                              tACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGt                                                                                                          >  2:378007/1‑300 (MQ=255)                                                                                                                                                                                 cgcgTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAACACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTaga                                                                                                       >  2:58394/1‑300 (MQ=255)                                                                                                                                                                                  gcgTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTagaa                                                                                                      <  1:39895/300‑1 (MQ=255)                                                                                                                                                                                   cctcgtcCATTAAGCGTACTCATATGGAACTTGTTGGCAAATCGCCAGTGATTGTTTTTTGTGATTCGGATATTGAAGCAGTGTTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTagaag                                                                                                     <  2:442684/298‑1 (MQ=255)                                                                                                                                                                                   catcttcCATTAAGCGTACTCATAGGGAACTTATTGGCAAAGCGCCAGTGGTTGTTTTGGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACGTTTGGTTATTACAATGCTGGACAGGATTGTANTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCCACGGAGCTTGGACCTTTAAGCTCGCTTGCGCATCTCGAACGCGTCGGCAAGGCAGTagaag                                                                                                     <  2:207202/295‑1 (MQ=255)                                                                                                                                                                                       gtcCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAGCGCGTCGGCAAGGCAGTAGAAGAGGc                                                                                                 >  2:459738/1‑300 (MQ=255)                                                                                                                                                                                        tcCATTAAGCGTACTCATATTGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTTAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCg                                                                                                <  2:220939/300‑1 (MQ=255)                                                                                                                                                                                          cATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGaa                                                                                              >  2:770934/1‑300 (MQ=255)                                                                                                                                                                                                    aCTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTAAGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGAAAggg                                                                                    >  2:271290/1‑300 (MQ=255)                                                                                                                                                                                                     cTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGc                                                                                   >  2:139281/1‑300 (MQ=255)                                                                                                                                                                                                     cTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGc                                                                                   <  1:740186/300‑1 (MQ=255)                                                                                                                                                                                                         tatGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACAt                                                                               >  2:421299/1‑300 (MQ=255)                                                                                                                                                                                                          agggAACTTTGTGGCGAATCGCCAGTGATTGTTTGCGATTATGCGGATAGTAGAGCATTGTTCTAAAGTGTACGTACATGTGGTTATTACAATGCTGGGCAGGATGGTACTGCGGCTTGTCGGATCTGCGCGCAGAAAGACATTTACGATACGCGGGTGG‑AAAAACTAGGTGCTGCAGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACTGAGCTTGGACCTTTAAGCTCGCTGGCGCAGCGCGAATGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATc                                                                              <  2:110344/298‑1 (MQ=255)                                                                                                                                                                                                                         aaatcgcCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGc                                                               <  1:465176/296‑1 (MQ=255)                                                                                                                                                                                                                                        gTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTAGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAg                                                >  2:328350/1‑300 (MQ=255)                                                                                                                                                                                                                                                       gATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGTACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCCGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTAc                                 <  1:139089/300‑1 (MQ=255)                                                                                                                                                                                                                                                       gATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTAc                                 <  1:137181/300‑1 (MQ=255)                                                                                                                                                                                                                                                       gATAGTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATNTGGTGCGCCAGATGACGAGTCTACGGNGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTAc                                 <  2:583778/300‑1 (MQ=255)                                                                                                                                                                                                                                                         tatTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTCCGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTAGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTAcca                               >  2:439996/1‑298 (MQ=255)                                                                                                                                                                                                                                                            tGAAGCAGTGTTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTNGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATgc                            <  2:740817/300‑1 (MQ=255)                                                                                                                                                                                                                                                                     ggTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgct                   >  1:530878/1‑300 (MQ=255)                                                                                                                                                                                                                                                                     ggTCGAAGGTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTAGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgct                   >  1:237334/1‑300 (MQ=255)                                                                                                                                                                                                                                                                     ggTCGAAGGTGTACGTACATTTGATTATTACAATGCTGGACAGGGTTGTACTGCGGCTTGTGGGGTCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACAGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGGAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACgct                   >  1:341280/1‑300 (MQ=255)                                                                                                                                                                                                                                                                            ggTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTgctggct            <  1:376182/300‑1 (MQ=255)                                                                                                                                                                                                                                                                            ggTGTACGTACATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATAAGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAGAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACNGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTnctggct            >  2:699633/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                 acgtacATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGgcgc       >  1:764451/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                      cATTTGGTTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGGAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAc  <  1:278400/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                      CCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGGCAGAGGCAAAACGGTGGGTGATCCGCTGACCGGTCATCCCAAAGTGCGGATGGTGTCGCTGACGGGCTCTATCGCCACCGGCGAGCACATCATCAGCCATACCGCGTCGTCCATTAAGCGTACTCATATGGAACTTGGTGGCAAAGCGCCAGTGATTGTTTTTGATGATGCGGATATTGAAGCAGTGGTCGAAGGTGTACGTACATTTGGCTATTACAATGCTGGACAGGATTGTACTGCGGCTTGTCGGATCTACGCGCAAAAAGGCATTTACGATACGCTGGTGG‑AAAAACTGGGTGCTGCGGTGGCAACGTTAAAATCTGGTGCGCCAGATGACGAGTCTACGGAGCTTGGACCTTTAAGCTCGCTGGCGCATCTCGAACGCGTCGGCAAGGCAGTAGAAGAGGCGAAAGCGACAGGGCACATCAAAGTGATCACTGGCGGTGAAAAGCGCAAGGGTAATGGCTATTACTATGCGCCGACGCTGCTGGCTGGCGCATTAC  >  NC_000913/1516109‑1516685

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