Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,013,607	T→A	V121E (GTG→GAG)	SAUSA300_RS04960 →	Serine protease HtrA‑like

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,013,607	0	T	A	100.0%	257.2 / NA	87	V121E (GTG→GAG)	SAUSA300_RS04960	Serine protease HtrA‑like
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (87/0);  total (87/0)

AATCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGTGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAG  >  NC_007793/1013510‑1013702                                                                                                  |                                                                                                aaTCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAgaga                                                                                               >  1:1233748/1‑100 (MQ=255) aaTCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAgaga                                                                                               >  1:3766103/1‑100 (MQ=255) aaTCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAgaga                                                                                               >  1:4131495/1‑100 (MQ=255)   tCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGaaa                                                                                             >  1:5003894/1‑100 (MQ=255)        aTAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGGTACAgagatcggaag                                                                                        >  1:1841616/1‑93 (MQ=255)                 aaGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGATAAATCAAAATATTATGa                                                                               >  1:39789/2‑100 (MQ=255)                  aGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAg                                                                              >  1:5588312/1‑100 (MQ=255)                   gtgtataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGc                                                                             >  1:4903772/4‑100 (MQ=255)                   gCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGc                                                                             >  1:1645452/1‑100 (MQ=255)                       ataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCaaaa                                                                         >  1:3999971/1‑100 (MQ=255)                       ataACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCaaaa                                                                         >  1:3816553/1‑100 (MQ=255)                           cAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCt                                                                     >  1:4931416/1‑100 (MQ=255)                                   cttattTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGAc                                                             >  1:2107569/2‑100 (MQ=255)                                         tcaCAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATc                                                       >  1:612047/3‑100 (MQ=255)                                         tAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTacc                                                       >  1:5317065/1‑98 (MQ=255)                                           aCAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCaa                                                     >  1:4927868/1‑100 (MQ=255)                                            cAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAAc                                                    >  1:2858078/1‑100 (MQ=255)                                                      aTCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAAcc                                          >  1:1096683/1‑100 (MQ=255)                                                      aTCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAAcc                                          >  1:1990977/1‑100 (MQ=255)                                                          aaaaaaCAAAATTCATTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATa                                      >  1:455789/1‑100 (MQ=255)                                                                        cttAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACt                        >  1:2284989/2‑100 (MQ=255)                                                                        aTTAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACt                        >  1:1879921/1‑100 (MQ=255)                                                                         ttAGATTCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTg                       >  1:3176184/1‑100 (MQ=255)                                                                           ccccccAGTGGACCAATATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGaa                     >  1:3933780/6‑100 (MQ=255)                                                                               tCAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATga                 >  1:1454412/1‑100 (MQ=255)                                                                                cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:2592393/1‑100 (MQ=255)                                                                                cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:316217/1‑100 (MQ=255)                                                                                cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:4238382/1‑100 (MQ=255)                                                                                cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:2514668/1‑100 (MQ=255)                                                                                cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:1875066/1‑100 (MQ=255)                                                                                cAGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgag                >  1:1819836/1‑100 (MQ=255)                                                                                 tgTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:5537977/2‑100 (MQ=255)                                                                                 cgTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:3527248/2‑100 (MQ=255)                                                                                 actGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:818922/3‑100 (MQ=255)                                                                                 aGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:912643/1‑100 (MQ=255)                                                                                 aGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:5263707/1‑100 (MQ=255)                                                                                 aGTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATgaga               >  1:357851/1‑100 (MQ=255)                                                                                  ttGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaa              >  1:1030114/2‑100 (MQ=255)                                                                                  gTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaa              >  1:2572016/1‑100 (MQ=255)                                                                                  gTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaa              >  1:552210/1‑100 (MQ=255)                                                                                  gTGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaa              >  1:3644582/1‑100 (MQ=255)                                                                                   tGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGGTAAAGTAGAATCAACTGAAATGAGaaa             >  1:5117514/1‑100 (MQ=255)                                                                                   tGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaaa             >  1:1779416/1‑100 (MQ=255)                                                                                   tGGACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGaaa             >  1:3708299/1‑100 (MQ=255)                                                                                    tgACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:2469942/2‑100 (MQ=255)                                                                                    tgACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:1523715/2‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGaaagatcggaa            >  1:128103/1‑92 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:4466892/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:3844390/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:625696/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:4039756/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:619050/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:4998515/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:3204015/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:5264748/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:3454985/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:118851/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:2401622/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:2834191/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:2847849/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:285750/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:2873998/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:3405650/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:3389580/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGAGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:1029043/1‑100 (MQ=255)                                                                                    ggACCAAGATACAGAGAAACCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:1065626/1‑100 (MQ=255)                                                                                    ataCCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAg            >  1:2579539/3‑100 (MQ=255)                                                                                     atccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGc           >  1:5692649/3‑100 (MQ=255)                                                                                      tgcAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:2925862/3‑100 (MQ=255)                                                                                      tccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:2750537/2‑100 (MQ=255)                                                                                      gccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACCGAAATGAGAAAGCt          >  1:1000377/2‑100 (MQ=255)                                                                                      aCCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGATAAAGCt          >  1:4848261/1‑100 (MQ=255)                                                                                      aCCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:2031623/1‑100 (MQ=255)                                                                                      aCCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:168609/1‑100 (MQ=255)                                                                                      aCCAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCt          >  1:1240841/1‑100 (MQ=255)                                                                                       ccAAGATACAGAGAAGTCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:806550/1‑100 (MQ=255)                                                                                       ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGGATCAACTGAAATGAGAAAGCTa         >  1:1606834/1‑100 (MQ=255)                                                                                       ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:4919699/1‑100 (MQ=255)                                                                                       ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:3150708/1‑100 (MQ=255)                                                                                       ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:392088/1‑100 (MQ=255)                                                                                       ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:1968744/1‑100 (MQ=255)                                                                                       ccAAGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTa         >  1:4056579/1‑100 (MQ=255)                                                                                         caGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAg       >  1:4122592/2‑100 (MQ=255)                                                                                          tgATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:5763179/2‑100 (MQ=255)                                                                                          aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:2511944/1‑100 (MQ=255)                                                                                          aGATACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGt      >  1:1196082/1‑100 (MQ=255)                                                                                           taTACAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGtt     >  1:1783069/2‑100 (MQ=255)                                                                                              aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:2124442/1‑100 (MQ=255)                                                                                              aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:3809847/1‑100 (MQ=255)                                                                                              aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:866299/1‑100 (MQ=255)                                                                                              aCAGAGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAg  >  1:101537/1‑100 (MQ=255)                                                                                                  |                                                                                                AATCAGGATAGTAAGTCAGCATATAACAGAGATCATTATTTAACAGACGATGTATCTAAAAAACAAAATTCATTAGATTCAGTGGACCAAGATACAGTGAAATCAAAATATTATGAGCAAAATTCTGAAGCGACTTTATCAACTAAATCAACCGATAAAGTAGAATCAACTGAAATGAGAAAGCTAAGTTCAG  >  NC_007793/1013510‑1013702

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