Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,777	0	G	T	40.7%	82.5 / 24.4	54	Q245K (CAA→AAA)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base G (10/22);  new base T (22/0);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 9.26e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGT  >  NC_000913/646629‑646917                                                                                                                                                       |                                                                                                                                             aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                                             <  2:12422/150‑1 (MQ=255) aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                                             <  1:262223/150‑1 (MQ=255)  atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAAGCCCCCCTTTTTCTATAAATTTTtcgc                                                                                                                                           >  2:504476/1‑151 (MQ=255)  atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGGCCGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTgcgc                                                                                                                                           <  1:116918/151‑1 (MQ=255)     aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAATGCCCCCCTTTTTCAATAAATTTTTCGcccc                                                                                                                                        >  2:978019/1‑149 (MQ=255)         cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATGCCCCCCTTTTTCAATAATGTTTTCGCCTcg                                                                                                                                       >  2:105297/1‑148 (MQ=255)           gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTcgcg                                                                                                                                     <  1:893141/148‑1 (MQ=255)             gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTa                                                                                                                                <  2:534012/151‑1 (MQ=255)             gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTa                                                                                                                                <  2:518716/151‑1 (MQ=255)                aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GGATTCCCCCCTTTTTCAATAATTGTTTTGCCCCGCGCCGTAGcc                                                                                                                             >  2:576193/1‑151 (MQ=255)                   aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgc                                                                                                                          <  1:685019/151‑1 (MQ=255)                   aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGC‑GAATGCCCCCCTTTTTCAATAATTTTTTGGCCTCGCGCTGGAGcc                                                                                                                             >  2:100095/1‑148 (MQ=255)                      cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCa                                                                                                                       <  1:406425/151‑1 (MQ=255)                      cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTG‑GAAAGCCCCCCTTTTTCAATAATGTTTGCGGCTCGCGCTGTa                                                                                                                                >  2:665179/1‑142 (MQ=255)                        ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATGCCCCCCTTTTTCAATAATTTTTGTGCCTCCCGCTGTAGCCAgcgc                                                                                                                        >  2:192622/1‑148 (MQ=255)                           aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                                                                   <  2:352867/150‑1 (MQ=255)                              gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcccccc                                                                                                               <  1:43285/151‑1 (MQ=255)                               aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATGCCCCCCTTTTTCTATAATGTTTTTGCCTCGCGCTGTTGCCACCGCCGGccccccc                                                                                                              >  1:444368/1‑150 (MQ=255)                                 cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCg                                                                                                            <  1:934465/151‑1 (MQ=255)                                  aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGGCGCCCCCCTTTTGTAAAAATATTTTAGCCTCGCGGTGTAGTCCGCGCGGGcccccccccc                                                                                                           >  1:756207/1‑147 (MQ=255)                                   gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAAAGCCCCCCTTTTTCAATAAATTTTTTGCCCCCCCCCGTAGCCCGCGCAGGccccccccccc                                                                                                          >  2:742690/1‑146 (MQ=255)                                     aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTG‑AAAGGCCCCCTTTTTGC‑AGAAAGTTTTCGCCCCCCGCCGCGGCCACCCCCGGccccccccccccc                                                                                                        >  1:167029/1‑144 (MQ=255)                                     aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGCCGGCT‑GAATGCCCCCTTTTTGCAACAATGGTTTAGCCTCCCGCCGTTGACTCCGCCCGcccccccctgccc                                                                                                        >  2:221380/1‑144 (MQ=255)                                               tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATTCCCCCCTTTTTTAATAAAGTTTTTGCCCCCCCCCGGTGCCCGCGCCGGcccccccccccccccc                                                                                                     >  1:574330/1‑134 (MQ=255)                                                 ccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATGCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCCGTAGCCAGCGCAGGGCCCCCCCCCCCCCCCCCGCTGCa                                                                                             >  2:656080/1‑150 (MQ=255)                                                   gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCCTTTTTTAATAATGTTTTCGCCCCCCGCCGTAGCCCGCGCCGGGCCCCCCCGCCCCCCCGCGATGCAAc                                                                                           >  1:382275/1‑150 (MQ=255)                                                    ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATGCCCCCCTTTTTTAAAAATGTTTTTGCCCCCCGCCGTAGCCCGCGCCGGGCCCCCCCCCCCCCCCCCGATGCaa                                                                                            >  2:615007/1‑148 (MQ=255)                                                        gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAAGGCCCCCTTTTTCAATAATGTTTTTGCCCTGCGCCGTAGCCCGCGCCGGGCCCCCTCCCCCCCCCGCGaag                                                                                               >  2:810201/1‑139 (MQ=255)                                                          cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGAATGCCCCCTTTTT‑CAATAATGTTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCACCCCCCGATGCaaagcttgt                                                                                     >  2:494478/1‑142 (MQ=255)                                                            gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAAGCCCCCCTTTTTGAATAATGTTTTCCCCCCGCCCTGTAGCCAGCGCAGGGCCCCCCCCCCACCCCGCGCTGCACAGCTGGTGTCg                                                                                 >  2:502026/1‑151 (MQ=255)                                                                  tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                                            <  1:813032/150‑1 (MQ=255)                                                                       cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAATGCCCCCTTTTTTCAATAATGTTTTTGCCCCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCACCCCGCGGTGCCACGTTGGGATTGCCGGTGATCg                                                                       >  2:836040/1‑150 (MQ=255)                                                                       cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATAGCCCCTTTTTGCAATAAATTTTTCGCCTCCCGCTGTAGCCAACGCCGGCCCCCCCCCCCCCCGCCCGGTGCAACCTTGGGATCGcccgtgtt                                                                         >  2:222017/1‑142 (MQ=255)                                                                         aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                                                                    <  2:611960/151‑1 (MQ=255)                                                                              cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GCAGTCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGCTGCAACCTTGGTGTTGCCGCTGATCGccccccgt                                                               >  1:76266/1‑145 (MQ=255)                                                                                cccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATGGCCCCCTTTTGTCATAAAGTTTTCGCCCCCCGCCGTTGCCCGCGCCGGCCCCCCTCGCCCCCCCGCGAGGCAAAGTTGGGATCGCCGTTGATTGCCCTCAGt                                                               >  2:367757/1‑149 (MQ=255)                                                                                  cTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATTCCCCCCTTTTTTCATAAAGTTTTCGCCCCGCCCTGTAGGCCGCGCACGCCCCCCTCCCCACCCCGCGAGGCCACGGTGGGATCGCCCtttgttgcc                                                                     >  2:113302/1‑134 (MQ=255)                                                                                    gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                                          <  1:977882/150‑1 (MQ=255)                                                                                      cGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATGGCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCCGTAGCCCGCGCCGGGCCCCCCCCCCCCCCCCCCATGCAACGGTGGGATCGCCGGTGTTCGCCCCCAGTAgg                                                            >  1:794355/1‑146 (MQ=255)                                                                                                gagaTAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                                             <  1:724086/151‑1 (MQ=255)                                                                                                 agaTAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCa                                            <  1:939477/151‑1 (MQ=255)                                                                                                           gaTCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc                                  <  2:343679/151‑1 (MQ=255)                                                                                                              cggc‑ggGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCCGGCCCCCCTCGCCACCCCGCGATGCAACGTTGGGATCGCCGTTGATCGCCATCAGTAGGGGCAAGGGATCGAGCCATGCCAGTTCAgg                               >  1:40644/1‑151 (MQ=255)                                                                                                                   gggggTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt                           <  2:713756/151‑1 (MQ=255)                                                                                                                    ggggTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa                          <  2:240913/151‑1 (MQ=255)                                                                                                                          c‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGGTGATCGCCATCAGTAGGAGCAAGGAATCGAGCAATGCCAGTTCAGGATCTAAcccc                     >  2:290243/1‑150 (MQ=255)                                                                                                                            gAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTg                   <  2:410426/151‑1 (MQ=255)                                                                                                                            gAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTg                   >  2:51165/1‑151 (MQ=255)                                                                                                                            gAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTg                   >  2:481128/1‑151 (MQ=255)                                                                                                                             aaTGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTg                   >  2:977093/1‑150 (MQ=255)                                                                                                                                gCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATcc               >  1:907532/1‑151 (MQ=255)                                                                                                                                gCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATc                <  1:373759/150‑1 (MQ=255)                                                                                                                                  ccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcag             <  2:982/151‑1 (MQ=255)                                                                                                                                      tttttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcagc            >  1:111221/1‑148 (MQ=255)                                                                                                                                        tttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCagag        <  2:395363/150‑1 (MQ=255)                                                                                                                                         ttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTg      >  1:36114/1‑151 (MQ=255)                                                                                                                                             aataatGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGt  <  2:502855/151‑1 (MQ=255)                                                                                                                                             aataatGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGt  <  1:100078/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                             AATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGCAGAGTGAGGT  >  NC_000913/646629‑646917

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