Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,370,616	A→G	56.8%	intergenic (‑204/‑72)	yjeH ← / → groS	putative transporter/Cpn10 chaperonin GroES, small subunit of GroESL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,370,616	0	A	G	56.8%	48.2 / 111.1	95	intergenic (‑204/‑72)	yjeH/groS	putative transporter/Cpn10 chaperonin GroES, small subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base G (29/25);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

AGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAT  >  NC_000913/4370520‑4370711                                                                                                 |                                                                                                aGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCa                                                                                              >  2:1670990/1‑100 (MQ=255) aGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCa                                                                                              >  1:1255854/1‑100 (MQ=255)  ggTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAg                                                                                             >  1:15597/1‑100 (MQ=255)  ggTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAg                                                                                             >  1:2069717/1‑100 (MQ=255)   gTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGc                                                                                            <  2:2153379/100‑1 (MQ=255)    tGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGcc                                                                                           >  1:1907401/1‑100 (MQ=255)      aaTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCgg                                                                                         <  1:1193179/100‑1 (MQ=255)       aTCAGGGTTTTCAGCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCggg                                                                                        <  2:435246/100‑1 (MQ=255)       aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCggg                                                                                        <  2:1295703/100‑1 (MQ=255)       aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCggg                                                                                        <  1:759410/100‑1 (MQ=255)       aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCggg                                                                                        >  1:1576284/1‑100 (MQ=255)       aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCggg                                                                                        <  1:1698739/100‑1 (MQ=255)       aTCAGAGTTTTCACTCGATTTTGTGCTGATCAGTATTTTTTTTCTTTTTCCCCCGTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCggg                                                                                        <  2:347475/100‑1 (MQ=255)         ccaggttcgcACCCGATTTTGTGCTGTTCAGAATTTTTCTTCTTTTTCCCCCTTTAAGGGGCGAAGCCTCATCCCCCTTTCTCTGGTCGCCAGCCggg                                                                                        <  2:21905/89‑1 (MQ=255)         cAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTTTTTTTTTTTTTTGAAGAAGCGAAGCTTCATCCCTTTTTCTCTGCTCGCCAGCCGGGaa                                                                                      >  2:16110/1‑100 (MQ=255)            ggTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAAcc                                                                                   >  1:1417611/1‑100 (MQ=255)            ggTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAAcc                                                                                   <  1:1299311/100‑1 (MQ=255)             gtgttCAGCAGATGTGGTGCGCATCAGGATTTTCCTTCTTTTTCCTCCTGGAAGGGGCGAAGCCTCTTCCCCCCTTCTCTGGTCGCCAGCCGGGa                                                                                       <  2:1866700/92‑1 (MQ=255)             gTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCa                                                                                  <  1:1375638/100‑1 (MQ=255)                 tCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTa                                                                              <  1:1573939/100‑1 (MQ=255)                 tCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTa                                                                              <  2:2069700/100‑1 (MQ=255)                  cACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTaa                                                                             <  1:760620/100‑1 (MQ=255)                   aCCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCTCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAg                                                                            >  2:1926495/1‑100 (MQ=255)                   aCCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTa                                                                              <  1:606535/98‑1 (MQ=255)                    cccGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGc                                                                           >  1:2298287/1‑100 (MQ=255)                    cccGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGc                                                                           >  2:1798659/1‑100 (MQ=255)                     ccGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCt                                                                          >  1:1042527/1‑100 (MQ=255)                           ttGTGCTGATCAGAATTTTTTTTCTTTTTGCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTCGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCg                                                                    <  2:1417603/100‑1 (MQ=255)                              tgCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCa                                                                 <  1:658280/100‑1 (MQ=255)                               gCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGCACCCACGTAAGCTCCGGCGTccc                                                                >  1:2334083/1‑98 (MQ=255)                               gCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAc                                                                >  1:742660/1‑100 (MQ=255)                                  gATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGTAGCCTCATCCCCATTTCTCTGGTCACCTGCCGGGAAACCACGTAAGCTCCGGCGTCAccc                                                              >  1:870192/1‑99 (MQ=255)                                  gATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCa                                                             >  1:759731/1‑100 (MQ=255)                                  gATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCa                                                             <  1:706104/100‑1 (MQ=255)                                      aGAATTTTTTTTCTGTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAAc                                                         <  2:1602086/100‑1 (MQ=255)                                       gAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                        >  2:1603700/1‑100 (MQ=255)                                         aTTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGa                                                      <  2:742657/100‑1 (MQ=255)                                            ttttttCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCGGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATAc                                                   <  2:1697847/100‑1 (MQ=255)                                              ttttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACgg                                                 >  1:606511/1‑100 (MQ=255)                                              ttttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACgg                                                 <  1:2161081/100‑1 (MQ=255)                                              ttttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACgg                                                 >  1:1339111/1‑100 (MQ=255)                                              ttttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCGGGCGTCACCCATAACAGATACg                                                  >  2:1368872/1‑99 (MQ=255)                                               tttctttttcCCCGTTGAAGGGGCGAAGCCTCATCCCCGTTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTACGGCGTCACCCATAACAGATACGGa                                                >  2:1454451/1‑100 (MQ=255)                                                ttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGAc                                               >  2:622351/1‑100 (MQ=255)                                                 tctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACt                                              >  2:1899382/1‑100 (MQ=255)                                                 tctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACt                                              <  1:685886/100‑1 (MQ=255)                                                    ttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtc                                           <  1:2202131/100‑1 (MQ=255)                                                     tttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtct                                          <  1:6987/100‑1 (MQ=255)                                                      ttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtctc                                         <  1:1959723/100‑1 (MQ=255)                                                        cccccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTCGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaa                                       >  2:1262126/1‑100 (MQ=255)                                                        cccccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTCGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaa                                       >  2:1262158/1‑100 (MQ=255)                                                         ccccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaaa                                      >  2:1004635/1‑100 (MQ=255)                                                         ccccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCaaa                                      >  1:1985650/1‑100 (MQ=255)                                                          cccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAg                                     <  2:1751183/100‑1 (MQ=255)                                                          cccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAg                                     >  1:63111/1‑100 (MQ=255)                                                          accTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAg                                     >  2:2173201/2‑100 (MQ=255)                                                           ccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAgg                                    <  2:1872677/100‑1 (MQ=255)                                                            cTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGga                                   >  1:266585/1‑100 (MQ=255)                                                              tGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgaga                                 >  1:398490/1‑100 (MQ=255)                                                               gcaGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGagag                                >  1:1592434/3‑100 (MQ=255)                                                               gAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGagag                                >  2:588612/1‑100 (MQ=255)                                                                aaGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGt                               >  1:1790743/1‑100 (MQ=255)                                                                 aGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGtt                              >  1:870724/1‑100 (MQ=255)                                                                  ggggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGCTa                             >  2:2128975/1‑100 (MQ=255)                                                                  ggggCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCTAAGGAGAGTTa                             <  2:1806837/100‑1 (MQ=255)                                                                    ggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGTCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           >  2:1514035/1‑100 (MQ=255)                                                                    ggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           <  2:14209/100‑1 (MQ=255)                                                                    ggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           <  1:724332/100‑1 (MQ=255)                                                                    ggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           <  1:2216539/100‑1 (MQ=255)                                                                    ggCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAAGCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           >  1:1591608/1‑100 (MQ=255)                                                                    ggCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           >  2:1025177/1‑100 (MQ=255)                                                                    ggCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                           >  1:1418163/1‑100 (MQ=255)                                                                     gCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCa                          <  1:1125086/100‑1 (MQ=255)                                                                       gAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAAt                        <  1:1343515/100‑1 (MQ=255)                                                                        aaGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATg                       <  1:1389863/100‑1 (MQ=255)                                                                          gCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGaa                     >  2:1768589/1‑100 (MQ=255)                                                                           ccTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAat                    >  1:295063/1‑100 (MQ=255)                                                                            cTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAat                    >  1:1357128/1‑99 (MQ=255)                                                                            cTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAat                    <  2:2084773/99‑1 (MQ=255)                                                                             tCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAatat                  >  2:1724894/1‑100 (MQ=255)                                                                              cATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATAtt                 >  2:1328716/1‑100 (MQ=255)                                                                              cATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATAtt                 <  2:303282/100‑1 (MQ=255)                                                                               aTCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTc                <  2:1907387/100‑1 (MQ=255)                                                                                tCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCg               <  1:387331/100‑1 (MQ=255)                                                                                  cccATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc             <  1:442924/100‑1 (MQ=255)                                                                                  cccATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc             >  1:1786477/1‑100 (MQ=255)                                                                                  cccATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc             <  1:588615/100‑1 (MQ=255)                                                                                      tttCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCAtt         <  1:1644583/100‑1 (MQ=255)                                                                                        tctcTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGc       <  2:2093662/100‑1 (MQ=255)                                                                                         ctctGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCa      <  2:132225/100‑1 (MQ=255)                                                                                            tGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGa   >  1:315292/1‑100 (MQ=255)                                                                                             ggTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAt  >  1:819819/1‑100 (MQ=255)                                                                                             ggTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAt  >  2:167762/1‑100 (MQ=255)                                                                                             ggTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAt  >  2:87550/1‑100 (MQ=255)                                                                                             ggTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAt  >  2:1230259/1‑100 (MQ=255)                                                                                                 |                                                                                                AGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGAT  >  NC_000913/4370520‑4370711

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