Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,667,023	C→T	Q186Q (CAG→CAA)	gadA ←	glutamate decarboxylase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,667,023	0	C	T	100.0%	90.2 / NA	44	Q186Q (CAG→CAA)	gadA	glutamate decarboxylase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/22);  total (22/22)

AGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAACTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATT  >  NC_000913/3666736‑3667238                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                        aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGgcgc                                                                                                                                                                                                             >  1:199759/1‑300 (MQ=17) aGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGgcgc                                                                                                                                                                                                             >  2:674766/1‑300 (MQ=17)           gACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGAtct                                                                                                                                                                                                   >  1:757682/1‑300 (MQ=17)                          cTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCcaaat                                                                                                                                                                                    >  1:762928/1‑297 (MQ=17)                          cTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCANGTCGATGTCGATACCGGTGTCGGCCTGGAATTNATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTAACGCAGCTCCACAt                                                                                                                                                                                    >  2:582824/1‑300 (MQ=17)                           tGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATTTCGATACCGGTGTCGGCCGGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCTGGGCGCATAGGGATCTCACGCAGCTCcac                                                                                                                                                                                      <  2:302956/297‑1 (MQ=12)                           tGAAGCACTGAGCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACAt                                                                                                                                                                                    <  2:757685/299‑1 (MQ=17)                           tGAAGCACGGATCGACTTCACACGCGGCAGGTGGAAGTCCCAGACGATATCCGGGGCGACGAACGGGGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTTTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGACACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACAt                                                                                                                                                                                    <  2:398447/299‑1 (MQ=255)                                   tGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAg                                                                                                                                                                            <  1:673358/299‑1 (MQ=17)                                                  gcgGCAGGCGGAAGTCCCAGACGATATCCGGTGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGcc                                                                                                                                                            <  1:475373/300‑1 (MQ=17)                                                       aGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCcagca                                                                                                                                                        >  1:763241/1‑299 (MQ=17)                                                               gTCCCAGACGATATCCTGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTAc                                                                                                                                               <  1:701392/300‑1 (MQ=17)                                                               gTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCCATCCTACGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGAGGGCGACTTTATGCAAGCAGGTTTGt                                                                                                                                                 >  2:80714/1‑298 (MQ=11)                                                               gTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTAc                                                                                                                                               <  1:638220/300‑1 (MQ=17)                                                               gTCCAAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTAc                                                                                                                                               >  1:418397/1‑300 (MQ=17)                                                                tCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGTATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTAcc                                                                                                                                              >  2:53497/1‑300 (MQ=17)                                                                    aGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGAc                                                                                                                                          <  2:57773/300‑1 (MQ=17)                                                                     gACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGAcc                                                                                                                                         >  1:674882/1‑300 (MQ=17)                                                                         atatCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGNCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGcac                                                                                                                                     <  1:562359/300‑1 (MQ=17)                                                                               ggggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCTGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAg                                                                                                                                <  1:218695/299‑1 (MQ=17)                                                                               ggggCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAg                                                                                                                                <  2:715495/299‑1 (MQ=17)                                                                                            gTGCCAGGAAGCCACCNCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGNTTTGTACCGGACCGCACANCAGGTTTGGTTTATcc                                                                                                                   >  1:669888/1‑299 (MQ=17)                                                                                            gTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGGGTTTCGGGTCCCTAAACAATTGAACGGGGCGCATAGGGATCTCACGCAGCTCCAAATCCCAGTAGCGGGCGAATTTATGCAGGCAGATTTGTAACGGAACGCACACCAGGTTTGGTTTAt                                                                                                                     >  2:59945/1‑297 (MQ=11)                                                                                            gTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATcc                                                                                                                   >  1:84271/1‑299 (MQ=17)                                                                                                cAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGTAGATTTGTACCGGACCGCACACTAGGTTTGGTTTATCCGtt                                                                                                                >  2:170419/1‑298 (MQ=17)                                                                                                          ccGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCa                                                                                                    >  1:177873/1‑300 (MQ=17)                                                                                                           cGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCa                                                                                                    >  1:592590/1‑299 (MQ=17)                                                                                                                      tCGATGTGCATTTCGATGTCGATACCGGGGTCGGCCTTGAATTTATCCGTCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCTGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGc                                                                                        <  1:210288/300‑1 (MQ=255)                                                                                                                      tCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATNAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGc                                                                                        <  2:557468/300‑1 (MQ=17)                                                                                                                      tCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGc                                                                                        <  1:207982/300‑1 (MQ=17)                                                                                                                               atgtcgatgtcgATACCGGTGTTGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGTACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTgcgc                                                                                  <  2:35520/297‑1 (MQ=17)                                                                                                                                 gtcgatgtcgatACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAAc                                                                              >  1:229187/1‑299 (MQ=17)                                                                                                                                                   gtCGACCTGGAATTTATCCAGCACATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATcc                                                            <  2:638179/299‑1 (MQ=14)                                                                                                                                                      ggCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCccgc                                                         <  1:545559/299‑1 (MQ=14)                                                                                                                                                                 tATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAgg                                             <  2:52147/300‑1 (MQ=14)                                                                                                                                                                         gcATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTCCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCNGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCCTACGTTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGgaa                                     >  2:252126/1‑300 (MQ=9)                                                                                                                                                                              gTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATCGGGATCTCACGCAGCTCCACATCCAAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCATCGGAAGAAc                                 >  2:275760/1‑299 (MQ=9)                                                                                                                                                                                     ggTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCAAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGATTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGgt                          >  2:289430/1‑299 (MQ=14)                                                                                                                                                                                      gTTGTGTGGACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACTACGCCGATGGTTTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAAAAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGTCGAATTTATGCCAGCAGATTTGTACCGTACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGgtgt                        <  2:680302/300‑1 (MQ=11)                                                                                                                                                                                      gTTGTGGGAACTCATAGTTACCGGTGTAGTTCACGCCGAAAGTCTGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGgtgt                        <  2:680289/300‑1 (MQ=14)                                                                                                                                                                                      gTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCAAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCAGCCGAGCATACAGGCCTCGGAAGAATCAATGgtgt                        >  2:36052/1‑300 (MQ=9)                                                                                                                                                                                           gggAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTg                   <  1:641686/300‑1 (MQ=14)                                                                                                                                                                                                 tCATATTTGCCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGGTGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCTGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAAc              <  2:740131/299‑1 (MQ=9)                                                                                                                                                                                                            ggTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAATTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCAtt  >  1:541706/1‑300 (MQ=11)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                        AGCGGAGCCAGACCGAATTTATGGCCTGAAGCACTGATCGATTTCACACGCGGCAGGCGGAAGTCCCAGACGATATCCGGGGCGACGAACGGTGCCAGGAAGCCACCGCTGGCAGCGTCGATGTGCATGTCGATGTCGATACCGGTGTCGGCCTGGAATTTATCCAGCGCATCGTGCAGCGGTTGTGGGAACTCATAGTTACCGGTGTAGGTCACGCCGAAAGTCGGCACCACGCCGATGGTGTTTTCGTCACAGGCTTCAATCATGCGTTTCGGGTCCATAAACAACTGACCGGGGCGCATAGGGATCTCACGCAGCTCCACATCCCAGTAGCGGGCGAATTTATGCCAGCAGATTTGTACCGGACCGCACACCAGGTTTGGTTTATCCGTTGGTTTGCCTGCAGCTTCCATACGCTTGCGCCAACGCCATTTCATCGCCATCCCGCCGAGCATACAGGCCTCGGAAGAACCAATGGTGTTGGTGCCAACGGCCTGACCATT  >  NC_000913/3666736‑3667238

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