Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,747	0	T	G	44.3%	100.7 / 13.1	61	T255P (ACC→CCC)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (4/30);  new base G (27/0);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 2.25e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTA  >  NC_000913/646599‑646893                                                                                                                                                      |                                                                                                                                                     aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGtt                                                                                                                                                     <  1:957127/150‑1 (MQ=255) aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTc                                                                                                                                                    <  1:695396/151‑1 (MQ=255)    atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑Ga                                                                                                                                                 <  2:248582/150‑1 (MQ=255)    atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑Ga                                                                                                                                                 <  1:756546/150‑1 (MQ=255)      ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAt                                                                                                                                               <  2:597755/150‑1 (MQ=255)       tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑g                                                                                                                                             <  2:956156/150‑1 (MQ=255)         cAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑Gccc                                                                                                                                          <  2:841586/151‑1 (MQ=255)         cAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAAT‑Gccc                                                                                                                                          >  2:373825/1‑151 (MQ=255)            aTTTATCAAGAGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAGT‑Gcccccc                                                                                                                                       >  1:741966/1‑150 (MQ=255)              ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCtt                                                                                                                                      <  2:494614/150‑1 (MQ=255)                aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCttttt                                                                                                                                   <  2:922139/151‑1 (MQ=255)                    aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCaa                                                                                                                               <  1:296177/151‑1 (MQ=255)                    aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCaa                                                                                                                               <  2:891781/151‑1 (MQ=255)                      tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAATGGCCCCCtttttca                                                                                                                                 >  2:549236/1‑146 (MQ=255)                        cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑TCCCCCcttttttcataa                                                                                                                            >  1:176275/1‑143 (MQ=255)                         ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAA‑GCCCCCTTTTTGtaataa                                                                                                                            >  2:314628/1‑149 (MQ=255)                          tAAATAATTAAATCAGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCaataat                                                                                                                           <  1:150973/149‑1 (MQ=255)                               aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTgc                                                                                                                     <  1:262223/150‑1 (MQ=255)                               aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTgc                                                                                                                     <  2:12422/150‑1 (MQ=255)                                atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAA‑GCCCCCCTTTTTCTATAAATTTTtcgc                                                                                                                   >  2:504476/1‑151 (MQ=255)                                atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGGCCGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTgcgc                                                                                                                   <  1:116918/151‑1 (MQ=255)                                   aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAAT‑GCCCCCCTTTTTCAATAAATTTTTCGcccc                                                                                                                >  2:978019/1‑149 (MQ=255)                                       cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAT‑GCCCCCCTTTTTCAATAATGTTTTCGCCTcg                                                                                                               >  2:105297/1‑148 (MQ=255)                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTcgcg                                                                                                             <  1:893141/148‑1 (MQ=255)                                           gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTa                                                                                                        <  2:534012/151‑1 (MQ=255)                                           gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTa                                                                                                        <  2:518716/151‑1 (MQ=255)                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GGAT‑TCCCCCCTTTTTCAATAATTGTTTTGCCCCGCGCCGTAGcc                                                                                                     >  2:576193/1‑151 (MQ=255)                                                 aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgc                                                                                                  <  1:685019/151‑1 (MQ=255)                                                 aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGC‑GAAT‑GCCCCCCTTTTTCAATAATTTTTTGGCCTCGCGCTGGAGcc                                                                                                     >  2:100095/1‑148 (MQ=255)                                                    cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCa                                                                                               <  1:406425/151‑1 (MQ=255)                                                    cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTG‑GAAA‑GCCCCCCTTTTTCAATAATGTTTGCGGCTCGCGCTGTa                                                                                                        >  2:665179/1‑142 (MQ=255)                                                      ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGAT‑GCCCCCCTTTTTCAATAATTTTTGTGCCTCCCGCTGTAGCCAgcgc                                                                                                >  2:192622/1‑148 (MQ=255)                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                                           <  2:352867/150‑1 (MQ=255)                                                            gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcccccc                                                                                       <  1:43285/151‑1 (MQ=255)                                                             aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGAT‑GCCCCCCTTTTTCTATAATGTTTTTGCCTCGCGCTGTTGCCACCGCCGGccccccc                                                                                      >  1:444368/1‑150 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCg                                                                                    <  1:934465/151‑1 (MQ=255)                                                                aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGGC‑GCCCCCCTTTTGTAAAAATATTTTAGCCTCGCGGTGTAGTCCGCGCGGGcccccccccc                                                                                   >  1:756207/1‑147 (MQ=255)                                                                 gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAAA‑GCCCCCCTTTTTCAATAAATTTTTTGCCCCCCCCCGTAGCCCGCGCAGGccccccccccc                                                                                  >  2:742690/1‑146 (MQ=255)                                                                   aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTG‑AAAG‑GCCCCCTTTTTGC‑AGAAAGTTTTCGCCCCCCGCCGCGGCCACCCCCGGccccccccccccc                                                                                >  1:167029/1‑144 (MQ=255)                                                                   aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGCCGGCT‑GAAT‑GCCCCCTTTTTGCAACAATGGTTTAGCCTCCCGCCGTTGACTCCGCCCGcccccccctgccc                                                                                >  2:221380/1‑144 (MQ=255)                                                                             tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGAT‑TCCCCCCTTTTTTAATAAAGTTTTTGCCCCCCCCCGGTGCCCGCGCCGGcccccccccccccccc                                                                             >  1:574330/1‑134 (MQ=255)                                                                               ccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAT‑GCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCCGTAGCCAGCGCAGGGCCCCCCCCCCCCCCCCCGCTGCa                                                                     >  2:656080/1‑150 (MQ=255)                                                                                 gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCCTTTTTTAATAATGTTTTCGCCCCCCGCCGTAGCCCGCGCCGGGCCCCCCCGCCCCCCCGCGATGCAAc                                                                   >  1:382275/1‑150 (MQ=255)                                                                                  ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGAT‑GCCCCCCTTTTTTAAAAATGTTTTTGCCCCCCGCCGTAGCCCGCGCCGGGCCCCCCCCCCCCCCCCCGATGCaa                                                                    >  2:615007/1‑148 (MQ=255)                                                                                      gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAA‑GGCCCCCTTTTTCAATAATGTTTTTGCCCTGCGCCGTAGCCCGCGCCGGGCCCCCTCCCCCCCCCGCGaag                                                                       >  2:810201/1‑139 (MQ=255)                                                                                        cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGAAT‑GCCCCCTTTTT‑CAATAATGTTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCACCCCCCGATGCaaagcttgt                                                             >  2:494478/1‑142 (MQ=255)                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAA‑GCCCCCCTTTTTGAATAATGTTTTCCCCCCGCCCTGTAGCCAGCGCAGGGCCCCCCCCCCACCCCGCGCTGCACAGCTGGTGTCg                                                         >  2:502026/1‑151 (MQ=255)                                                                                                tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                    <  1:813032/150‑1 (MQ=255)                                                                                                     cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAAT‑GCCCCCTTTTTTCAATAATGTTTTTGCCCCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCACCCCGCGGTGCCACGTTGGGATTGCCGGTGATCg                                               >  2:836040/1‑150 (MQ=255)                                                                                                     cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGAT‑AGCCCCTTTTTGCAATAAATTTTTCGCCTCCCGCTGTAGCCAACGCCGGCCCCCCCCCCCCCCGCCCGGTGCAACCTTGGGATCGcccgtgtt                                                 >  2:222017/1‑142 (MQ=255)                                                                                                       aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                                            <  2:611960/151‑1 (MQ=255)                                                                                                            cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GCAG‑TCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGCTGCAACCTTGGTGTTGCCGCTGATCGccccccgt                                       >  1:76266/1‑145 (MQ=255)                                                                                                              cccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAT‑GGCCCCCTTTTGTCATAAAGTTTTCGCCCCCCGCCGTTGCCCGCGCCGGCCCCCCTCGCCCCCCCGCGAGGCAAAGTTGGGATCGCCGTTGATTGCCCTCAGt                                       >  2:367757/1‑149 (MQ=255)                                                                                                                cTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAT‑TCCCCCCTTTTTTCATAAAGTTTTCGCCCCGCCCTGTAGGCCGCGCACGCCCCCCTCCCCACCCCGCGAGGCCACGGTGGGATCGCCCtttgttgcc                                             >  2:113302/1‑134 (MQ=255)                                                                                                                  gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                  <  1:977882/150‑1 (MQ=255)                                                                                                                    cGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAT‑GGCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCCGTAGCCCGCGCCGGGCCCCCCCCCCCCCCCCCCATGCAACGGTGGGATCGCCGGTGTTCGCCCCCAGTAgg                                    >  1:794355/1‑146 (MQ=255)                                                                                                                              gagaTAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                     <  1:724086/151‑1 (MQ=255)                                                                                                                               agaTAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCa                    <  1:939477/151‑1 (MQ=255)                                                                                                                                         gaTCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc          <  2:343679/151‑1 (MQ=255)                                                                                                                                            cggc‑ggGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCCGGCCCCCCTCGCCACCCCGCGATGCAACGTTGGGATCGCCGTTGATCGCCATCAGTAGGGGCAAGGGATCGAGCCATGCCAGTTCAgg       >  1:40644/1‑151 (MQ=255)                                                                                                                                                 gggggTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt   <  2:713756/151‑1 (MQ=255)                                                                                                                                                  ggggTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa  <  2:240913/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTA  >  NC_000913/646599‑646893

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