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

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAGGATCTAA  >  NC_000913/646599‑646894                                                                                                                                                     |                                                                                                                                                       aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACTCTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTc                                                                                                                                                      <  1:752330/151‑1 (MQ=255)     tttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCg                                                                                                                                                     >  2:547070/1‑148 (MQ=255)      ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                                  <  1:986460/150‑1 (MQ=255)       tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCATg                                                                                                                                                 >  1:511091/1‑150 (MQ=255)           aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑ccccc                                                                                                                                           <  2:825579/151‑1 (MQ=255)             tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑ccccc                                                                                                                                           >  2:470055/1‑149 (MQ=255)               tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAAG‑CCCCCcttt                                                                                                                                       >  2:595347/1‑151 (MQ=255)                aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCttttt                                                                                                                                      <  1:1022054/151‑1 (MQ=255)                 tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTGTGAATG‑CCCCCcttttt                                                                                                                                     >  2:641328/1‑150 (MQ=255)                 tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCctttt                                                                                                                                      >  1:1129697/1‑150 (MQ=255)                   aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCcttttt                                                                                                                                     >  2:350383/1‑148 (MQ=255)                    aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCcttt                                                                                                                                       >  1:81483/1‑146 (MQ=255)                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCaataat                                                                                                                              <  2:816363/149‑1 (MQ=255)                          tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCaataat                                                                                                                              <  2:399686/149‑1 (MQ=255)                            aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCTATAATtttt                                                                                                                          >  2:692220/1‑151 (MQ=255)                             ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTg                                                                                                                         <  1:258805/151‑1 (MQ=255)                             ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCTTTTTTCAATAATttttt                                                                                                                         >  2:379986/1‑150 (MQ=255)                              taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATG‑CCCCCTTTTTTCCATAATttttt                                                                                                                         >  2:522725/1‑149 (MQ=255)                              taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAG‑GCCCCCTTTTTCAAAAATGtttt                                                                                                                         >  1:612014/1‑149 (MQ=255)                                    aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCTTTTTTCAATAATGTTTTCGCCTc                                                                                                                   >  2:302069/1‑150 (MQ=255)                                          tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAG‑CCCCCTTTTTGCCATAAAGATTTCGCCCCccgc                                                                                                               >  1:49539/1‑148 (MQ=255)                                          tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAG‑GCCCCCTTTTGCAATAAAGTTTTTGCCTccc                                                                                                                 >  2:157772/1‑144 (MQ=255)                                           gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGt                                                                                                            <  2:902574/150‑1 (MQ=255)                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAg                                                                                                          <  2:705450/149‑1 (MQ=255)                                                 aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGTAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAgc                                                                                                    <  1:376796/151‑1 (MQ=255)                                                 aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAG‑CCCCCCTTTTGCTATAAAGTTTTCGCCCCGCGCCGTAG‑CCAgc                                                                                                    >  2:1113523/1‑151 (MQ=255)                                                    cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATG‑GCCCCCTTTTTTAATAATGTTTTCGCCCCGCGCTGTAG‑CCa                                                                                                      >  1:421171/1‑146 (MQ=255)                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATG‑CCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCCGTAG‑CCAGCGCAGGccc                                                                                            >  2:1132719/1‑151 (MQ=255)                                                          aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGcccc                                                                                           <  2:837463/151‑1 (MQ=255)                                                            gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑GCCCCCTTTTTTAATAAAGTTTTGGCCCCCCGCCGTTG‑GCAGCGCAGGgccccc                                                                                         >  1:1100845/1‑151 (MQ=255)                                                                      tCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTGC‑ATAATGTTTTCGCCCCCCGCCGTAG‑CCCGCGCAGGCCCCCCCCGCcccccc                                                                               >  1:701585/1‑149 (MQ=255)                                                                       cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATG‑CCCCCCTTTTTCAATAATGTTTTTGCCCCCCGCTGTAG‑CCCGCGCAGGCCCCCCCCCCCAcccc                                                                               >  2:87128/1‑148 (MQ=255)                                                                             tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTTAATAAAATTTTCGCCCCCCGCTGTAG‑CCAGCGCAGGcccccccccccccccccc                                                                             >  2:592486/1‑134 (MQ=255)                                                                               ccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCaa                                                                      <  2:81488/151‑1 (MQ=255)                                                                                cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAAc                                                                     <  2:89172/151‑1 (MQ=255)                                                                                    tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGCGGAGCGGCTG‑CCCCCCTTTTTTAGTAATGCTTGTGCCGCTCGTTGTTGACCGCCGCA‑GCCCCCCCCCCCACcgcgtcaagccacagt                                                                  >  1:523811/1‑138 (MQ=255)                                                                                        cTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCTTTTTTCAATAATGTTTTCGCCCCCCGCCGTTG‑CCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCAAAGTTgg                                                                >  1:578624/1‑148 (MQ=255)                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                            <  2:732422/150‑1 (MQ=255)                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCg                                                           <  1:260396/151‑1 (MQ=255)                                                                                               ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCTTTTTTTAATAATGTTTTCGCCTCCCGCCGTTG‑ACAGCGCAGGCCCCCCCCCCCACCCCCCGATGCAACGTTGGGATCGcc                                                         >  1:971249/1‑148 (MQ=255)                                                                                                tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                      <  2:960347/150‑1 (MQ=255)                                                                                                  gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGAt                                                   <  1:470032/151‑1 (MQ=255)                                                                                                    tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                                  <  2:553371/150‑1 (MQ=255)                                                                                                     cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGCGGGTCCAATG‑CCCCCCTTTTTCCATAATTTTTTCGCCTCCCGCTGTCG‑CCACCGCAGGCCCCCCCCCCCCCCCCGCCGTGCCACGCTGGTATCCCCCTTggt                                                   >  1:121509/1‑146 (MQ=255)                                                                                                       aTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCa                                              <  2:719101/151‑1 (MQ=255)                                                                                                         acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCCAATG‑CCCCCCTTTTTCAATAAAGTTTTTGCCTCCCCCTGTAG‑GCCGCGCAGGCCCCCCCCCCCCCCCCCCGGGTCAACGCTGGTATTGCCGGTGGTTGcc                                               >  1:143176/1‑148 (MQ=255)                                                                                                           acTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAg                                          <  2:182201/151‑1 (MQ=255)                                                                                                             tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCCTTTTTTAATAATGTTTTCGCCCCCCGCCGTAG‑CCCGCGCAGGCCCCCCCCCCCCCCCCGCGGTGCAACGTTGGTATTGCCGGTGATCGCCCTc                                            >  2:392654/1‑147 (MQ=255)                                                                                                                 tGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCAATAATTTTTTCGCCTCCCGCCGTAG‑CCCGCGCAGGCCCCCCCCGCCCCCCCGCGCGGCAACGGTGGGATCGCCGGTGATCGCCATCAGTAgg                                      >  2:248032/1‑149 (MQ=255)                                                                                                                  gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGc                                   <  2:589702/151‑1 (MQ=255)                                                                                                                   tCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTTGAGG‑CCCCCCTTTTTCTGCAATGATTTCTCCCCCCGCCGGGG‑TCTGCGCAGGCGCCCCCCCTCCCCCCCCGGTGCGACGGTAGTATCCCCGTTCAtgcccaccaatatgt                                     >  1:766055/1‑134 (MQ=255)                                                                                                                     gAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAg                                <  1:909855/151‑1 (MQ=255)                                                                                                                     gAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATG‑CCCCCCTTTTTCAATAAAGTTTTTGCCCCCCGCCGTAG‑CCAGCGCAGGCCCCCCCCGCCACCCCCCGGTGCAACCGTGGTATCCCCCCTGATCGCCATCAGTAGGAg                                    >  2:508950/1‑147 (MQ=255)                                                                                                                        cTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATT‑CCCCCTTTTTTCAATAAAGTTTTCCCCCCGCCCTGTAG‑CCAGCGCGGGCACCCCTCGCCACACCGCGGTGCAACGTTGGTATTGGCGTTGATCGCCCTCCGTAGGAGCaa                                 >  2:416491/1‑147 (MQ=255)                                                                                                                         tGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATG‑CCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCTGTAG‑CCAGCGCCGGCCCCCCCCGCCACCCCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAggga                             >  2:749168/1‑148 (MQ=255)                                                                                                                            cGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                         <  1:413966/151‑1 (MQ=255)                                                                                                                                gaTAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCaa                     <  1:596718/151‑1 (MQ=255)                                                                                                                                      cGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑g              <  1:443968/151‑1 (MQ=255)                                                                                                                                        agaTCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑Gtt            <  2:547193/151‑1 (MQ=255)                                                                                                                                            cggcggGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAg         <  2:473059/150‑1 (MQ=255)                                                                                                                                             ggcggGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAgg        >  1:489118/1‑150 (MQ=255)                                                                                                                                               cggGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAGGATc     <  2:961449/151‑1 (MQ=255)                                                                                                                                                gggggTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAGGATCt    <  2:611960/151‑1 (MQ=255)                                                                                                                                                gggggTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAGGATCt    >  2:541741/1‑151 (MQ=255)                                                                                                                                                gggggTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAGGATCt    <  2:193667/151‑1 (MQ=255)                                                                                                                                                gggggTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCAACCGCGCGAGGCAACGTGGGTATCCCCGTTATTCCCCATCGGTAGGAGCAAGGAATGGGGCAA‑GCCAGGTTCAGGATc     >  1:407193/1‑150 (MQ=255)                                                                                                                                                 ggggTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAGGATCt    >  2:79718/1‑150 (MQ=255)                                                                                                                                                  gggTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAGGATCTaa  <  1:960103/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAG‑CCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCA‑GTTCAGGATCTAA  >  NC_000913/646599‑646894

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