Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	628,652	0	T	G	44.8%	78.7 / 13.3	58	G34G (GGT→GGG)	entA	2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (10/22);  new base G (26/0);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 1.19e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTT  >  NC_000913/628504‑628798                                                                                                                                                     |                                                                                                                                                   ccGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGttt                                                                                                                                                  <  1:701794/151‑1 (MQ=255)   gACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGTGGAAGTTAAAGGgttt                                                                                                                                                 >  2:591784/1‑150 (MQ=255)        tCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGCGGGAAAGGTACAGGTTTTGGTc                                                                                                                                             >  1:445722/1‑149 (MQ=255)          gACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATc                                                                                                                                             >  2:734630/1‑147 (MQ=255)           aCGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGGGCGGAAGGTACAGGGTTTGATCaa                                                                                                                                           >  2:901474/1‑148 (MQ=255)           aCGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGCGCGGAAGTTACAGGGTTTGATCAAGCg                                                                                                                                        >  2:86625/1‑151 (MQ=255)              ccTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGATACAGGTTTTGATCAAGCGTTc                                                                                                                                     >  1:319353/1‑151 (MQ=255)                tggtggAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGTTCAc                                                                                                                                   >  1:128173/1‑151 (MQ=255)                    ggAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTTATCAAGCGTTCACTc                                                                                                                                 >  2:292329/1‑149 (MQ=255)                      aaGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGGTACAGGGTTTGTTCAAGCGGTCCCCCAgg                                                                                                                              >  2:440646/1‑150 (MQ=255)                       aGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAg                                                                                                                            <  1:658358/151‑1 (MQ=255)                          tACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGGGCGAAAGGTACAGGGTTTTATCAAGAGTTCACTCAGGAGCa                                                                                                                          >  2:402441/1‑150 (MQ=255)                          tACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGGAAGGTACAGGTTTTGATCAAACGTTCACTCAGGAg                                                                                                                            >  1:775591/1‑148 (MQ=255)                          tACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGGTACAGGGTTTGATCAAGAGTTCACTCAgg                                                                                                                              >  2:506600/1‑146 (MQ=255)                           aCTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCaa                                                                                                                         >  1:742609/1‑150 (MQ=255)                              ctcCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAata                                                                                                                       <  2:929487/149‑1 (MQ=255)                                    cgAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGGGAAAGGTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCctt                                                                                                                >  2:903891/1‑150 (MQ=255)                                            aataatGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGaa                                                                                                       <  2:857723/151‑1 (MQ=255)                                             ataatGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAg                                                                                                      <  2:200328/151‑1 (MQ=255)                                                atGGGTTTCCCCGGTAAAAAAGTGTTGGTAACCGGCGTTAGTAAAAGTTTTGGCTACGCCCCGGCGCTGGCTTTTTTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTg                                                                                                    <  2:791242/150‑1 (MQ=255)                                                atGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGGAAGGTACAGGGTTTGATCAAGCGTTCACTCAGGAGGAATATCCCTTTGCGACCGACGTg                                                                                                    >  2:651875/1‑150 (MQ=255)                                                 tGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGTTTACTCAGGAGCAATATCCCTTTGCGGCCGAAGGgat                                                                                                  >  2:705950/1‑151 (MQ=255)                                                         aGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTg                                                                                          <  1:373996/151‑1 (MQ=255)                                                          gCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGGTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGGGGTGGGTGtt                                                                                           >  2:600770/1‑149 (MQ=255)                                                            ggTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGCAAAATACAGGTTTTGATCAAGCGGTCACCCAGGAGGAATATCCCCTTGCGACCGGAAGGAGGGATGtttc                                                                                         >  2:1017169/1‑147 (MQ=255)                                                              tAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGGAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGGTGTTGcc                                                                                        >  1:754632/1‑148 (MQ=255)                                                               aaaaaTGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGGGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGAc                                                                                     >  1:644035/1‑150 (MQ=255)                                                                aaaaTGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTTGTCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGc                                                                                   >  2:503482/1‑151 (MQ=255)                                                                  aaTGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTg                                                                                 <  1:439730/151‑1 (MQ=255)                                                                   aTGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTgc                                                                                <  2:756189/151‑1 (MQ=255)                                                                      tCTGGGTAACCGGCGCAGGTAAAGGTATCGTCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTgcgc                                                                              >  1:439151/1‑150 (MQ=255)                                                                         gggTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAg                                                                            <  2:532686/149‑1 (MQ=255)                                                                          ggTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTc                                                                         <  2:1065779/151‑1 (MQ=255)                                                                             aaCCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGGAAAATACAGGGTTTGGTTAAGCGTTTACTTAGGAGGAATATCCCTTTGGGACCGCAGTGATGGATGGTGGCGACGCTGCGCAGGgcgc                                                                       >  2:801625/1‑150 (MQ=255)                                                                              aCCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGGcgcg                                                                      >  2:495821/1‑150 (MQ=255)                                                                                cGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCaa                                                                   <  2:419393/151‑1 (MQ=255)                                                                                cGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCCGCGCAGGgcgcgc                                                                     >  2:288178/1‑149 (MQ=255)                                                                                       ggTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGGTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAgtg                                                                >  2:189574/1‑147 (MQ=255)                                                                                          aaaGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGGTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGGTGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAg                                                          >  1:488366/1‑150 (MQ=255)                                                                                                aTCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGt                                                   <  1:343869/151‑1 (MQ=255)                                                                                                aTCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGGTGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGt                                                   >  1:705622/1‑151 (MQ=255)                                                                                                  cGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGtt                                                  <  2:89971/150‑1 (MQ=255)                                                                                                  cGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGtt                                                  <  1:966893/150‑1 (MQ=255)                                                                                                  cGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCCACGCTGCGCAGGGCGCGCCAGTGTGTCAGCGACTGtt                                                  >  1:240665/1‑150 (MQ=255)                                                                                                   ggCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCCTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGGCCCGCAAGTGTGTCAGCGACTGtt                                                  >  2:74359/1‑149 (MQ=255)                                                                                                      tACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGGTCACTCAGGAGCAATATCCCTTTGCGACCGAAGGGATGGATGTTGCCGACGCCGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTTGCTg                                             >  1:12256/1‑151 (MQ=255)                                                                                                       aCGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGa                                            <  1:270990/151‑1 (MQ=255)                                                                                                       aCGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTAAAGGGTTTTATCAAGCGTTCACTCCGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCCGCGCAGGGCGCGCAAGTGTGTCAGCGGCTGTTAGCTg                                             >  2:321784/1‑150 (MQ=255)                                                                                                          ccACGGCGCTGGCGTTTGTTGAGGCGGGAGCGGAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGCAGTGATGGGTGTTGCCCACGCTGCGCAGGGCGCGCAAGTGTGTCAGCGACTGTTAGCTGaaa                                          >  2:558804/1‑150 (MQ=255)                                                                                                           cACGGCGCTGGCGTTTGTTGAGGCGGGAGCGGAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGGCTGTTAGCTGAAAc                                         >  1:1047562/1‑150 (MQ=255)                                                                                                                cgcTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCg                                   <  2:217533/151‑1 (MQ=255)                                                                                                                       gTTTGTTGAGGCGGGGGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGAc                            <  1:367357/151‑1 (MQ=255)                                                                                                                       gTTTGTTGAGGCGGGAGCGGAAGGTACAGGTTTTGATCAAGCGTTCAATCAGGGGGAATATCCCTTTTCGACCGAAGTGATGGATGTTGCCCACGCTGCGCAGGTCGCGCAAGTGTGTCCGGGACTGTTAGCTGAAACGGGGCGACTgg                              >  2:189682/1‑149 (MQ=255)                                                                                                                       gTTTGTTGAGGCGGGAGCGAAAGTTACAGGGTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGGGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGa                             >  1:639480/1‑150 (MQ=255)                                                                                                                          tgttgAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGAcgc                          <  2:535495/150‑1 (MQ=255)                                                                                                                           gttgAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGAcgcg                         <  2:369387/150‑1 (MQ=255)                                                                                                                                        cGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATgcggc           <  2:775468/151‑1 (MQ=255)                                                                                                                                            aGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGa       <  2:492764/151‑1 (MQ=255)                                                                                                                                              ttACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGaa      >  1:526653/1‑150 (MQ=255)                                                                                                                                                 cAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAAtttt  <  2:847834/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CCGACCATCGACGCCTGGTGGAAGCTACTCTCCCGCGAGGTGAAATAATGGATTTCAGCGGTAAAAATGTCTGGGTAACCGGCGCAGGTAAAGGTATCGGCTACGCCACGGCGCTGGCGTTTGTTGAGGCGGGAGCGAAAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTT  >  NC_000913/628504‑628798

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