Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,359,827	0	C	T	45.3%	20.1 / 108.6	75	G69S (GGC→AGC)	cadB	putative lysine/cadaverine transporter
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (18/16);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

ATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATT  >  NC_000913/4359684‑4359977                                                                                                                                                |                                                                                                                                                       aTTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCttgt                                                                                                                                                 <  1:28694/151‑1 (MQ=255)   ttAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCttgttg                                                                                                                                               >  1:1099428/1‑151 (MQ=255)     aaTACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCg                                                                                                                                             >  2:1021006/1‑151 (MQ=255)      aTACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCgg                                                                                                                                            >  2:320170/1‑151 (MQ=255)      aTACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCgg                                                                                                                                            >  1:915536/1‑151 (MQ=255)            ggAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTCCCTTGTTGCGGGttttt                                                                                                                                      <  1:1235556/151‑1 (MQ=255)                       gAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTc                                                                                                                           >  1:996570/1‑151 (MQ=255)                        aaaGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCg                                                                                                                          >  1:301780/1‑151 (MQ=255)                            ataagataCAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGc                                                                                                                       >  2:11489/1‑150 (MQ=255)                                 ataCAGCGGTAATACCAATCGCCAGGTTACCAAGCCAGTTAGCATGGTAATAAAGAACACCTGTCTGTAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATatac                                                                                                                 <  1:276002/151‑1 (MQ=255)                                   aCAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATatacat                                                                                                               >  2:1247314/1‑151 (MQ=255)                                     aGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATatacatac                                                                                                             <  2:107674/151‑1 (MQ=255)                                        ggTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGcc                                                                                                          >  1:723934/1‑151 (MQ=255)                                         gTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCa                                                                                                         <  1:208919/151‑1 (MQ=255)                                          tAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGcc                                                                                                          >  1:42922/1‑149 (MQ=255)                                           aaTACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGc                                                                                                       <  2:371022/151‑1 (MQ=255)                                                cAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACAt                                                                                                  >  2:425613/1‑151 (MQ=255)                                                  aTCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTg                                                                                                <  1:769917/151‑1 (MQ=255)                                                  aTCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTg                                                                                                <  1:703429/151‑1 (MQ=255)                                                    cGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGc                                                                                               >  1:886813/1‑150 (MQ=255)                                                      ccAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCAcc                                                                                            >  1:427470/1‑151 (MQ=255)                                                       cAGGTTACCCAACCAGTTAGCATGGAAATAAAGAACACCTGTCTGAAAACCAAATGCAGGAGAACTTTATCCGGCATAAGCAATCGGGCTACCTTGTTGCGTGTTTTTTGATGCCAGTCGGGCATATACATCTGCCATCGACATTGCAc                                                                                             >  1:886074/1‑149 (MQ=255)                                                       cAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCAcc                                                                                            >  1:37223/1‑150 (MQ=255)                                                          gTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCaata                                                                                        <  2:258808/151‑1 (MQ=255)                                                            tACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCaataa                                                                                       <  1:114464/150‑1 (MQ=255)                                                             aCCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAataata                                                                                     <  2:857837/151‑1 (MQ=255)                                                               cAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCaataat                                                                                      >  2:414824/1‑148 (MQ=255)                                                                aaTCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATag                                                                                    >  1:978870/1‑149 (MQ=255)                                                                 aTCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAgaga                                                                                 <  1:67495/151‑1 (MQ=255)                                                                  tCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATag                                                                                    >  1:497021/1‑147 (MQ=255)                                                                    cAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATa                                                                               >  1:522123/1‑150 (MQ=255)                                                                     aGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAAt                                                                             <  2:1058858/151‑1 (MQ=255)                                                                        tAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATcc                                                                           <  2:745012/150‑1 (MQ=255)                                                                              ggTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCa                                                                    <  1:539130/151‑1 (MQ=255)                                                                               gTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAAcccc                                                                     >  1:1291945/1‑149 (MQ=255)                                                                                taataaAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAAcccc                                                                     >  2:1271826/1‑148 (MQ=255)                                                                                 aataaAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAAcccc                                                                     >  1:102331/1‑147 (MQ=255)                                                                                  ataaAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGAt                                                                 >  1:969088/1‑150 (MQ=255)                                                                                  aaacagaagaCCTGTCTGAAAGCTAAATGCAGGGGAAATTTCGCTGGCATAAGCAATTGGTCTACCTTGATGCGGGTTTTTTGCTGCCAGTCGGGCATATACATACGCCAGCTACATTGCACCCATAATAGTGATAATCCAACCCCAGAt                                                                 >  1:1151727/5‑150 (MQ=255)                                                                                   taaAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGAt                                                                 >  2:1233682/1‑149 (MQ=255)                                                                                      aGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCaa                                                            <  2:265367/151‑1 (MQ=255)                                                                                       gAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAAt                                                           >  2:1213600/1‑151 (MQ=255)                                                                                          cacCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATacc                                                        >  1:452474/1‑151 (MQ=255)                                                                                           acCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATacc                                                        <  2:648332/150‑1 (MQ=255)                                                                                            ccTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATaccac                                                      >  1:621925/1‑151 (MQ=255)                                                                                               gTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCg                                                    <  1:214131/150‑1 (MQ=255)                                                                                                   gAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGAt                                                  <  2:1042755/148‑1 (MQ=255)                                                                                                   gAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGAt                                                  <  1:543380/148‑1 (MQ=255)                                                                                                    aaacCCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATCATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACtt                                              >  1:801573/1‑151 (MQ=255)                                                                                                      aaCCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGc                                            <  2:701658/151‑1 (MQ=255)                                                                                                        ccAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTa                                          <  2:1099467/151‑1 (MQ=255)                                                                                                        ccAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCt                                           >  1:1278937/1‑150 (MQ=255)                                                                                                         cAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAg                                         <  1:145625/151‑1 (MQ=255)                                                                                                            aTGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGtt                                      <  1:1258389/151‑1 (MQ=255)                                                                                                            aTGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGtt                                      >  1:481239/1‑151 (MQ=255)                                                                                                            aTGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGtt                                      <  1:881302/151‑1 (MQ=255)                                                                                                              gCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTc                                     <  2:913590/150‑1 (MQ=255)                                                                                                               cAGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGc                                   >  2:546183/1‑151 (MQ=255)                                                                                                                aGGGGAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCa                                  <  1:13602/151‑1 (MQ=255)                                                                                                                    gAAATTTCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGta                              <  1:911974/151‑1 (MQ=255)                                                                                                                    gAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGt                               >  1:688840/1‑150 (MQ=255)                                                                                                                         ttCTCCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTaataat                         >  1:427274/1‑151 (MQ=255)                                                                                                                          tctcCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATg                        >  2:1231945/1‑151 (MQ=255)                                                                                                                           ctcCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGc                       >  2:863711/1‑151 (MQ=255)                                                                                                                            tcCGGCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGc                       >  2:361687/1‑150 (MQ=255)                                                                                                                             ccGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATg                        <  2:874787/148‑1 (MQ=255)                                                                                                                              cGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAAt                    <  1:213065/151‑1 (MQ=255)                                                                                                                               ggCATAAGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATa                   <  2:886853/151‑1 (MQ=255)                                                                                                                                    aaGCAATTGGGCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGc               >  2:437643/1‑150 (MQ=255)                                                                                                                                    aaGCAATCGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATa                   <  1:1142514/146‑1 (MQ=255)                                                                                                                                    aaGAAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCGAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATTGCCATACCACCGATACTTGCTATGATCGCAGGTAATAATGCAATACCGc               >  1:710890/1‑150 (MQ=255)                                                                                                                                            gggCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatca      <  2:411641/151‑1 (MQ=255)                                                                                                                                            gggCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatca      <  1:1020968/151‑1 (MQ=255)                                                                                                                                            gggCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatc       <  2:1095156/150‑1 (MQ=255)                                                                                                                                             ggCTACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatcat     <  1:1017759/151‑1 (MQ=255)                                                                                                                                               ccACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCatat   >  1:1036162/1‑151 (MQ=255)                                                                                                                                                cACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATAtt  >  1:976236/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                       ATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATT  >  NC_000913/4359684‑4359977

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