Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,652,855	(G)6→7	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	98.9%	312.3 / ‑3.4	88	coding (356/846 nt)	sseA	3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (36/51);  total (36/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943                                                                                                      |                                                                                                aGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCt                                                                                                  <  1:723595/100‑1 (MQ=255)     tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg                                                                                             <  1:585542/100‑1 (MQ=255)     tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggg                                                                                             >  2:445266/1‑100 (MQ=255)        gCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑ggggggg                                                                                          <  1:1072927/100‑2 (MQ=255)         cACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggggt                                                                                         >  1:303425/1‑98 (MQ=255)          aCCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggggtg                                                                                        <  2:1416368/100‑4 (MQ=255)           ccggATTGTCTATGACGGAGGTAATCTTTTCTCAGCCCGACGAGCATGGTGGGTGCGGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggggtgg                                                                                       <  2:587421/97‑5 (MQ=255)           ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggggtgg                                                                                       >  1:2285341/1‑96 (MQ=255)           ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑gggggggtg                                                                                        >  1:2700943/1‑96 (MQ=255)             ngATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGAc                                                                                     >  2:2265298/2‑100 (MQ=255)              gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGTTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACt                                                                                    <  1:1207629/100‑1 (MQ=255)              gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACt                                                                                    <  2:45294/100‑1 (MQ=255)                 tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CGCCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGc                                                                                 <  2:119376/100‑1 (MQ=255)                 tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGc                                                                                 <  2:32827/100‑1 (MQ=255)                  gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCa                                                                                >  2:516027/1‑100 (MQ=255)                    cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTCGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAgg                                                                              >  2:279626/1‑100 (MQ=255)                      aTGCCGAAGGTAATCGTTCCTCCGAGACAGGCGCTTGGTGGATGCTGCGCCAGCTCT‑TTGGAGGGAAAGTGGATAGTCT‑GGGGGGTTGGACTTGCAGGCt                                                                            <  2:1270320/100‑1 (MQ=255)                      aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTCGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCt                                                                            <  2:1779028/100‑1 (MQ=255)                      aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGc                                                                             >  2:1181364/1‑99 (MQ=255)                      aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCt                                                                            <  1:2468153/100‑1 (MQ=255)                       tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTg                                                                           <  1:818005/100‑1 (MQ=255)                        gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTgg                                                                          <  2:1113192/100‑1 (MQ=255)                         aCGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGc                                                                         <  2:1847786/100‑1 (MQ=255)                              ggTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAgcgc                                                                    >  1:370126/1‑100 (MQ=255)                               gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGcgcg                                                                   <  2:1130943/100‑1 (MQ=255)                               gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGcgcg                                                                   <  1:1104715/100‑1 (MQ=255)                                tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCga                                                                  >  2:384569/1‑100 (MQ=255)                                 aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                 >  2:1266584/1‑100 (MQ=255)                                 aaTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgat                                                                 >  2:1823996/1‑100 (MQ=255)                                   tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  2:2781507/1‑100 (MQ=255)                                   tCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatga                                                               >  2:2002418/1‑100 (MQ=255)                                    cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCgatgat                                                              >  2:2058973/1‑100 (MQ=255)                                     ttttCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATc                                                             >  1:1917318/1‑100 (MQ=255)                                      tttCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATct                                                            <  1:445267/100‑1 (MQ=255)                                      tttCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATct                                                            <  1:1181368/100‑1 (MQ=255)                                       ttCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctg                                                           <  2:415347/100‑1 (MQ=255)                                       ttCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctg                                                           <  1:1452447/100‑1 (MQ=255)                                       ttCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTGTGGTGTAGAGAAAGTGGCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctg                                                           <  2:415361/100‑1 (MQ=255)                                        tctcAGCCCCACGGGCATGGTGGATGCGGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          <  2:2417808/100‑1 (MQ=255)                                        tctcAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgc                                                          <  1:1868195/100‑1 (MQ=255)                                         ctcAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                         >  2:769537/1‑100 (MQ=255)                                          tcAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                        <  1:1623341/100‑1 (MQ=255)                                             gCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  2:540791/100‑1 (MQ=255)                                             gCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  2:381744/100‑1 (MQ=255)                                              ccccACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                    <  2:1012137/100‑1 (MQ=255)                                              ccccACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                    <  2:1994150/100‑1 (MQ=255)                                               cccACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                   <  1:2606465/100‑1 (MQ=255)                                                ccACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                  >  2:1493528/1‑100 (MQ=255)                                                ccACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                    >  1:342919/1‑98 (MQ=255)                                                 cACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 >  1:2201808/1‑100 (MQ=255)                                                 cACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                 >  2:2218233/1‑100 (MQ=255)                                                  aCGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaaga                                                <  1:2781528/100‑1 (MQ=255)                                                  aCGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                     <  2:1596031/95‑1 (MQ=255)                                                     aGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:935928/100‑1 (MQ=255)                                                     aGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                             <  2:935823/100‑1 (MQ=255)                                                      gCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                            >  2:1832083/1‑100 (MQ=255)                                                      gCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                            >  2:1721321/1‑100 (MQ=255)                                                      gCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                            >  1:830520/1‑100 (MQ=255)                                                        aTGGTGNATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          <  1:2218246/100‑1 (MQ=255)                                                        aTGGTGGATGCTGCG‑CACCTTTGGTGTAGATAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTgn                                          <  1:2265311/100‑2 (MQ=255)                                                        aTGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          <  1:591794/100‑1 (MQ=255)                                                        aTGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                          <  1:1643772/100‑1 (MQ=255)                                                         tggtggATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCa                                         >  1:1271259/1‑100 (MQ=255)                                                           gtggATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  2:459248/100‑1 (MQ=255)                                                           gtggATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                       <  2:510739/100‑1 (MQ=255)                                                               aTGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                   <  2:2700923/100‑1 (MQ=255)                                                                    gcg‑cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              >  2:572162/1‑100 (MQ=255)                                                                    gcg‑cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              <  1:2058987/100‑1 (MQ=255)                                                                    gcg‑cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGcc                              >  2:435467/1‑100 (MQ=255)                                                                     cg‑cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                             <  1:2613214/100‑1 (MQ=255)                                                                     cg‑cACCTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                             <  2:1999150/100‑1 (MQ=255)                                                                      g‑cACGTTTGGTGTAGGGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGGAGTAGAGCTGCCgg                            <  1:147243/100‑1 (MQ=255)                                                                          ccTTTGATGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                         <  2:2267987/100‑1 (MQ=255)                                                                           cTTTGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                        <  1:101959/100‑1 (MQ=255)                                                                             ttGGTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgag                      <  1:516029/100‑1 (MQ=255)                                                                               ggTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    <  1:103724/100‑1 (MQ=255)                                                                               ggTGTAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                    <  2:2524848/100‑1 (MQ=255)                                                                                 tgtAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                  <  2:1355786/100‑1 (MQ=255)                                                                                  gtAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 <  2:342918/100‑1 (MQ=255)                                                                                  gtAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                 >  2:910481/1‑100 (MQ=255)                                                                                   tAGAGAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa                <  2:1271254/100‑1 (MQ=255)                                                                                    agagAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  2:764550/1‑100 (MQ=255)                                                                                    agagAAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa               >  1:428206/1‑100 (MQ=255)                                                                                     gagaAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGGGCTGCCGGAAGGAGAGTTTAAc              >  2:2216951/1‑100 (MQ=255)                                                                                     gagaAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAc              <  2:418242/100‑1 (MQ=255)                                                                                       gaAAGTGTCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcgc            >  1:479531/1‑100 (MQ=255)                                                                                            tgtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGt       <  1:2676597/100‑1 (MQ=255)                                                                                             gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:828683/1‑100 (MQ=255)                                                                                             gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:552501/1‑100 (MQ=255)                                                                                             gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      <  2:455280/100‑1 (MQ=255)                                                                                             gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  1:2315388/1‑100 (MQ=255)                                                                                             gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:2697085/1‑100 (MQ=255)                                                                                             gtCGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGtt      >  2:2869272/1‑100 (MQ=255)                                                                                               cGATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTa    >  1:923563/1‑100 (MQ=255)                                                                                                gATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa   >  1:1111775/1‑100 (MQ=255)                                                                                                gATTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGAGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa   >  1:843916/1‑100 (MQ=255)                                                                                                 aTTCTGGGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAt  <  2:2577712/100‑1 (MQ=255)                                                                                                      |                                                                                                AGGATAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCG‑CACCTTTGGTGTAGAGAAAGTGTCGATTCT‑GGGGGG‑TGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAAT  >  NC_000913/2652750‑2652943

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