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.4%	88.0 / 19.9	58	D55G (GAT→GGT)	putP	proline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (1/33);  new base G (24/0);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.95e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCG  >  NC_000913/1079320‑1079608                                                                                                                                                     |                                                                                                                                                ccGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGata                                                                                                                                               <  2:17416/151‑1 (MQ=255)    aTGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGa                                                                                                                                            <  1:211983/151‑1 (MQ=255)           tGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAgg‑gg                                                                                                                                       >  1:443014/1‑148 (MQ=255)             aCATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑ggctggc                                                                                                                                  <  1:237590/151‑1 (MQ=255)               aTTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGTGGCGGattataaagggtgg                                                                                                                                    >  1:1078734/1‑136 (MQ=255)                      gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATAAGAGC‑GGCTGGCT‑GTTAAttg                                                                                                                        <  2:397429/151‑3 (MQ=255)                        cTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGTGGGTATTAGC‑GGCGGGCT‑GTTGATtggg                                                                                                                      >  1:343051/1‑151 (MQ=255)                          atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATgggg                                                                                                                      <  1:26589/149‑1 (MQ=255)                           tatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGA‑gggtggggggttattgg                                                                                                                        >  1:293502/1‑136 (MQ=255)                            atCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTg                                                                                                                   <  2:1098276/150‑1 (MQ=255)                            atCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTg                                                                                                                   <  2:259515/150‑1 (MQ=255)                            atCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGC‑GGGCGGGT‑GTTTATgggg                                                                                                                      >  1:244427/1‑147 (MQ=255)                                ttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGc                                                                                                                  <  2:1135793/147‑1 (MQ=255)                                 tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCggg                                                                                                              <  2:344729/150‑1 (MQ=255)                                  ggCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGTCGGGTATGAGC‑GGCTGGGT‑GTTTATTGGGGTGGCggggg                                                                                                            >  1:785583/1‑149 (MQ=255)                                      tGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGG‑GGCTGGGT‑GTTTATGGGGGTGCCGGGGGCTg                                                                                                         >  1:930818/1‑150 (MQ=255)                                       gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGCGCGCTG‑tt                                                                                                      <  2:520971/151‑1 (MQ=255)                                       gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGTGC‑GGGCGGCTGGTTATTGGGGTTGCggggg                                                                                                             >  2:396554/1‑145 (MQ=255)                                          attgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑ttttt                                                                                                   <  2:853432/151‑1 (MQ=255)                                           ttgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGC‑GGGTGGCC‑GGTTATGGGGGTGGCGGGGGGTG‑ttttt                                                                                                   >  2:1147112/1‑150 (MQ=255)                                            tgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGG‑GGGGGGGT‑GTTTATGGGGGGGGCGGGCGGTG‑ttttttt                                                                                                 >  1:618738/1‑149 (MQ=255)                                             gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTtttct                                                                                                 <  2:291198/150‑1 (MQ=255)                                              attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGG‑GGGCGGCT‑GGTTATGGGGGTGGCGGGGGCTG‑Gtttttttt                                                                                               >  1:407332/1‑151 (MQ=255)                                                  ggTTTATCGGCTGGCGATCAACGAAAAACTTTGACGACGATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCgg                                                                                           <  2:461204/151‑1 (MQ=255)                                                    tttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGAGG‑GGGTGGCT‑GTTTATGGGGGTGCCGGGGGCTG‑TTTTTTTTTTCgggg                                                                                         >  1:454000/1‑150 (MQ=255)                                                     ttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGATGAGC‑GGCTGGGC‑GTTTATGGGGGTGGCGGGCGGTG‑GTTTTTTTTTCgggg                                                                                         >  1:131906/1‑149 (MQ=255)                                                      tATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGAt                                                                                        <  1:1164826/150‑1 (MQ=255)                                                         cGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTcc                                                                                    <  2:752862/151‑1 (MQ=255)                                                              ggCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGaaa                                                                                <  1:355582/150‑1 (MQ=255)                                                                  aTCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTgg                                                                           <  1:866322/151‑1 (MQ=255)                                                                    cAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGa                                                                          <  1:565158/150‑1 (MQ=255)                                                                     aaCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGAt                                                                         <  2:1031343/150‑1 (MQ=255)                                                                         aaaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGGGG‑GGGGGGGC‑GTTTATTGGGGTGTCGGGGGCGG‑TTTTTTTTTCCGGGGGTTCCGAAAGAAGGATCg                                                                       >  1:636281/1‑148 (MQ=255)                                                                           aaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCAtt                                                                  <  1:265919/151‑1 (MQ=255)                                                                              cTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGG‑GGGTGGGT‑GTTTATTGGGGTGCCGGGGGGTG‑TTTTTTTTTTCGGGGTTTTCGGAAACTGGGTCGCCCTTgg                                                                >  1:1125100/1‑150 (MQ=255)                                                                                ttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGc                                                               <  1:906672/149‑1 (MQ=255)                                                                                  gacgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGa                                                           <  2:346671/151‑1 (MQ=255)                                                                                   acgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGCGGGCGGGCG‑GTT‑ATGGGGGTGGCGGGGGCGG‑GTTTTTTTTTTTGGGTTTTCGTAAGGTGGGTCggccctgggccgtc                                                          >  2:908284/1‑138 (MQ=255)                                                                                         atatTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTa                                                      <  1:550729/149‑1 (MQ=255)                                                                                            ttCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGG‑GGGTGGGT‑GTTTATGGGGGTGCCGGGCGCTG‑TTTTTTTTTTCCGGGTTTTCGCAAGCTGGGTCGCCATTGGCCTGACATTTGgcg                                                  >  1:436500/1‑150 (MQ=255)                                                                                            ttCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGGATGAGG‑GGGGGGGC‑GTTAATGGGGGTGGCGGGCGCTG‑GTTTTTTTTTCGGGGGTTTCGGAAGCTGGGTTGCCCTTGGCCCGACCTTAggg                                                   >  1:570944/1‑148 (MQ=255)                                                                                             tCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGgcgc                                                 <  2:660971/150‑1 (MQ=255)                                                                                                gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTgg                                             <  1:1171866/151‑1 (MQ=255)                                                                                                gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTg                                              <  1:619119/150‑1 (MQ=255)                                                                                                gggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGAGC‑GGGTGGCG‑GTTTATGGGGGTGCCGGGGGGCG‑TTTTTTTTTTCCGGGTTTTCGGAAGCTGGAACGCCCTTTGGCCGACATTAAgcgggggg                                             >  1:1012392/1‑146 (MQ=255)                                                                                                  gCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGAt                                           <  1:25035/151‑1 (MQ=255)                                                                                                    ggTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTTTGTGC‑GGGTGGCTGGTTATTGGGTTGGCGGGGGGCGG‑TTTTTTTTTTTGGG‑TTTTCGAAAACTGGGTTGCCACTTGCCTGTGCTTAGgcggggggg                                            >  2:926340/1‑142 (MQ=255)                                                                                                         tAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGA‑GGGTGGGT‑GGTTATTGGGGTGGCGGGGGCGGTTTTTTTTTTCCGGGTTTTCGGAAA‑ATGGATCGCCCTTTGCCCGACATTTGGCGCGGGGGTTAAATgg                                    >  1:1133453/1‑151 (MQ=255)                                                                                                          aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTgg                                    <  1:1137843/150‑1 (MQ=255)                                                                                                          aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGG‑GGGTGGGT‑GTTAATGGGGGGGGCGGGGGGGG‑TTTTTTTTTTCGGGGGTTTCGAAAACGGGATCGCCACTGGCCTGACATTAGGGGCGGGGGTTAACTgg                                    >  1:335039/1‑150 (MQ=255)                                                                                                           gTCTTGGGGCATGATTGACGGCATTAACGGCGGGTGCGTCGGATATGGTC‑GGCTGGCT‑GTTAATGGGGTTGGCGGGGGCTG‑TTTGTCTTTCCGGGATTTCCGAAAGCTGGATTGCCATTAGCCTGACATTAGGCGTGTGGATTAACTGGaa                                  <  2:340793/151‑1 (MQ=255)                                                                                                              ttGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGAGC‑GGCTGGCT‑GTTAATGGGGGTGGCGGGGGGTG‑TTTTTTTTTTCCGGGTTTTTGCAAGCTGGATCGCCATTTGCCTGACATTAGGGGCGTGGGTTAACTggga                                  >  1:1089207/1‑146 (MQ=255)                                                                                                                    cATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTTTGAGG‑GGCTGGCG‑GTTTATGTGGGTGGCGGGCGCCG‑TTTTTTTTTTCGGGGGTTTCGGAAACGGGATCGCCCCTTGCCTGGGCTTAGGCGCGGGGGTTTACTGGAAGCgggtgg                          >  1:183001/1‑150 (MQ=255)                                                                                                                       tCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCg                       <  2:148884/150‑1 (MQ=255)                                                                                                                        cGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCggg                     <  1:542711/151‑1 (MQ=255)                                                                                                                                   tATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtg           <  1:889432/150‑1 (MQ=255)                                                                                                                                    aTCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtgt          <  2:812854/150‑1 (MQ=255)                                                                                                                                         cggGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGTGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGCTTTACGAAAGCTGGATCGCCATTGGCCTGACATTACGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGGTGCGTGTGCa       >  2:396007/1‑148 (MQ=255)                                                                                                                                          ggggGCGCCGGGTATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAACGCTGGATCGCCATTGGCCTGACCTTACGCGCGTGGATTAACTGGCAGCCGGCGGCCGGGCGGTTGCGGGTGCATAcc   >  2:213983/1‑151 (MQ=255)                                                                                                                                           ggTGCGTCGGATATGAGC‑GGCTGGCT‑GGTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTGCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATTCCg  <  2:28056/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGC‑GGCTGGCT‑GTTAATGGGGTTGCCGGGCGCTG‑TTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCG  >  NC_000913/1079320‑1079608

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