Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,703,598	C→T	17.2%	R237R (CGG→CGA)	ydgJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,703,598	0	C	T	17.2%	226.0 / 34.4	99	R237R (CGG→CGA)	ydgJ	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (52/30);  new base T (13/4);  total (65/34)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

AGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTC  >  NC_000913/1703501‑1703695                                                                                                  |                                                                                                  aGTCTTCCTGCGGTAGACGCGCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCg                                                                                                 <  1:1955677/100‑1 (MQ=255)    cTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGt                                                                                                >  1:4225819/1‑98 (MQ=255)     ttCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGc                                                                                              <  1:3235560/99‑1 (MQ=255)     ttCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCt                                                                                             <  1:3587005/100‑1 (MQ=255)       ccTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGa                                                                                           <  1:3379149/100‑1 (MQ=255)        cTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGcatgg                                                                                          <  1:2365351/100‑5 (MQ=255)        cTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGAc                                                                                          >  2:3847394/1‑100 (MQ=255)        cTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGAc                                                                                          <  1:1400297/100‑1 (MQ=255)        cTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGAc                                                                                          <  1:1420863/100‑1 (MQ=255)         tGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGGTATATCGTGcatg                                                                                           <  1:1744711/98‑4 (MQ=255)           cGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCa                                                                                       >  2:2780049/1‑100 (MQ=255)            ggTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGACAGAGGCCATATTTGAAATAACTGCCTGGGGATGCATGCACGAAGTACCGTGCTGACt                                                                                         >  2:3947668/1‑97 (MQ=255)             gTAGACGCTCGCCATTTGTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCTCATAACTGCCTCGGTATCCATGCACGGTATATCGTGcatgggtccc                                                                                     <  1:373254/100‑10 (MQ=255)              tAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCg                                                                                    >  1:4048564/1‑100 (MQ=255)               aGACGCTCGCCATTTTTCGGACGTTCTGCCGGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCgg                                                                                   <  1:4304371/100‑1 (MQ=255)               aGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCgg                                                                                   <  1:2046210/100‑1 (MQ=255)               aGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCgg                                                                                   <  2:319537/100‑1 (MQ=255)                  cccTCGCCATTTTTCAGACGTTCTTCCTGGGGCTCGAGGCCATATTTCACATACCTCCCTCGGGATCCATGCACGTTATACCGTGCTGACTCAGCGGCTg                                                                                <  1:112577/98‑1 (MQ=255)                    cTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGcc                                                                              >  1:4019660/1‑100 (MQ=255)                    cTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGcc                                                                              <  1:3519677/100‑1 (MQ=255)                    cTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGcc                                                                              <  1:2513150/100‑1 (MQ=255)                     tCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCa                                                                             <  1:4399128/100‑1 (MQ=255)                       gCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGGGGCTGCCAGc                                                                           >  2:946959/1‑100 (MQ=255)                        ccATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCa                                                                          >  1:1553421/1‑100 (MQ=255)                          aTTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATg                                                                        >  2:1750170/1‑100 (MQ=255)                              ttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTAc                                                                    <  1:3506927/100‑1 (MQ=255)                              ttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTAc                                                                    <  1:3386453/100‑1 (MQ=255)                              ttCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTAc                                                                    >  1:541544/1‑100 (MQ=255)                                     gTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCCGCATGGTACCGTGTaa                                                             >  2:3499881/1‑100 (MQ=255)                                      ttctccCTGTGGATCAAGGCCATATTCCACATACCTGCCACGAGCCCCATGCACGATATACCGTGCTCACTCCTCGGCTGCCAGCATGGTACCGTGTaa                                                             <  1:3013243/99‑1 (MQ=255)                                        cttcCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAAt                                                          >  2:3946627/1‑100 (MQ=255)                                         ttcCTGTGGATCGAGGCCATATGTCACATGACTGCCGCGGGATCCGTGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCATGTaaaa                                                           <  1:3766004/98‑1 (MQ=255)                                            cTGTGGATCGAGGCCATATTTCACATAACTGCCGCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTATAAAATGAc                                                       >  1:2889351/1‑99 (MQ=255)                                             tgtgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTc                                                     >  2:1717918/1‑100 (MQ=255)                                             tgtgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTc                                                     >  2:3763934/1‑100 (MQ=255)                                              gtgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTCCCGTGTAAAATGACTcg                                                    >  1:4053412/1‑100 (MQ=255)                                               tgGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgc                                                   >  1:1923059/1‑100 (MQ=255)                                                 gATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATCCCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccg                                                 >  1:1456411/1‑100 (MQ=255)                                                 gATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccg                                                 >  1:945105/1‑100 (MQ=255)                                                 gATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccg                                                 >  1:4281303/1‑100 (MQ=255)                                                  ntCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGCAAAATGACTcgccac                                                >  1:2669253/2‑98 (MQ=255)                                                  aTCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTcgccgc                                                >  2:3349271/1‑100 (MQ=255)                                                     gAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCtgt                                             >  1:673473/1‑100 (MQ=255)                                                     gAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCtgt                                             <  2:1923066/100‑1 (MQ=255)                                                       ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTgg                                           >  1:112990/1‑100 (MQ=255)                                                       ggCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTgg                                           >  2:35516/1‑100 (MQ=255)                                                        gCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTCACTCAGCGGCTGCCAGCATGGTACCCTGCAAAATGACTCGCCGCt                                               >  1:2119017/1‑95 (MQ=255)                                                          cATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGGATGGTACCGGGTAAAATGATTCGCCGCTGTgg                                           >  1:2598724/1‑97 (MQ=255)                                                            tatTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAgg                                      <  1:2142065/100‑1 (MQ=255)                                                             atTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGa                                     <  1:833123/100‑1 (MQ=255)                                                              tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAATATGACTCGCCGCTGTGGATAGGAc                                    >  1:3263091/1‑100 (MQ=255)                                                              tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGAc                                    >  2:3246962/1‑100 (MQ=255)                                                              tttCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGAc                                    >  2:3626561/1‑100 (MQ=255)                                                                 cacaTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAg                                 <  1:1890581/100‑1 (MQ=255)                                                                    aTAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGa                                <  1:2874551/98‑1 (MQ=255)                                                                     tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATgg                             >  2:1700036/1‑100 (MQ=255)                                                                     tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATgg                             >  1:3596555/1‑100 (MQ=255)                                                                     tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                              >  1:3718174/1‑99 (MQ=255)                                                                     tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATg                              >  1:4068178/1‑99 (MQ=255)                                                                     tAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGGGTAAAATGACTCGCCGCCGTGGATAGGACa                                   >  1:3713639/1‑94 (MQ=255)                                                                         tGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTg                         >  1:3836509/1‑100 (MQ=255)                                                                          gCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTgg                        >  2:3979506/1‑100 (MQ=255)                                                                          gCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTgg                        >  1:1808172/1‑100 (MQ=255)                                                                           ccTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGa                       <  2:4019690/100‑1 (MQ=255)                                                                            cTCGGGATCCATGCACGATATACCGTGCTGACTCAGTGGCTGCCAGCATGGTACCCTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTg                         <  1:2413126/97‑1 (MQ=255)                                                                             tCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGaaa                     >  1:1590765/1‑100 (MQ=255)                                                                             tCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGa                       <  1:2738944/98‑1 (MQ=255)                                                                               gggATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAaata                   >  2:918302/1‑100 (MQ=255)                                                                               gctgagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAaata                   >  2:3604573/10‑100 (MQ=255)                                                                               gctgagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACGCGCCGCTGTGGATAGGACAAGATGGCGTGGAaata                   >  1:2462578/10‑100 (MQ=255)                                                                                ggATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGAATGGTACCGTGTAAAATAACTCGCCAGTGTGGATCGGACAAGGTGGCGTAGAaata                   >  2:2569618/1‑99 (MQ=255)                                                                                ctgagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTCAAATGACTCGCCGCTGTGGATAGGACAAAATGGCGTGGAaataa                  >  1:4320748/9‑100 (MQ=255)                                                                                ctgagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAaataa                  >  2:336121/9‑100 (MQ=255)                                                                                 gATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAaataat                 >  2:1409568/1‑100 (MQ=255)                                                                                  gagtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATc                >  2:3460941/7‑100 (MQ=255)                                                                                  aTCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATc                <  2:1553427/100‑1 (MQ=255)                                                                                   tCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg               >  1:1023631/1‑100 (MQ=255)                                                                                   tCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg               <  1:3786790/100‑1 (MQ=255)                                                                                   agtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg               >  2:495916/6‑100 (MQ=255)                                                                                   agtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg               >  2:494647/6‑100 (MQ=255)                                                                                   agtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCg               >  2:1572273/6‑100 (MQ=255)                                                                                   agtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATAGCGTGGAaat                    >  1:2738054/6‑95 (MQ=255)                                                                                    gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATTGTACCGTGTAAAATGACTCGCCGCTGGGGATAGGACAAGATGGCGTGGAAATAATCg               >  1:2444502/5‑99 (MQ=255)                                                                                    gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCgg              >  1:960091/5‑100 (MQ=255)                                                                                    gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCgg              >  2:258772/5‑100 (MQ=255)                                                                                    gtcagCACGGTATATCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGGGGAAATAATCgg              >  2:2776740/5‑100 (MQ=255)                                                                                    ccATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTAACGTGTAAAAGGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCgg              <  1:1050865/100‑1 (MQ=255)                                                                                      aTGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTc            <  1:1087290/100‑1 (MQ=255)                                                                                       ngCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCg           >  2:2198970/2‑100 (MQ=255)                                                                                        gCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGa          >  1:1291993/1‑100 (MQ=255)                                                                                        gCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGa          >  2:765128/1‑100 (MQ=255)                                                                                         cACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGAc         >  2:2104904/1‑100 (MQ=255)                                                                                         cACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGAc         >  2:3377849/1‑100 (MQ=255)                                                                                          aCGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACt        >  2:2223712/1‑100 (MQ=255)                                                                                          aCGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACt        >  2:2223628/1‑100 (MQ=255)                                                                                           cGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTg       <  1:2484916/100‑1 (MQ=255)                                                                                              tataCCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTgcgc    >  1:1985935/1‑100 (MQ=255)                                                                                              tataCCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTgc      <  1:1711568/98‑1 (MQ=255)                                                                                                 aCCGTGCTGACTCAGCAGCTGCCAGCATGGTACCGTGTCAAATGACGCGCCTCTGTGGATAGGACAAGATGACGTGAAAATAATCGGTCGACTGCGCTc  >  1:4317017/1‑99 (MQ=255)                                                                                                  |                                                                                                  AGTCTTCCTGCGGTAGACGCTCGCCATTTTTCAGACGTTCTTCCTGTGGATCGAGGCCATATTTCACATAACTGCCTCGGGATCCATGCACGATATACCGTGCTGACTCAGCGGCTGCCAGCATGGTACCGTGTAAAATGACTCGCCGCTGTGGATAGGACAAGATGGCGTGGAAATAATCGGTCGACTGCGCTC  >  NC_000913/1703501‑1703695

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