Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,571,973:1	+G	7.1%	coding (73/1401 nt)	gadB ←	glutamate decarboxylase B, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,571,973	1	.	G	7.1%	223.7 / 11.9	56	coding (73/1401 nt)	gadB	glutamate decarboxylase B, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base . (23/29);  new base G (2/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATT  >  NC_000913/1571689‑1572256                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            aaCGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑ggta                                                                                                                                                                                                                                                                                         >  1:221446/1‑287 (MQ=255) aaCGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGaa                                                                                                                                                                                                                                                                                         >  1:571840/1‑289 (MQ=255)               ttttCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCAcc                                                                                                                                                                                                                                                                           >  2:417750/1‑289 (MQ=255)               ttttCGGCCCAGGCGCATGCAACAGATCGGCAACCAGATTTACGCAACGCAGGTCGATGGCTGCGGATTGCAGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCAcc                                                                                                                                                                                                                                                                           <  1:617048/289‑1 (MQ=255)                            cgcaTGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGaa                                                                                                                                                                                                                                                                 <  1:466128/286‑1 (MQ=255)                                                       tACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACt                                                                                                                                                                                                                                   <  2:74314/289‑1 (MQ=255)                                                              cgcaGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGcttct                                                                                                                                                                                                                            >  2:51325/1‑289 (MQ=255)                                                                   gTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTAt                                                                                                                                                                                                                         >  2:10750/1‑288 (MQ=255)                                                                   gTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATcc                                                                                                                                                                                                                       >  1:605369/1‑289 (MQ=255)                                                                    tCGTTGGCTGCGGATGTCGTATATTCTTCTTGGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTAt                                                                                                                                                                                                                         <  2:705231/286‑1 (MQ=255)                                                                      gATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCAt                                                                                                                                                                                                                     >  1:203807/1‑289 (MQ=255)                                                                            tgcggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTaaa                                                                                                                                                                                                               <  2:605369/288‑1 (MQ=255)                                                                             gcggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAAATTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACt                                                                                                                                                                                                                                   >  2:209736/1‑267 (MQ=255)                                                                             gcggattgcggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAAATTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACt                                                                                                                                                                                                                                   <  1:209736/267‑1 (MQ=255)                                                                                  ttgcggatATTCTTCTTTGTCGATCCAGTTTTTTTTAATGGATAAATCCATCAATTTGTGGCCATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTa                                                                                                                                                                                                        <  2:418213/289‑1 (MQ=255)                                                                                   tggggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGGTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTaa                                                                                                                                                                                                       <  2:620719/286‑1 (MQ=255)                                                                                      ggatATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTaaaa                                                                                                                                                                                                     <  2:508762/288‑1 (MQ=255)                                                                                               tcttTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATc                                                                                                                                                                                           >  1:346678/1‑289 (MQ=255)                                                                                                   tGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCAt                                                                                                                                                                                       >  1:25980/1‑289 (MQ=255)                                                                                                      cGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTaa                                                                                                                                                                                    <  2:25980/289‑1 (MQ=255)                                                                                                      cGATCCAGTTTTGGTTAATGGATAAATCCATCAATTTGTGGACATTGTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTaa                                                                                                                                                                                    <  2:141153/289‑1 (MQ=255)                                                                                                        aTCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTaaaa                                                                                                                                                                                  >  2:641454/1‑289 (MQ=255)                                                                                                                 ttGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACC‑AAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGAtt                                                                                                                                                                        <  1:275794/289‑1 (MQ=255)                                                                                                                        tGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGa                                                                                                                                                                  >  2:450319/1‑289 (MQ=255)                                                                                                                            tAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATaaa                                                                                                                                                              <  1:417750/289‑1 (MQ=255)                                                                                                                                  cATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGc                                                                                                                                                        <  1:51325/289‑1 (MQ=255)                                                                                                                                    tCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCaa                                                                                                                                                      <  1:450319/289‑1 (MQ=255)                                                                                                                                    tCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCaa                                                                                                                                                      >  1:399390/1‑289 (MQ=255)                                                                                                                                         ttGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCa                                                                                                                                                 <  1:80992/289‑1 (MQ=255)                                                                                                                                                cATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACt                                                                                                                                          <  2:346678/289‑1 (MQ=255)                                                                                                                                                cATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACt                                                                                                                                          <  2:28974/289‑1 (MQ=255)                                                                                                                                                               aGGACTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           <  2:212206/289‑1 (MQ=255)                                                                                                                                                                  tCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           <  2:203807/287‑1 (MQ=255)                                                                                                                                                                      gCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCt                                                                                                                    >  2:316285/1‑289 (MQ=255)                                                                                                                                                                                gCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttattt                                                                                                          >  2:96448/1‑289 (MQ=255)                                                                                                                                                                                       tCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGTGGGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTg                                                                                                    <  1:10750/289‑1 (MQ=255)                                                                                                                                                                                               gCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTa                                                                                            >  2:682782/1‑288 (MQ=255)                                                                                                                                                                                                cGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           <  1:84346/256‑1 (MQ=255)                                                                                                                                                                                                cGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           >  2:84346/1‑256 (MQ=255)                                                                                                                                                                                                cGTTGCCATCAAGATATAAGTGGTCATTGAGAATCTGGAATGCGACAGCGTCGCGCATTTCGTGCAGCGGGAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTATTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATAGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTAc                                                                                           <  2:251675/288‑1 (MQ=255)                                                                                                                                                                                                    gCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGaa                                                                                      <  1:316285/289‑1 (MQ=255)                                                                                                                                                                                                          aaGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATAGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTaa                                                                                                                                 <  2:101848/240‑1 (MQ=255)                                                                                                                                                                                                          aaGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATAGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTaa                                                                                                                                 >  1:101848/1‑240 (MQ=255)                                                                                                                                                                                                                 aaTTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCa                                                                         >  2:28305/1‑289 (MQ=255)                                                                                                                                                                                                                       tCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaaa                                                                   >  1:197976/1‑289 (MQ=255)                                                                                                                                                                                                                        cATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGaaaa                                                                  >  1:369979/1‑289 (MQ=255)                                                                                                                                                                                                                             aTAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTaa                                                                                                                                 >  2:372051/1‑221 (MQ=255)                                                                                                                                                                                                                             aTAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTaa                                                                                                                                 <  1:372051/221‑1 (MQ=255)                                                                                                                                                                                                                                     gAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGaa                                                                                                                                                                                                                                                                 <  2:408534/85‑1 (MQ=255)                                                                                                                                                                                                                                     gAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGaa                                                                                                                                                                                                                                                                 >  1:408534/1‑85 (MQ=255)                                                                                                                                                                                                                                       aTGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGaa                                                   >  1:299358/1‑289 (MQ=255)                                                                                                                                                                                                                                             cATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           <  2:104787/211‑1 (MQ=255)                                                                                                                                                                                                                                             cATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACa                                                                                                                           >  1:104787/1‑211 (MQ=255)                                                                                                                                                                                                                                                        aTTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTa                                  <  1:641454/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                       cGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACt                                                                                                                                          <  2:577945/154‑1 (MQ=255)                                                                                                                                                                                                                                                                                       cGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACt                                                                                                                                          >  1:577945/1‑154 (MQ=255)                                                                                                                                                                                                                                                                                        gATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCAtt  <  1:28305/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            AACGGCCTGACCATTTTTCGGCGCAGGCGCATGCCACAGATCGGCAACCATATTTACGCAACGCAGGTCGATGGCTGCGGATTGCGGATATTCTTCTTTGTCGATCCAGTTTTTGTTAATGGATAAATCCATCAATTTGTGGACATTTTCGTCGTCCCAGGTCTGGCAGAAAGTGGCCAGGTTCTGACGAGCGTTGCCATCAAGATATAATTCGTCATTGATAATCTGGAATGCGACATCGTCGCGCATTTCGTGCAGCGGAAAACGTTTTGATTCTGCGATAGT‑GGAAATAGACTTCGCACCAAAACGTGAATCGAGTAGTTCCGACCTTAAATCCGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATAATATTTTTTCATTGAAAAACAATACAATATGAAATTCTTGGGTGGTGGTAAGGTGTTTTATGCTGTTATTTTTATGCGCATT  >  NC_000913/1571689‑1572256

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