Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,946,790	T→G	100%	G180G (GGT→GGG)	yebB →	DUF830 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,946,790	0	T	G	96.3%	262.7 / 0.1	81	G180G (GGT→GGG)	yebB	DUF830 family protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (39/39);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

CGAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAG  >  NC_000913/1946693‑1946884                                                                                                  |                                                                                               cGAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGgtt                                                                                              >  1:846200/1‑100 (MQ=255) cGAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGgtt                                                                                              <  1:1722528/100‑1 (MQ=255) cGAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGgtt                                                                                              <  1:1722537/100‑1 (MQ=255)   aaTGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTgg                                                                                            <  2:2638969/100‑1 (MQ=255)    aTGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGgt                                                                                           <  2:2438898/100‑1 (MQ=255)     tGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGgtg                                                                                          <  1:2805862/100‑1 (MQ=255)     tGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGgtg                                                                                          <  1:2946661/100‑1 (MQ=255)      gCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGgtgt                                                                                         <  2:2727492/100‑1 (MQ=255)      gCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGCTGgtg                                                                                          >  2:916935/1‑99 (MQ=255)      gCAAAACTCACTTTCTGGAAATTCTGGTTCCTAGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGgtgt                                                                                         <  1:732245/100‑1 (MQ=255)       cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGtt                                                                                        >  2:2971190/1‑100 (MQ=255)       cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCAAGCCAGTTTGTGGCATCATCCGGGGTTGgtgt                                                                                         >  1:1335777/1‑99 (MQ=255)        aaaaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTg                                                                                       >  2:2687510/1‑100 (MQ=255)          agcTCACTTTCTGGAAATTCTGGTTCCTGGGTTATATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGAt                                                                                     <  1:3046648/98‑1 (MQ=255)           aCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGAtt                                                                                    >  2:1933310/1‑100 (MQ=255)             tCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCa                                                                                  <  2:430576/100‑1 (MQ=255)                     tGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggagg                                                                          >  2:2779380/1‑100 (MQ=255)                     tGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggag                                                                           >  1:733440/1‑99 (MQ=255)                      ggAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCATGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggaggg                                                                         <  2:1554320/100‑1 (MQ=255)                      ggAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggaggg                                                                         >  2:645815/1‑100 (MQ=255)                      ggAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggaggg                                                                         >  1:857920/1‑100 (MQ=255)                      ggAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATAATCCGGGGTTGGTGTTGATTCATGCggaggg                                                                         >  1:1134796/1‑100 (MQ=255)                           ttCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTgg                                                                    <  2:788218/100‑1 (MQ=255)                            tCTGGTTTCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGa                                                                   <  2:1624421/100‑1 (MQ=255)                               ggTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAAc                                                                >  1:35946/1‑100 (MQ=255)                                gTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACg                                                               <  2:757448/100‑1 (MQ=255)                                 ttCCTGGGTTCTATTCCGTGGGAGCGTAAAACCTTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGc                                                              <  1:2357057/100‑1 (MQ=255)                                  tCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGcc                                                             >  1:2577445/1‑100 (MQ=255)                                  tCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGcc                                                             <  1:2971223/100‑1 (MQ=255)                                  tCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGcc                                                             <  1:1933333/100‑1 (MQ=255)                                   ccTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCAAAGGTTGATGTTGATTCATGCGGAGGAAGTGGAAACGCCt                                                            >  1:72666/1‑100 (MQ=255)                                    cTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTc                                                           >  2:1828163/1‑100 (MQ=255)                                    cTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTc                                                           >  2:2962263/1‑100 (MQ=255)                                      gggTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAg                                                         >  1:2360908/1‑100 (MQ=255)                                       ggTTCTATTCCGTGGGAGCGTAAAATCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGc                                                        <  1:2827705/100‑1 (MQ=255)                                        gTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGcc                                                       <  1:850853/100‑1 (MQ=255)                                           cTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGa                                                    >  2:2186138/1‑100 (MQ=255)                                            tATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGaa                                                   >  1:429124/1‑100 (MQ=255)                                            tATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGaa                                                   <  2:2410325/100‑1 (MQ=255)                                            tATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGa                                                    <  1:2927051/99‑1 (MQ=255)                                             aTTCCGTGGGAGCGTGAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                  >  1:1155511/1‑100 (MQ=255)                                             aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                  >  2:2610346/1‑100 (MQ=255)                                             aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                  <  2:2421059/100‑1 (MQ=255)                                             aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                  <  2:1696787/100‑1 (MQ=255)                                              ttCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACt                                                 <  2:2862548/100‑1 (MQ=255)                                               tCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTg                                                >  2:816950/1‑100 (MQ=255)                                               tCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTg                                                >  2:985039/1‑100 (MQ=255)                                               tCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTg                                                <  1:1729772/100‑1 (MQ=255)                                                 cGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGAc                                              <  1:995977/100‑1 (MQ=255)                                                 cGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGAc                                              >  1:2485855/1‑100 (MQ=255)                                                          gTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGt                                     >  1:1405481/1‑100 (MQ=255)                                                           tAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGta                                    <  1:2743613/100‑1 (MQ=255)                                                           tAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGta                                    <  2:2140377/100‑1 (MQ=255)                                                           tAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGta                                    >  1:2486730/1‑100 (MQ=255)                                                           tAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACAGAGGCGGta                                    >  1:1364176/1‑100 (MQ=255)                                                             aaaCCGTCACGCCAGCCAGTTTGTGGCATCAACCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCGCAGCCGGAACTGACCGAGGCGGtata                                  >  2:333022/1‑100 (MQ=255)                                                               aCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAAc                                <  2:1657959/100‑1 (MQ=255)                                                                 cGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACtt                              >  1:393248/1‑100 (MQ=255)                                                                  gttaCGCCAGCCAGTTTGTGGCATCATCCTGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTa                             <  1:949223/97‑1 (MQ=255)                                                                  gTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTa                             >  2:1569725/1‑100 (MQ=255)                                                                  gTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTa                             >  2:2803163/1‑100 (MQ=255)                                                                  gTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGAGGAAACACCTCAGCCGGAACTGACCGAGGCGGTATAACtaa                             >  1:420705/1‑98 (MQ=255)                                                                     aCGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACg                          >  2:1180407/1‑100 (MQ=255)                                                                         cagccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAg                       <  1:731299/99‑1 (MQ=255)                                                                         cagccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGt                      <  1:2530460/100‑1 (MQ=255)                                                                             cagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGcc                  <  1:1562444/100‑1 (MQ=255)                                                                              agTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGATACGCCTCAGCCGGAACTGACCGAGGCGGCATAACTTAACGCAGTCGcct                 <  2:1848030/100‑2 (MQ=255)                                                                                 ttGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCctct              <  2:857916/100‑1 (MQ=255)                                                                                   gtgGCATCATCCGGGGTTGGTGTTGATTCATGCGGTGGGAGTGGAAACATCTCAGTCGGAACTGACCGAGTAGTTATAACTTAACGCAATCGCCCTCTCg            >  2:2230089/1‑100 (MQ=255)                                                                                   gtgGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg            >  1:3036047/1‑100 (MQ=255)                                                                                   gtgGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg            <  2:1699002/100‑1 (MQ=255)                                                                                   gtgGCATCATCCGGGGATGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTAAGCCGGAACTGACCGAGGCAGTATAACTTAACGCAGTCGCCCTCACg            >  2:1205363/1‑100 (MQ=255)                                                                                    tgGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGc           <  2:1049121/100‑1 (MQ=255)                                                                                      gCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCa         >  2:853005/1‑100 (MQ=255)                                                                                      gCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCa         <  2:1130391/100‑1 (MQ=255)                                                                                      gCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGAGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCa         >  2:1866714/1‑100 (MQ=255)                                                                                       catcatCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAg        <  2:1203719/100‑1 (MQ=255)                                                                                        atcatcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAgg       >  2:2433631/1‑100 (MQ=255)                                                                                        atcatcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAgg       <  1:2908370/100‑1 (MQ=255)                                                                                        atcatcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAgg       <  1:1896104/100‑1 (MQ=255)                                                                                           atcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc    >  1:1252329/1‑100 (MQ=255)                                                                                            tcCGGGGTTCGTTTTGATTCATGCGAAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCa   <  2:1364164/100‑1 (MQ=255)                                                                                             ccGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAg  <  1:767144/100‑1 (MQ=255)                                                                                             ccGGGGTTGGTGTTGATTCATGAGGAGGGAGTGTAAACGCCTCAGCCTGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCTCCAGGt      >  2:549513/1‑96 (MQ=255)                                                                                               ggggTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACAGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAg  >  2:1390554/1‑98 (MQ=255)                                                                                               ggggTTGGTGTTGATTCATGAGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCa   <  2:2486701/97‑1 (MQ=255)                                                                                                  |                                                                                               CGAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAG  >  NC_000913/1946693‑1946884

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: ‑