D3A :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,101,334	A→G	S75S (AGT→AGC)	SAUSA300_RS05410 ←	divalent metal cation transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,101,334	0	A	G	85.6%	1035.2 / 124.9	547	S75S (AGT→AGC)	SAUSA300_RS05410	divalent metal cation transporter
Reads supporting (aligned to +/- strand): ref base A (22/57); new base G (3/465); total (25/522)
Fisher's exact test for biased strand distribution p-value = 3.26e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-04

Only 100 of 550 total aligned reads displayed.

TAGGTCTTGATAAATAATGTCGTGTCATTTGAGCTAAGTCCATACCTGTTGCTATTCCTAATCTCACTGTCATGCTTTGAAGTAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATCCAGTTTCCGGGATCCATGTAACCAACAGCGACT  >  NC_007793/1101237‑1101433

caggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCt                                                                                                   <  1:3833604/90‑1 (MQ=255)

caggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCt                                                                                                   <  1:4369157/90‑1 (MQ=255)

caggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCt                                                                                                   <  1:4667878/90‑1 (MQ=255)

caggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCt                                                                                                   <  1:4248730/90‑1 (MQ=255)

caggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCt                                                                                                   <  1:1660205/90‑1 (MQ=255)

caggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCt                                                                                                   <  1:4923201/90‑1 (MQ=255)

caggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCt                                                                                                   <  1:1227488/90‑1 (MQ=255)

caggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCt                                                                                                   <  1:1183446/90‑1 (MQ=255)

caggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCt                                                                                                   <  1:5226366/90‑1 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCTGTTGCTATTCCTAATCTTACAGTCATACTGTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:6177315/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:314333/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:5895614/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:5431341/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:3957797/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:5655233/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:3622381/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:2452074/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:3599086/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:4573560/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:5082801/92‑3 (MQ=255)

  ggcttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtag                                                                                                 <  1:4794963/92‑3 (MQ=255)

    cttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtaggg                                                                                               <  1:5559185/94‑5 (MQ=255)

    cttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGag                                                                                               <  1:604939/94‑2 (MQ=255)

    cttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGaa                                                                                               <  1:156643/94‑1 (MQ=255)

    cttatttaaaaaaTGACGTGTCATTTGTGCTAAATCCATAACAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCtaggg                                                                                               <  1:6108205/94‑5 (MQ=255)

        tttaaaaaaTGACGTGTCATTTGTGCTAAGTCCATACCAGTTGCTATTCCTAATCTTCCAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGAtta                                                                                           <  1:5598134/98‑3 (MQ=255)

        tttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGAtta                                                                                           <  1:1716980/98‑3 (MQ=255)

        tttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGAtta                                                                                           <  1:4048606/98‑3 (MQ=255)

        tttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGAtta                                                                                           <  1:2568799/98‑3 (MQ=255)

        tttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGAtta                                                                                           <  1:224495/98‑3 (MQ=255)

        tttaaaaaaTGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACGTTGTAACAGCATAGCAGATAAGCTAGAGAtta                                                                                           <  1:6532285/98‑3 (MQ=255)

                atGACGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACg                                                                                   <  1:4838100/100‑1 (MQ=255)

                  gacGTGTCATTTGTGCTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGaa                                                                                 >  1:1850852/3‑100 (MQ=255)

                                 ttAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:6370366/99‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATTCTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:4458517/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGATATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:5180429/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:3875748/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:1019785/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:3437372/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:1627972/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:1884514/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:2039743/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:3019337/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:642346/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:130597/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:2750018/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTGTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:6321818/100‑1 (MQ=255)

                                 cTAAATCCATACCAGTTGCTATTCCTAATCTTACAGACATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGcc                                                                  <  1:5440589/100‑1 (MQ=255)

                                   aaatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCat                                                                <  1:4117564/97‑1 (MQ=255)

                                   aaatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCat                                                                <  1:5486109/97‑1 (MQ=255)

                                   aaatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCat                                                                <  1:6500136/97‑1 (MQ=255)

                                   aaatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCat                                                                <  1:824628/97‑1 (MQ=255)

                                   aaatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCat                                                                <  1:2832418/97‑1 (MQ=255)

                                   aaatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCat                                                                <  1:5252029/97‑1 (MQ=255)

                                   aaatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCat                                                                <  1:321531/97‑1 (MQ=255)

                                   aaatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCat                                                                <  1:327894/97‑1 (MQ=255)

                                   aaatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCat                                                                <  1:6091393/97‑1 (MQ=255)

                                    aatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCata                                                               <  1:5016900/98‑1 (MQ=255)

                                    aatCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCata                                                               <  1:5946143/98‑1 (MQ=255)

                                    aGTCCATACCTGTTGCTATTCCTAATCTCACTGTCATGCTTTGAAGTAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCatg                                                               <  1:570839/100‑2 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:2144758/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:900215/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:5389496/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:1760917/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:339726/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:2216854/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:2490934/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:5852156/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:297136/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:3059950/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatat                                                              <  1:3166104/99‑1 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatag                                                              <  1:700843/99‑2 (MQ=255)

                                     atCCATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCatag                                                              <  1:1049190/99‑2 (MQ=255)

                                       ccATACCTGTTGCTATTCCTAATCTCACTGTCATGCTTTGAAGTAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTg                                                            >  1:5788037/1‑100 (MQ=255)

                                       ccATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTg                                                            <  1:4692843/100‑1 (MQ=255)

                                       ccATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTg                                                            <  1:6583611/100‑1 (MQ=255)

                                       ccATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTg                                                            <  1:1591951/100‑1 (MQ=255)

                                       ccATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTg                                                            <  1:6136139/100‑1 (MQ=255)

                                       ccATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTg                                                            <  1:3363060/100‑1 (MQ=255)

                                       ccATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTg                                                            <  1:5997339/100‑1 (MQ=255)

                                       ccATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTg                                                            <  1:3747221/100‑1 (MQ=255)

                                       ccATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTg                                                            <  1:4315232/100‑1 (MQ=255)

                                       ccATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTg                                                            <  1:4629306/100‑1 (MQ=255)

                                        cATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTgg                                                           <  1:4161585/100‑1 (MQ=255)

                                        cATACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTgc                                                           <  1:196963/100‑2 (MQ=255)

                                         aTACCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTgct                                                          <  1:6253137/100‑3 (MQ=255)

                                           aCCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTGTgct                                                        <  1:2915950/100‑2 (MQ=255)

                                           aCCAGTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTGTgcc                                                        <  1:5310511/100‑2 (MQ=255)

                                             cagTTGCTATTCCTAATCTTACAGTCATACTTTGTAACAGCATAGCAGATAAGCTAGAGATTAAGATTACGAATAACAAGGTATAGCCATATTgggggag                                                      <  1:2102358/98‑5 (MQ=255)

                                                  tgctcttccgaTCTCACTGTCATGCTTTGAAGTAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCtt                                                 <  1:4091889/90‑1 (MQ=255)

                                                               tcACTGTCATGCTTTGAAGCAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATCGCCATATTGGGCGCCACCTTGCATTGATGTTAt                                    >  1:3491147/1‑100 (MQ=255)

                                                                cACTGTCATGCTTTGAAGTAACATTGCTGGTAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATc                                   >  1:4219398/1‑100 (MQ=255)

                                                                cACTGTCATGCTTTGAAGTAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATc                                   >  1:3774498/1‑100 (MQ=255)

                                                                 aCTGTCATGCTTTGAAGTAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATcc                                  <  1:2806740/100‑1 (MQ=255)

                                                                 aCTGTCATGCTTTGAAGTAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATcc                                  <  1:1461377/100‑1 (MQ=255)

                                                                         gCTTTGAAGTAACATTGCTGGTAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATCCAGTTTCCg                          <  1:5755486/100‑1 (MQ=255)

                                                                             tGAAGTAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATCCAGTTTCCGGGaa                      <  1:2626748/100‑2 (MQ=255)

                                                                                          cTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATCCAGTTTCCGGGATCCATGTAACCagg         <  1:1073911/100‑3 (MQ=255)

                                                                                                 aCTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATCCAGTTTCCGGGATCCATGTAACCAACAGCGagt  <  1:2357040/100‑3 (MQ=255)

                                                                                                 aCTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATCCAGTTTCCGGGATCCATGCAACCAACAGCGAcg  <  1:112340/100‑2 (MQ=255)

TAGGTCTTGATAAATAATGTCGTGTCATTTGAGCTAAGTCCATACCTGTTGCTATTCCTAATCTCACTGTCATGCTTTGAAGTAACATTGCTGATAAACTTGAAATAAGAATTACGAATAGCAAAGTATAGCCATATTGGGCGCCACCTTGCATTGATGTTATCCAGTTTCCGGGATCCATGTAACCAACAGCGACT  >  NC_007793/1101237‑1101433

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