Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,793,897	Δ1 bp	100%	coding (213/552 nt)	btuE ←	glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,897	0	T	.	100.0%	298.2 / NA	64	coding (213/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (36/28);  total (36/28)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCG  >  NC_000913/1793753‑1794041                                                                                                                                                 |                                                                                                                                                   gACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgttc                                                                                                                                                <  1:904312/149‑5 (MQ=255)    tGCGGTAGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAg                                                                                                                                            >  2:1113973/1‑149 (MQ=255)        gTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCATTCCTGGGTTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAAc                                                                                                                                        >  1:838455/1‑149 (MQ=255)         tCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACt                                                                                                                                       <  2:771207/149‑1 (MQ=255)           ggCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCGTGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCCGAAACTgg                                                                                                                                     >  2:1130913/1‑149 (MQ=255)            gcgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGt                                                                                                                                    <  2:327166/149‑1 (MQ=255)            gcgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGt                                                                                                                                    >  1:333015/1‑149 (MQ=255)             cgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTGCTTCGACATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGtt                                                                                                                                   <  2:810373/149‑1 (MQ=255)             cgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGtt                                                                                                                                   <  1:497257/149‑1 (MQ=255)                 ggggCAATCAATTTTTGATCCAGCGGACGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACCTCGGGAACGTCACCCCCCATGTTGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCAc                                                                                                                               <  1:265956/147‑1 (MQ=255)                        tcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATc                                                                                                                       >  1:465141/1‑149 (MQ=255)                                 gATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTAATGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GTAATCCCAGCACCa                                                                                                              >  1:826757/1‑149 (MQ=255)                                  aTACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGTGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCAt                                                                                                             <  2:1118206/149‑1 (MQ=255)                                    aCAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATaa                                                                                                           >  1:1062738/1‑149 (MQ=255)                                      accGGATGGCGTCCTTCGCCATTAACTTTATTCTTACTGAACATCGGGAACGTCACCCCCCATGTGTTGGTACAGTAAGTTTTAATCCCTTCATCGCTGCCCGGGTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATaaaa                                                                                                         <  2:125685/147‑1 (MQ=255)                                      aGCGGTTGGCGTCCTTCGCCATTAACATCAAGGCTACTGAACATAGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑AGTTGCAGAAACTGGTTGCACG‑GGAATCCCCGCACCATaaaa                                                                                                         >  2:535852/1‑149 (MQ=255)                                      aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATaaaa                                                                                                         >  2:821940/1‑149 (MQ=255)                                         ggATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCt                                                                                                      <  2:887073/149‑1 (MQ=255)                                            tGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCtcga                                                                                                   <  2:295174/149‑1 (MQ=255)                                             ggCGTCCTTCGCCATTAACTTCAATCTTGCTGCACATCGGGCACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCtcgat                                                                                                  <  2:129703/149‑1 (MQ=255)                                                  ccctCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGAcc                                                                                             <  1:310719/146‑1 (MQ=255)                                                     tCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCCCGATCGACCCAg                                                                                          >  2:837048/1‑149 (MQ=255)                                                        ccATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCt                                                                                       >  1:1059464/1‑149 (MQ=255)                                                           ttAACTTGAATCTTACTGAACATCGGGAACGTCCCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGGATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTc                                                                                    <  2:338451/149‑1 (MQ=255)                                                               cTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                                                                >  1:997462/1‑149 (MQ=255)                                                                         aCTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCaa                                                                      <  1:1109296/149‑1 (MQ=255)                                                                         aCTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCaa                                                                      <  2:880739/149‑1 (MQ=255)                                                                          cTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAAc                                                                     >  1:967220/1‑149 (MQ=255)                                                                           tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACt                                                                    <  2:261462/149‑1 (MQ=255)                                                                           tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACt                                                                    >  1:907590/1‑149 (MQ=255)                                                                           tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCAGCCCGGGTC‑TGTTCCAGAAACCGGTTGCACG‑GGAAACCCAGCACCATACAACCTCGATCGACCCGGGCTTTCTGAATCTACTCCAACt                                                                    >  2:119509/1‑149 (MQ=255)                                                                             aaCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGc                                                                  <  1:821940/149‑1 (MQ=255)                                                                               cATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTc                                                                >  1:939984/1‑149 (MQ=255)                                                                                aTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCa                                                               >  1:1110084/1‑149 (MQ=255)                                                                                      caCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTgc                                                        <  2:1032828/148‑1 (MQ=255)                                                                                        cGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTgcgg                                                      >  2:188119/1‑149 (MQ=255)                                                                                         gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTCCCCATA‑TTgcggc                                                     >  1:1026396/1‑149 (MQ=255)                                                                                         gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TCATCCCG‑AGTTTGATCCTCGCAGAGTCCCCGGACCATAAAACAACGATCGGCGTAGTGTCTCTGATTATTGTCCAACTGCTCATA‑GTgcgtc                                                     >  1:844207/1‑147 (MQ=255)                                                                                                cccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGc                                              <  2:967220/149‑1 (MQ=255)                                                                                                  cATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGccc                                            >  2:64014/1‑148 (MQ=255)                                                                                                   aTGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTACCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGAGGCGTTAAGCcac                                           >  1:824858/1‑149 (MQ=255)                                                                                                          ggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGa                                    <  1:1113973/149‑1 (MQ=255)                                                                                                           gtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCTTA‑TTGCGGCGTTAAGCCACACTTTGAg                                   >  2:721443/1‑149 (MQ=255)                                                                                                            tACATTAAGTTTTCATCTCTTCATCGCCGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGGTCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAgg                                  <  2:480182/149‑1 (MQ=255)                                                                                                             aCAGTAAGTTGTAATCTCCTCAGCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATACAACCTCGATCGTCCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGc                                 <  2:173173/149‑1 (MQ=255)                                                                                                             aCAGGAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGc                                 <  1:185020/149‑1 (MQ=255)                                                                                                                 tAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGTTTTCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGTCGACa                             >  1:1135531/1‑149 (MQ=255)                                                                                                                 tAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACa                             >  1:112598/1‑149 (MQ=255)                                                                                                                      tttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGTTATTgtc                        >  2:805454/1‑147 (MQ=255)                                                                                                                         aaTCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAAt                     >  1:567705/1‑149 (MQ=255)                                                                                                                            ctctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa                  <  1:1136619/149‑1 (MQ=255)                                                                                                                            ctctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa                  <  1:973369/149‑1 (MQ=255)                                                                                                                            cgctTCATCGCTGCCCGGTGC‑CGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGCTCGACCCAGGCTTTCTTAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa                  <  1:582749/147‑1 (MQ=255)                                                                                                                             tctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGAcaaccaac                 >  2:13771/1‑149 (MQ=255)                                                                                                                             tctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAc                 >  2:1012995/1‑149 (MQ=255)                                                                                                                             tctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAc                 >  1:155238/1‑149 (MQ=255)                                                                                                                              ctTCATCGCTGCCCGGTTC‑TGTTACAGAAACTGGGTGCACG‑GGAATCCCAGCAGCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa                >  1:635386/1‑149 (MQ=255)                                                                                                                                tCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCGCGATCGACCCAGGCTTTCTGAATATTCTCCGACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACACTCAACAGc              >  2:314822/1‑149 (MQ=255)                                                                                                                                 cATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCCCCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGca             >  1:580786/1‑149 (MQ=255)                                                                                                                                 cATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGca             <  2:838455/149‑1 (MQ=255)                                                                                                                                  aTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGATGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGGCAATCCACACTTTGAGGCGACATTGACAATCAACAGcac            >  1:601930/1‑149 (MQ=255)                                                                                                                                           cGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGc   <  1:866276/149‑1 (MQ=255)                                                                                                                                           cGGTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACCGCT‑ATACTTAGGGCTAGAATCCACACGTTGGGGCGACATGGACAATCAACAGCACATTACCGGc   >  2:812993/1‑149 (MQ=255)                                                                                                                                            ggTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGCCACTCAACAGCACATTACCGGCg  >  1:725624/1‑149 (MQ=255)                                                                                                                                            ggTTC‑TGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCg  <  2:731605/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACG‑GGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATA‑TTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCG  >  NC_000913/1793753‑1794041

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