Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,258,974	0	T	G	41.2%	117.4 / 12.1	68	G411G (GGT→GGG)	dinF	oxidative stress resistance protein; putative MATE family efflux pump; UV and mitomycin C inducible protein
Reads supporting (aligned to +/- strand):  ref base T (11/29);  new base G (28/0);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.90e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGA  >  NC_000913/4258827‑4259121                                                                                                                                                     |                                                                                                                                                    ccAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGtttt                                                                                                                                                  <  2:561105/151‑1 (MQ=255)    gATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGgttt                                                                                                                                                  >  2:657766/1‑148 (MQ=255)    gATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTgcg                                                                                                                                               <  1:466765/151‑1 (MQ=255)     aTTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCACCGGtttt                                                                                                                                                  >  2:1104389/1‑147 (MQ=255)              ctggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTTCGCTGACGcccc                                                                                                                                     >  1:1195985/1‑151 (MQ=255)               tggctgACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCCCCCCAGGGTTTTCGCTGACGcccc                                                                                                                                     >  1:656091/1‑150 (MQ=255)                    gACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGCGGGTGGCC‑CCGCAGGGGTTGCGCTGACGCTCCTTACGc                                                                                                                               >  2:790670/1‑151 (MQ=255)                    gACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTAc                                                                                                                                 <  2:1126783/149‑1 (MQ=255)                     aCCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCCCAGGGTTTGCGCTGCCGCTCCTTTCGc                                                                                                                               >  2:846948/1‑150 (MQ=255)                        gCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGGTTTTCGCTGACGCTCCTTACgccgcc                                                                                                                           >  2:957500/1‑151 (MQ=255)                        gCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGGGGCCGCCGCAGGGTGTGCGCTGCCGCCCCTTACg                                                                                                                                >  2:816390/1‑146 (MQ=255)                           atcttatctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGGCGCTCCTTACGCTGCCGTg                                                                                                                        >  2:818838/1‑151 (MQ=255)                                atctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGTGGCTgg                                                                                                                   >  1:233104/1‑151 (MQ=255)                                 tctGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCCTGCGGTCCTTCCGCCGCCGTGGCTggg                                                                                                                  >  2:744093/1‑151 (MQ=255)                                    ggCAGGTGATTTTGCCGGTGGTCGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTaa                                                                                                               <  2:1020083/151‑1 (MQ=255)                                       aGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGCGGCCGCCGCAGGGTTTGCGCTGACCCTCCTTTCGCTGCCGTGGCTGGGTaa                                                                                                               >  1:1186223/1‑148 (MQ=255)                                          tGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCa                                                                                                            >  1:797536/1‑148 (MQ=255)                                           gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGc                                                                                                         <  1:227807/150‑1 (MQ=255)                                           gATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCt                                                                                                        <  2:958673/151‑1 (MQ=255)                                             ttttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCttt                                                                                                      >  1:231097/1‑151 (MQ=255)                                              tttGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCttt                                                                                                      >  1:1146961/1‑150 (MQ=255)                                                tGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGGTTTTCGCGGACGCCCCTTTCGCTGCCGGGGGTGGGGAATCATGCtttt                                                                                                     >  2:408053/1‑148 (MQ=255)                                                tGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTtgt                                                                                                    <  2:340461/150‑1 (MQ=255)                                                tGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGGAATCATGCTTtgt                                                                                                    >  1:1014309/1‑150 (MQ=255)                                                tGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAACGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCCTACGCTGCCGTGGGTGGGGAATCATGCTTtttg                                                                                                   >  2:101708/1‑151 (MQ=255)                                                    ggtggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGGTTTTCGCTGACGCTCCTTTCGCTGCCGTGGCTGGGGAATTAAGCTTTGtgg                                                                                                  >  2:1105663/1‑148 (MQ=255)                                                       ggtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGtggctggc                                                                                             <  1:1085091/150‑1 (MQ=255)                                                        gtTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCAt                                                                                           <  2:874990/151‑1 (MQ=255)                                                           ggCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCCCCGCCGGGTTTTCGCCGACGCTCCTTACGCCGCCGCGGCTGGGGAATCATGCTTTGTGGCTGGCAtt                                                                                          >  2:171279/1‑149 (MQ=255)                                                            gCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGACGCCCCTTTCGCCGCCGTGGCTGGGGAATCATGCTTTGTGGCTGGCATTaa                                                                                        >  2:110205/1‑150 (MQ=255)                                                             cGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAAcc                                                                                      <  1:557732/151‑1 (MQ=255)                                                              gTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCg                                                                                     <  2:235051/151‑1 (MQ=255)                                                               tCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCGCGGCTGGGGAATCATGCTTTGTGGCTGGCATTAAcc                                                                                      >  1:733259/1‑149 (MQ=255)                                                                   gtgtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTCGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCttt                                                                                <  2:141911/151‑1 (MQ=255)                                                                     gtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTc                                                                               <  2:218559/150‑1 (MQ=255)                                                                     gtTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACCCTGCCGTGGCTGGGGAAACCTGCTTTGTGGCTGGCCTTTAACGTCttt                                                                                >  1:9786/1‑149 (MQ=255)                                                                       tATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTg                                                                             <  1:436818/150‑1 (MQ=255)                                                                       tATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTg                                                                             <  2:414848/150‑1 (MQ=255)                                                                       tATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTg                                                                             <  2:592812/150‑1 (MQ=255)                                                                       tATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTgg                                                                            >  1:343251/1‑151 (MQ=255)                                                                        aTCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGGTTTGCGCTGACGCTCCTTACGCTGCCCTGGCTGGGGAATCATGCTTTGTGGCTGGCATTAACCGTCttttttg                                                                            >  1:687269/1‑146 (MQ=255)                                                                           tgctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCCCCGCAGGGTTTGCGCTGAAGCTCCCTACGCTGGCGTGGGTGGGGAATCATGTTTTGGGGCTGGCCTTTAACGTCttt                                                                                >  2:1168372/1‑143 (MQ=255)                                                                            gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCg                                                                          <  1:207606/148‑1 (MQ=255)                                                                            gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCg                                                                          <  1:428329/148‑1 (MQ=255)                                                                            gctgGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGtt                                                                        <  2:802455/150‑1 (MQ=255)                                                                               ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCCGGGTTTGCGCTGACGCTCCTTACGCCGCCGGGGCTGGGGAATCATGCTTTTTGGGTGGCCTTAAACCTCTTTTTTGCggt                                                                        >  2:1111932/1‑145 (MQ=255)                                                                               ggACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTgc                                                                      <  2:1072976/149‑1 (MQ=255)                                                                                 aCGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGGGGCCGCCGCCCGTTTTTCGCTGAAGCTCCTTACGCTGCCGTGGCGGGGTAAATCTACTTTGGGGCTTGCAATAACCGTCTTTTTCGCGTTTCGCgg                                                                  >  1:791351/1‑151 (MQ=255)                                                                                   ggCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCGCAGGTTTTTTGGTGACGCCCCTTACGCTGCCGTGGCTGGGTAATCATCCTTTTTGGCCGGCGTTAACCGGCTTTCTTGCGGGGCGCcgg                                                                 >  2:161618/1‑150 (MQ=255)                                                                                   ggCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGc                                                                <  1:28125/151‑1 (MQ=255)                                                                                       tGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTTCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAAACGTCTTTTTGGCGTTGCGCGGGCttt                                                             >  1:58086/1‑150 (MQ=255)                                                                                        gTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCCCAGGGTTTGCGCTGACGCCCCTTACGCTGCCGTGGCTGGGGAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGGGCGGGCttttt                                                           >  2:652954/1‑149 (MQ=255)                                                                                           tataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTctct                                                         <  1:980586/150‑1 (MQ=255)                                                                                            ataGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTg                                                        <  1:883517/150‑1 (MQ=255)                                                                                                  gcAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCCCCGCCGGGTTTGCGCCTACGCTCCCTACGCTGCCGCGGCCGGGTAAACATTCTTTTTGGGTGGCACTAACCGGCCTTCTGGCGGTGCGCGCGCCTTCTCTggcg                                                     >  2:226853/1‑147 (MQ=255)                                                                                                  gcAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGGGGGCGCCGCCGGGTTTGCGCCGACGCTCCTTATGCCGCCGTGGGTGGGGAATCATGCTTTGTGGCCGGCATTAAACGTCTTTTTTGGGGTGCGCGGGGTTTCTCTtgcgggtt                                                 >  2:204087/1‑148 (MQ=255)                                                                                                    aaCGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCCGGTTTTGCGCTGACCCCCCTTACGCTGCCGGGGCTGGGGAATTATGCTTTTTGGCCGGGCTTAACCGGCTTTTTTGGGTTGCGCGGGGTTTTTCTgggggccctt                                               >  2:1089781/1‑147 (MQ=255)                                                                                                       gcgTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTgg                                            <  1:1016124/151‑1 (MQ=255)                                                                                                         gTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGGCGCCCCCGGGTTTTCGCCGACGCTCCTTTAGCTGCCGCGGCTGGGGAATCATGCTTTGTGGCGGGGCTTAAACGGCTTTTTGGCGTTTCGGGGGGTTTCTCTGGCGGCCTTCTtgcg                                          >  2:88327/1‑151 (MQ=255)                                                                                                          tGCCACCGAAATGCGTAACAGTATGGC‑GGTTGCCCCCCCAGGGTTTTCGCTGACGCCCCTTACGCCGCCGTGGCCGGGGAATCCTGCTTTTTGGGTGGCATTTACCGTCTTTCTGGCGGTTCGCGGGGTTTTTCTTGGGGCCATc                                               >  2:1061565/1‑145 (MQ=255)                                                                                                             cACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGCGGCTGGGGAAACATTCTTTTTGGCTGGCATTAACCGTCTTTTTTGCGTTGCGCGGGGTTTCTCTGGCGGGTATTTGgcggcg                                       >  2:823798/1‑150 (MQ=255)                                                                                                                 gAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCCCCTTACGCTGCCGGGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGGCTTTTTGGCGTTGCGCGGGGTTTTTCTGGCGGCTATCTGGCGGCGTCAc                                   >  1:1126377/1‑150 (MQ=255)                                                                                                                  aaaTGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTg                                 <  2:981177/151‑1 (MQ=255)                                                                                                                            cAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATgg                        >  1:340487/1‑150 (MQ=255)                                                                                                                              gTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTAc                     >  2:514310/1‑151 (MQ=255)                                                                                                                                aTGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCt                   <  2:963339/151‑1 (MQ=255)                                                                                                                                  ggC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTgg                 <  1:818859/151‑1 (MQ=255)                                                                                                                                  ggC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTgg                 <  1:389391/151‑1 (MQ=255)                                                                                                                                             gccgcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCa        >  2:520312/1‑150 (MQ=255)                                                                                                                                                gcAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGt    <  2:797516/151‑1 (MQ=255)                                                                                                                                                   ggTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGa  <  2:915833/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCAGATTCAGCAGCTGGCTGACCGCTATCTTATCTGGCAGGTGATTTTGCCGGTGGTTGGCGTCTGGTGTTATCTGCTGGACGGCATGTTTATAGGCGCAACGCGTGCCACCGAAATGCGTAACAGTATGGC‑GGTGGCCGCCGCAGGTTTTGCGCTGACGCTCCTTACGCTGCCGTGGCTGGGTAATCATGCTTTGTGGCTGGCATTAACCGTCTTTCTGGCGTTGCGCGGGCTTTCTCTGGCGGCTATCTGGCGGCGTCACTGGCGCAATGGTACCTGGTTTGCCGCAACGTGA  >  NC_000913/4258827‑4259121

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