Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,093,359	C→T	24.4%	P162S (CCG→TCG)	ycdT →	diguanylate cyclase, membrane‑anchored

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,093,359	0	C	T	24.4%	146.4 / 41.4	90	P162S (CCG→TCG)	ycdT	diguanylate cyclase, membrane‑anchored
Reads supporting (aligned to +/- strand):  ref base C (25/42);  new base T (11/11);  total (36/54)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCA  >  NC_000913/1093261‑1093457                                                                                                   |                                                                                                   tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGcc                                                                                                   <  1:2041975/100‑1 (MQ=255) tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGcc                                                                                                   <  1:2042114/100‑1 (MQ=255) tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGcc                                                                                                   <  1:513159/100‑1 (MQ=255) tGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGcc                                                                                                   <  2:1025661/100‑1 (MQ=255)    cATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGAt                                                                                                >  2:1718936/1‑100 (MQ=255)      ttGGTGAATAAAGTTAGTGAAATCACAAAACAGCGCAATTTATTTTCCAAAAAAATTACTGTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAg                                                                                              <  2:2236996/100‑1 (MQ=255)      ttGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAg                                                                                              >  1:348294/1‑100 (MQ=255)       tGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGt                                                                                             <  2:1153879/100‑1 (MQ=255)       tGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGt                                                                                             >  2:2244776/1‑100 (MQ=255)        ggTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTTTTTTTTTTTTGCCGATAGtt                                                                                            >  2:1736844/1‑100 (MQ=255)        ggTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGtt                                                                                            <  2:582438/100‑1 (MQ=255)        ggTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGtt                                                                                            >  1:1067471/1‑100 (MQ=255)           gAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCt                                                                                         <  2:28420/100‑1 (MQ=255)           gAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCt                                                                                         >  1:1021787/1‑100 (MQ=255)           gAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCt                                                                                         >  2:7239/1‑100 (MQ=255)             aTAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGTTGCTca                                                                                       <  2:2186566/100‑1 (MQ=255)             aTAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTT‑GGGGGTCGATAGTTGCt                                                                                         >  2:13133/1‑97 (MQ=255)             aTAAAGTTACAGAAAACACAAAACTGCGCAATATATTTTTCAAAAAAATGACTTTGTGCATAAGTTTTTTTTTTGTTTTTGGGGGGCCGATAGTTGCTca                                                                                       <  2:1892093/100‑1 (MQ=255)              tAAAGTTAGTGAAATCACAAAACAGCGCGATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTcac                                                                                      <  1:611042/100‑1 (MQ=255)              tAAAGTTAGTGAAAACACAAAACAGCGTAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGTGTCCGATAGTTGCTcac                                                                                      >  2:40921/1‑100 (MQ=255)              tAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTcac                                                                                      <  1:1558643/100‑1 (MQ=255)              tAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTcac                                                                                      <  1:1922248/100‑1 (MQ=255)              tAAAGTTAGTGAAAACACAAAACAGCGAAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGTTGCTca                                                                                       <  1:59446/99‑1 (MQ=255)               aaaGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTcaca                                                                                     <  2:703886/100‑1 (MQ=255)                aaGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACAt                                                                                    <  2:871114/100‑1 (MQ=255)                 agttagttAAATCACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACAt                                                                                    <  1:770291/99‑1 (MQ=255)                    tAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACtt                                                                                >  2:1165005/1‑100 (MQ=255)                    tAGAGAAAACACAAAACGGCGCAATTTATTTTGCAGAAAGGTGAGTTTGTGCATAAGTTTGTTTTTTGGTTTTGGGGGGCCGATAGTTGCTc                                                                                        <  2:2059842/92‑1 (MQ=255)                      gTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGTTGCTCACATACTTTc                                                                              >  2:849891/1‑100 (MQ=255)                      gTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTCTGTGCATAAGTTTGTTTTTTGTTTTTGGGAGCCCGATAGTTGCTCACATACTTTc                                                                              >  2:445460/1‑100 (MQ=255)                        gAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCt                                                                             <  2:74767/99‑1 (MQ=255)                         aaaaCACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAg                                                                           <  2:599627/100‑1 (MQ=255)                           aaCACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTc                                                                         <  2:1449624/100‑1 (MQ=255)                             cacaAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGTTGCTCACATACTTTCTAGTCAt                                                                       >  1:1683704/1‑100 (MQ=255)                             cacaAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCAt                                                                       <  1:1592878/100‑1 (MQ=255)                              acaAAACAGCGTAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCCAGTCAtt                                                                      <  1:1173905/100‑1 (MQ=255)                              acaAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCAtt                                                                      <  2:1067469/100‑1 (MQ=255)                               caAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTa                                                                     <  2:604926/100‑1 (MQ=255)                               caAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTa                                                                     <  1:2224902/100‑1 (MQ=255)                                aaaaCAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCAtt                                                                      <  2:907058/98‑1 (MQ=255)                                 aaaCAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATg                                                                   <  1:1251709/100‑1 (MQ=255)                                    cAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATgag                                                                 <  1:124642/99‑1 (MQ=255)                                     aGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGagag                                                               <  2:751729/100‑1 (MQ=255)                                      gcgcAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGc                                                              >  1:419618/1‑100 (MQ=255)                                      gcgcAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGc                                                              >  1:275673/1‑100 (MQ=255)                                      gcgcAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTAATTATGagag                                                               >  1:1826744/1‑99 (MQ=255)                                        gcAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTa                                                            >  1:442017/1‑100 (MQ=255)                                         cAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTAc                                                           >  1:1499333/1‑100 (MQ=255)                                          aaTTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTa                                                            <  1:423314/98‑1 (MQ=255)                                          aaTTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACa                                                          <  2:1880104/100‑1 (MQ=255)                                            tttatttTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATgaga                                                                <  1:161759/92‑1 (MQ=255)                                               atttTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTa                                                     >  1:552103/1‑100 (MQ=255)                                                ttttCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTac                                                    <  2:2222626/100‑1 (MQ=255)                                                ttttCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTac                                                    <  1:1408088/100‑1 (MQ=255)                                                 tttCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTTATTATGAGAGCTACAATTTac                                                    <  2:1138107/99‑1 (MQ=255)                                                  ttCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTacac                                                  >  1:202933/1‑100 (MQ=255)                                                   ctcAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTAcaca                                                 <  1:1153773/98‑1 (MQ=255)                                                     cAATAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACAtt                                               >  2:2001215/1‑100 (MQ=255)                                                     cAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTTTTGTGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACAtt                                               >  2:624930/1‑100 (MQ=255)                                                     cAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTAcacaag                                               >  2:1153769/1‑98 (MQ=255)                                                     cAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACAtt                                               >  1:595293/1‑100 (MQ=255)                                                     cAAAAAAATGACATTGAGCATAAGATTGTTTTTTGTTTTTGAGGGATCAATAATTACTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACAtt                                               >  2:785297/1‑100 (MQ=255)                                                         aaaaTGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCg                                           <  1:918485/100‑1 (MQ=255)                                                         aaaaTGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCg                                           <  2:1001724/100‑1 (MQ=255)                                                         aaaaTGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCg                                           >  1:1938421/1‑100 (MQ=255)                                                         aaaaTGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCg                                           <  1:496260/100‑1 (MQ=255)                                                           aaTGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAg                                         <  2:1146585/100‑1 (MQ=255)                                                                cTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAAc                                    <  1:315849/100‑1 (MQ=255)                                                                 tttGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAAcc                                   <  1:2129906/100‑1 (MQ=255)                                                                  ttttgCATAATTTTTTTTTTTGTTTTTGGTGTGCCGATATTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAAcc                                   <  2:2229054/96‑1 (MQ=255)                                                                  ttGTGCATAAGTTTGTTTTTTGTTATTGGGTGGCAGATAGTTGCTCACATAATTTCTAGTCATTATTAGAGCTACAATTTACACATTTCCGAGTTAACCa                                  >  2:1203650/1‑100 (MQ=255)                                                                   tgtgCATAAGTTTGTTATATGTAATAGGGGGGTCGATAGTTACTCACATACTTTCTAGTAATTATGAGAGTTACAATTTACACATTGCAGAGTTAACCaa                                 >  2:304431/1‑100 (MQ=255)                                                                    gtgCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAAt                                <  1:1983646/100‑1 (MQ=255)                                                                     tgCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATg                               <  1:18616/100‑1 (MQ=255)                                                                      gCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGa                              >  1:1124835/1‑100 (MQ=255)                                                                         tAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGaaa                            >  1:294675/1‑99 (MQ=255)                                                                           aGTTTGTTTTTTGTTTTTGGGGGGTCAATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATg                         >  1:1366397/1‑100 (MQ=255)                                                                           aGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATg                         <  2:419619/100‑1 (MQ=255)                                                                           aGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATg                         <  2:1826732/100‑1 (MQ=255)                                                                            gtttgtttTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATgg                        <  1:1432311/100‑1 (MQ=255)                                                                               tgttttttgtttttGGGGGGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCa                     <  2:1716651/100‑1 (MQ=255)                                                                               tgttttttgtttttGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTc                      >  1:1151966/1‑99 (MQ=255)                                                                                 ttttttgtttttGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAgg                   >  1:593533/1‑100 (MQ=255)                                                                                 ttttttgtttttGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAgg                   <  1:1912789/100‑1 (MQ=255)                                                                                  tttttgtttttgGGGGGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAggt                  <  1:287136/100‑1 (MQ=255)                                                                                   ttttgtttttgGGGGGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAggtg                 <  2:927547/100‑1 (MQ=255)                                                                                    tttgtttttgGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAggtgg                <  1:1135359/100‑1 (MQ=255)                                                                                       gtttttgGGGGGTCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCTAGTTAACCAATGAAAAAGGTCAGGTGGttt             >  2:1935554/1‑100 (MQ=255)                                                                                        tttttgGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTg            <  1:878030/100‑1 (MQ=255)                                                                                           ttgGGTGGACGATAGTTGCTCACATACTTTCTAGTTATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGaa         <  2:63762/100‑1 (MQ=255)                                                                                             ggggggCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGaaa        >  1:483384/1‑99 (MQ=255)                                                                                                  gCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTc   <  2:2041091/99‑1 (MQ=255)                                                                                                  gCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCa  >  1:271152/1‑100 (MQ=255)                                                                                                  gCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCa  <  1:1166430/100‑1 (MQ=255)                                                                                                   |                                                                                                   TGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCA  >  NC_000913/1093261‑1093457

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