Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,079,468	0	A	G	41.7%	98.2 / 13.9	60	D55G (GAT→GGT)	putP	proline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (5/30);  new base G (25/0);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 5.49e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACA  >  NC_000913/1079320‑1079617                                                                                                                                                     |                                                                                                                                                         ccGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGata                                                                                                                                                        <  1:645836/151‑1 (MQ=255)      gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGat                                                                                                                                                         >  1:77993/1‑145 (MQ=255)       ttGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGGGCg                                                                                                                                                  >  1:739781/1‑151 (MQ=255)                    gtgtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTaa                                                                                                                                    <  1:235386/150‑1 (MQ=255)                       tCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGAATGGACGgct‑ggat‑ggttattg                                                                                                                                 >  2:45901/1‑143 (MQ=255)                            atCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGATGAGGGGCT‑GGGC‑GGTTATGGGGTTg                                                                                                                            >  1:597303/1‑150 (MQ=255)                               tttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGTGGCGGGGGATTTGGGCggct‑ggcgg                                                                                                                                         >  1:434971/1‑133 (MQ=255)                                ttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCg                                                                                                                         <  1:527050/149‑1 (MQ=255)                                ttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGGGGGG‑GGGG‑GTTTATTGGGGTGGCggg                                                                                                                       >  1:309425/1‑151 (MQ=255)                                    cATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑Ct                                                                                                                  <  1:174436/151‑1 (MQ=255)                                    cATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑Ct                                                                                                                  <  1:501414/151‑1 (MQ=255)                                    cATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGATATGAGCGGGT‑GGCT‑GTTTATGGGGGTGGCggggg                                                                                                                     >  2:214740/1‑147 (MQ=255)                                     aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATGAGCGGCG‑GGCTGGTTATTGGGGTGGCggggg                                                                                                                      >  2:768670/1‑147 (MQ=255)                                       gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGtt                                                                                                               <  2:307107/151‑1 (MQ=255)                                       gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGt                                                                                                                <  2:609652/150‑1 (MQ=255)                                         tatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTTTGAGCGGCT‑GGGT‑GTTTTTGGGGGTGGCGGGGG‑CCGtttt                                                                                                             >  1:722973/1‑151 (MQ=255)                                             gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTtttctt                                                                                                         <  2:985129/151‑1 (MQ=255)                                             gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTtttctt                                                                                                         <  1:596648/151‑1 (MQ=255)                                             gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTtttctt                                                                                                         <  2:475373/151‑1 (MQ=255)                                              attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATTAGCGGGC‑GGGT‑GTTTATGGGGGTGGCGGGCG‑CTGTtttttttt                                                                                                        >  1:904011/1‑151 (MQ=255)                                                         cGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGCAAGGCT‑CGCG‑GTTTTTGGGGTGGCAGGGGG‑GCGTtttttttttctcgggt                                                                                                >  2:78662/1‑140 (MQ=255)                                                             tGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGTCGGGTATGAGCGGGT‑GGGG‑GTTTATGGGGGTGCCCGGGG‑CCGTTTTTTTTTTCCGGGTTTTCGaa                                                                                          >  1:665964/1‑150 (MQ=255)                                                             tGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGGT‑GGGT‑GGTTATTGGGGTGGCGGGGG‑CTGTTTTTTTTTTCCGGGTTTTCGaa                                                                                          >  1:208530/1‑150 (MQ=255)                                                                cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACTTCATTATCGGCGCGTTCGTCGTATATCATCGGCT‑GTCT‑TTTAATGGCTTTACCTGACG‑CTGTTTTTCTTTCCCTCATTTCCCAAATCt                                                                                      <  2:104339/151‑1 (MQ=255)                                                                cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGCT‑GGGT‑GTTAATGGGGGTGCCGGGGGCCGGTTTTTTTTTCGGGGGTTTCGGaa                                                                                          >  1:127573/1‑148 (MQ=255)                                                                cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGCC‑GGCTGGTTATTGGGGTTCCCGGG‑G‑GTGGTTTTTTTTTTGGGGTTTCCGaaaac                                                                                       >  2:384172/1‑148 (MQ=255)                                                                 gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTg                                                                                     <  1:591871/151‑1 (MQ=255)                                                                 gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTg                                                                                     <  2:248912/151‑1 (MQ=255)                                                                  aTCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTgg                                                                                    <  1:947674/151‑1 (MQ=255)                                                                   tCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGGGGGGGT‑GGGT‑GGTGATGGGGGTGGCGGGGG‑GGGGTTTTTTTTTCgggggttttgaaaacg                                                                                      >  2:855974/1‑135 (MQ=255)                                                                   tCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGGGGGGCT‑GGGT‑GGTTATGGGGGTGGCGGGGG‑CTGTTTTTCTTTTCGGGGTTTTCGGAAACTgg                                                                                    >  1:701931/1‑150 (MQ=255)                                                                    cAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGCGGGT‑GGCT‑GTTTATTGGGGTGGCGGGCG‑CTGTTTTTTTTTTCGGGGGTTTCGAAAGCTgg                                                                                    >  2:647222/1‑149 (MQ=255)                                                                      aCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCg                                                                                <  1:63126/151‑1 (MQ=255)                                                                      aCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCg                                                                                <  1:995388/151‑1 (MQ=255)                                                                       cGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGc                                                                               <  2:179391/151‑1 (MQ=255)                                                                          aaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTGGGGTATGAGGGGGT‑GGCT‑GTTTATTGGGGTGGCGGGGG‑CTGTTTTTTTTTTCGGGGGTTTCGAAAACTGGGTCGccc                                                                             >  1:569105/1‑149 (MQ=255)                                                                            aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGTTGAGCGGGTGGGGG‑GTTCTGGGGGTTGCGGGGGG‑GGGTTTTTTTTTCGGGGTTTTTC‑TAAGCTGGGTCGccct                                                                            >  2:381440/1‑147 (MQ=255)                                                                               tttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTg                                                                          <  2:721300/148‑1 (MQ=255)                                                                                 tGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGGGCGGCT‑GGCT‑GTTTATGGGGGTGGCGGGGG‑GGGTTTTTTTTTTCGGGGTTTTCGCAAGCAGGGGCGcccctggg                                                                        >  2:541519/1‑142 (MQ=255)                                                                                    cgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGGGGGT‑GGGT‑GTTTATGGGGGTGGCGGGCG‑GTGTTTTTTTTTTCGGGGTTTTCGAAAACTGGAttggc                                                                              >  1:965846/1‑135 (MQ=255)                                                                                     gacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGTTTGATAGGCT‑GGGT‑GTTTATTGGGTTGTCGGGGG‑CGGTTTTTTTTTTCGGGGTTTTCCAAAAACGGATGGCCCTTTGCCTGAc                                                                   >  1:811604/1‑149 (MQ=255)                                                                                        tataTTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGGT‑GGCT‑GTTAATGGGGGTGGCGGGGG‑CTGTTTTTTTTTTCGGGGTTTTCGAAAACTGGGTCGCCATTGGCCTGACAtt                                                                >  1:381279/1‑149 (MQ=255)                                                                                           atTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGGGGGT‑GGGT‑GTTTATTGGGGTGCCGGGCG‑CTGTTTTTTTTTTCGGGGTTTTCGAAAGCTGGGTCGCCCTTTGCCTGACATTAgggg                                                           >  1:342839/1‑149 (MQ=255)                                                                                             tCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGgcgc                                                          <  2:370386/150‑1 (MQ=255)                                                                                              cTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGCGGGG‑GGGT‑GGTTATGGGGGTTCCGGGCG‑GTGGTTTTTTTTTCGGGGGTTTCGAAAGCTGGGTCGCCCTTGGGCTGGCAtttagcggg                                                         >  2:679022/1‑143 (MQ=255)                                                                                                gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTgg                                                      <  1:92065/151‑1 (MQ=255)                                                                                                  gCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGAt                                                    <  2:906706/151‑1 (MQ=255)                                                                                                    ggTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTa                                                  <  2:868541/151‑1 (MQ=255)                                                                                                     gTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATGAGCGGCT‑GGGT‑GTTTATGGGGGTGGCCGGGG‑CGGGTTTTTTTTCCGGGATTTTCGAAAAGTGGAGCGCCCTCGGGCTGACATTCGGGGCGGGGGtt                                                   >  1:443/1‑149 (MQ=255)                                                                                                        gTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTg                                              <  2:791765/151‑1 (MQ=255)                                                                                                           gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGaa                                           <  2:806657/151‑1 (MQ=255)                                                                                                                 ggCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGGGCGGCC‑GGGT‑GTTTATTGGGGTGGCGGGGG‑CTGGTTTTTTTTTCGGGGTTTTCGAAAGCTGGATCGCCAATGGGCTGGCATTAGGGGCGGGGGTTTACttgga                                           >  2:640490/1‑141 (MQ=255)                                                                                                                   ccATTCGTGACGGCATTATCGGCGGGTGCGGCGGATATGAGCGGCT‑GGCT‑GTTTATGGGGGTGCCGGGGG‑CTGTTTTTTTTTTCGGGGTTTTCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGGGGGTTTACTGGAAGCtggggg                                   >  1:324768/1‑151 (MQ=255)                                                                                                                    cATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGc                                  <  1:799762/151‑1 (MQ=255)                                                                                                                           gACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCg                            <  2:874449/150‑1 (MQ=255)                                                                                                                                 aTTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgt                     <  1:668045/151‑1 (MQ=255)                                                                                                                                   tATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtg                    <  2:578637/150‑1 (MQ=255)                                                                                                                                    aTCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtgt                   <  1:276182/150‑1 (MQ=255)                                                                                                                                       ggcggGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCAt               >  2:483540/1‑151 (MQ=255)                                                                                                                                         cggGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAc             >  1:174319/1‑151 (MQ=255)                                                                                                                                                    gATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaaca  <  2:465416/151‑1 (MQ=255)                                                                                                                                                    gATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaaca  <  2:173153/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                         CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCT‑GGCT‑GTTAATGGGGTTGCCGGGCG‑CTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACA  >  NC_000913/1079320‑1079617

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