Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,652,855:1	(G)6→7	100%	coding (362/846 nt)	sseA →	3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,652,849	1	.	G	97.5%	292.0 / 0.5	81	coding (356/846 nt)	sseA	3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (34/45);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652750‑2652942                                                                                                     |                                                                                               aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCt                                                                                                 >  2:2282340/1‑100 (MQ=255)     tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg                                                                                           >  1:246944/1‑100 (MQ=255)     tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gg                                                                                           >  1:1224589/1‑100 (MQ=255)      aaccACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑ggg                                                                                          <  2:1947742/97‑1 (MQ=255)      aaGCACCTGATTGGCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑ggg                                                                                          >  2:1459145/1‑100 (MQ=255)       aGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gggg                                                                                         >  1:1699962/1‑100 (MQ=255)        gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑ggggg                                                                                        <  2:151495/100‑2 (MQ=255)        gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑ggggg                                                                                        >  2:38052/1‑99 (MQ=255)        gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑ggggg                                                                                        >  1:66798/1‑99 (MQ=255)        gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑ggggg                                                                                        >  1:110816/1‑99 (MQ=255)        gCACCTGATTGTCTATGACGAAGGTAATATTTTCTCAGCCCCATCGGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGGCGGTTAT‑gg‑ggg                                                                                          >  2:2170310/1‑98 (MQ=255)          aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gggggtg                                                                                      <  1:663204/100‑4 (MQ=255)           ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gggggtgg                                                                                     <  1:310256/100‑5 (MQ=255)           ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gggggtgg                                                                                     <  1:310744/100‑5 (MQ=255)           ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gggggtgg                                                                                     <  2:2171671/100‑5 (MQ=255)           ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gg‑gggggtgg                                                                                     >  1:1121935/1‑96 (MQ=255)                ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTg                                                                                 >  1:18804/1‑100 (MQ=255)                 tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGGGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGc                                                                                <  2:1620710/100‑1 (MQ=255)                  gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCa                                                                               >  2:1795227/1‑100 (MQ=255)                  gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCa                                                                               <  1:2108555/100‑1 (MQ=255)                  gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCa                                                                               <  2:1177676/100‑1 (MQ=255)                     tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGc                                                                            <  2:705072/100‑1 (MQ=255)                         aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGc                                                                        <  1:1692030/100‑1 (MQ=255)                           gAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAg                                                                      <  1:1860195/100‑1 (MQ=255)                           gAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAg                                                                      <  1:1618886/100‑1 (MQ=255)                              ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGGGGACTTGCAGGCTGGCAg                                                                      >  2:2336453/1‑97 (MQ=255)                               gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGcgcg                                                                  <  2:429886/100‑1 (MQ=255)                                tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCga                                                                 >  2:752216/1‑100 (MQ=255)                                 aaTCTTTTCTCAGCCACACGAGCATGGTGGAGGCTGCGCACCTTTGGTGTAGTGGAAGAGTCGAGTAT‑GGTGGGGTGGACTTGCAGGCTgg                                                                         >  2:68426/1‑91 (MQ=255)                                  aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCgatg                                                               <  1:152555/100‑1 (MQ=255)                                    cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCgatga                                                              >  2:85738/1‑99 (MQ=255)                                     ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCgatgatn                                                            <  2:1155978/100‑2 (MQ=255)                                     ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                            >  2:584889/1‑100 (MQ=255)                                     ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATc                                                            <  1:2267641/100‑1 (MQ=255)                                         ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                       <  2:661613/100‑1 (MQ=255)                                          tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                       <  2:2331331/100‑1 (MQ=255)                                            aGCCCCACGAGCATGGTGGATGCTGCGCCCCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     <  1:872399/100‑1 (MQ=255)                                            aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     >  1:2146664/1‑100 (MQ=255)                                            aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     >  1:2141218/1‑100 (MQ=255)                                             gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  2:1200117/100‑1 (MQ=255)                                             gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  1:868921/100‑1 (MQ=255)                                               cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                  >  2:1931116/1‑100 (MQ=255)                                               cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                  <  1:1175685/100‑1 (MQ=255)                                               cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                  >  1:1641036/1‑100 (MQ=255)                                                ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                 >  2:165975/1‑100 (MQ=255)                                                 cACGAGCATGGTGGATGCTGCGCACCTGTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                <  2:119195/100‑1 (MQ=255)                                                   cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                              >  2:1317134/1‑100 (MQ=255)                                                    gAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGaagaag                                             <  1:333626/100‑1 (MQ=255)                                                     aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:2207218/100‑1 (MQ=255)                                                     aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:414732/100‑1 (MQ=255)                                                     aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:416701/100‑1 (MQ=255)                                                         tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                        >  1:1116448/1‑100 (MQ=255)                                                         tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                        >  1:2302941/1‑100 (MQ=255)                                                          ggtggATGCTGTGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAAAAGAAGGTGCAg                                       <  2:31817/100‑1 (MQ=255)                                                          ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                       >  2:32964/1‑100 (MQ=255)                                                           gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                      <  2:1593599/100‑1 (MQ=255)                                                            tggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTa                                     <  2:1540185/100‑1 (MQ=255)                                                             ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGGTGATCTGCTGTTAGAAGAAGGTGCAGTag                                    <  2:2176166/100‑1 (MQ=255)                                                              gATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTaga                                   <  1:1764297/100‑1 (MQ=255)                                                               aTGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                  >  1:2220707/1‑100 (MQ=255)                                                                   tgcGCACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                              >  2:641185/1‑100 (MQ=255)                                                                    gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                             <  1:1813973/100‑1 (MQ=255)                                                                    gcgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGc                              <  2:363866/99‑1 (MQ=255)                                                                     cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            >  1:1710250/1‑100 (MQ=255)                                                                     cgcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            <  1:1062286/100‑1 (MQ=255)                                                                      gcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                           <  1:1842536/100‑1 (MQ=255)                                                                      gcACCTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                           <  2:66798/100‑1 (MQ=255)                                                                         ccTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                        >  2:1051366/1‑100 (MQ=255)                                                                          cTTTGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                       <  2:894224/100‑1 (MQ=255)                                                                            ttGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                     <  1:1156355/100‑1 (MQ=255)                                                                            ttGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                     <  1:2103059/100‑1 (MQ=255)                                                                             tGGTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                     <  1:2031714/99‑1 (MQ=255)                                                                              ggTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   <  2:882302/100‑1 (MQ=255)                                                                              ggTGTAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   >  2:1048238/1‑100 (MQ=255)                                                                               gtgtAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                  <  2:206199/100‑1 (MQ=255)                                                                                tgtAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                 >  1:524114/1‑100 (MQ=255)                                                                                 gtAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                <  2:2169846/100‑1 (MQ=255)                                                                                  tAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa               >  1:1762029/1‑100 (MQ=255)                                                                                  tAGAGAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa               <  1:67345/100‑1 (MQ=255)                                                                                   agagAAAGTGTCGATTCTGGG‑GGGGTGGACTTGGAGGCTGGCATCGCGATGATCTGCTGGTAGGAGAAGGTGGAGTAGAGCTGCCGGAAGGAGAGgttca              >  2:1029861/1‑98 (MQ=255)                                                                                   agagAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGCGTTTaa              >  2:1723528/1‑100 (MQ=255)                                                                                   agagAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              >  2:1419656/1‑100 (MQ=255)                                                                                   agagAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              >  1:1291296/1‑100 (MQ=255)                                                                                   agagAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              >  2:170627/1‑100 (MQ=255)                                                                                   agagAAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              >  1:2196167/1‑100 (MQ=255)                                                                                    gagaAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc             >  1:1285344/1‑100 (MQ=255)                                                                                    gagaAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGAGCGATGATCTGCTGTTATAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              >  2:353454/1‑99 (MQ=255)                                                                                    gagaAAGTGGCGATTCTGGG‑GGGGTTGACTTGGAGGCTGGCAGCGCGATGAGCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc             <  2:231492/100‑1 (MQ=255)                                                                                     agaAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg            <  1:2312961/100‑1 (MQ=255)                                                                                     agaAAGTGTCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg            <  1:1409608/100‑1 (MQ=255)                                                                                          gtgtCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg       <  2:2329392/100‑1 (MQ=255)                                                                                           tgtCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt      >  1:1625846/1‑100 (MQ=255)                                                                                            gtCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt     <  1:869421/100‑1 (MQ=255)                                                                                             tCGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt    <  2:83315/100‑1 (MQ=255)                                                                                              cGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAGGCCGCGTTTa   >  1:2342955/1‑100 (MQ=255)                                                                                              cGATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa   >  1:1080056/1‑100 (MQ=255)                                                                                               gATTCTGGG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa  >  2:2278704/1‑100 (MQ=255)                                                                                                     |                                                                                               AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GG‑GGGGTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652750‑2652942

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