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

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646599‑646896                                                                                                                                                      |                                                                                                                                                      aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTc                                                                                                                                                     <  1:258390/151‑1 (MQ=255)  atatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCg                                                                                                                                                    <  2:184137/151‑1 (MQ=255)  atatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCg                                                                                                                                                    <  1:125085/151‑1 (MQ=255)  atatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCg                                                                                                                                                    <  1:25283/151‑1 (MQ=255)  atatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCg                                                                                                                                                    <  2:507291/151‑1 (MQ=255)    atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGaa                                                                                                                                                  <  1:373225/151‑1 (MQ=255)      ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAAt                                                                                                                                                 <  1:397205/150‑1 (MQ=255)      ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGaa                                                                                                                                                  >  1:129747/1‑149 (MQ=255)           aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGGATTccccc                                                                                                                                           >  2:491408/1‑151 (MQ=255)            aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCt                                                                                                                                          <  2:773087/151‑1 (MQ=255)            aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGATGccccc                                                                                                                                           >  2:385994/1‑150 (MQ=255)               tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCttt                                                                                                                                        <  2:573298/150‑1 (MQ=255)               tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGAATGGCCCCcttt                                                                                                                                       >  2:105230/1‑151 (MQ=255)                   aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATGCCCCCttttt                                                                                                                                      >  2:9162/1‑148 (MQ=255)                    aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCaa                                                                                                                                  <  1:485112/151‑1 (MQ=255)                      tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGTGTTGGATTCCCCCttttt                                                                                                                                      >  1:780033/1‑145 (MQ=255)                        cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTGGTAAGCCCCCCTTTTTCaataa                                                                                                                               >  1:721030/1‑150 (MQ=255)                           aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATGCCCCCcttttttaaaaa                                                                                                                               >  1:524394/1‑140 (MQ=255)                            aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGttt                                                                                                                          <  2:352944/151‑1 (MQ=255)                               aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                        <  2:460158/150‑1 (MQ=255)                               aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                        <  2:191344/150‑1 (MQ=255)                               aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgc                                                                                                                        <  1:636718/150‑1 (MQ=255)                               aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGATGCCCCCCTTTTTTAATAATtttt                                                                                                                          >  2:856753/1‑148 (MQ=255)                                atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCCTTTTTCCATAATGtttt                                                                                                                         >  2:1142929/1‑147 (MQ=255)                                  tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGATGGCCCCCTTTTTTAATAAAGTTTGTGcc                                                                                                                     >  1:130519/1‑150 (MQ=255)                                   aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTc                                                                                                                   <  2:334211/151‑1 (MQ=255)                                   aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGAAGCCCCCCTTTTTTAATAAAGTTTTCGcccc                                                                                                                   >  2:434636/1‑149 (MQ=255)                                       cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGAATGGCCCCTTTTTTGAATAATTTTTTCGCCCcgc                                                                                                                 >  2:44750/1‑149 (MQ=255)                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAAAGTCCCCCTTTTTCAATAATGTTTTCGCCCCccgc                                                                                                               >  2:849361/1‑149 (MQ=255)                                          tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGAATGCCCCCCCTTTTTAAAAATATTTTCGCCCcgc                                                                                                                 >  1:1158267/1‑146 (MQ=255)                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGAATGCCCCCCTTTTTCAATAATATTTTTGCCTCCCGCCGTAg                                                                                                          >  2:862705/1‑149 (MQ=255)                                               aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATGCCCCCCTTTTTTAATAATGTTTTCGCCCCGCGCCGTAggc                                                                                                        >  2:1157099/1‑148 (MQ=255)                                                 aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgc                                                                                                     <  1:463518/151‑1 (MQ=255)                                                    cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATTCCCCCCTTTTGCAATAATGTTTTTGCCCCCCGCTggtgcc                                                                                                        >  2:176485/1‑140 (MQ=255)                                                        tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGATGGCCCCCTTTTTCTATAAAGTTTGGGCCTCGCGCCGGTGCCCGCGCAGGcc                                                                                              >  1:33440/1‑151 (MQ=255)                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCAAGGCCCCCTTTTTG‑ACAAAAGTTTTTTCCTCGCGCCGTGGGCCGCGCACGccc                                                                                             >  1:611029/1‑151 (MQ=255)                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                                              <  1:825209/150‑1 (MQ=255)                                                                  tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATGCCCCCCCTTTTCCATAAAGTTTTCGCCTCCCGCCGTAGCCCGCGCCGGcccccccccccc                                                                                    >  2:652760/1‑145 (MQ=255)                                                                   aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGATGCCCCCTTTTTTCAATAATGTTTTTGCCCCGCGCTGTAGCCCGCGCAGGCCCCCCCCGcccc                                                                                   >  2:500884/1‑149 (MQ=255)                                                                   aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACCCCCTGTCGAACTGGCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTATGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGGCCGCCTCGcct                                                                                    <  1:1169077/150‑2 (MQ=255)                                                                      tCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTGGGATGCCCCCTTTTTTTAAAAAAGTTTTTGGCTCCCGCCGTGGCCAGCGCCGGCCCCCCTccgccccc                                                                                 >  2:820037/1‑143 (MQ=255)                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATGCCCCCCTTTTTCAAAAATGTTTTCGCCCCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCCCcccgc                                                                              >  2:1103854/1‑151 (MQ=255)                                                                            cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATGCCCCCCTTTTTTAAAAATGTTTTTGCCCCCCGCTGTAGCCAGCGCCGGCCCCCCCCCCCACCCCGCGATg                                                                          >  2:691056/1‑151 (MQ=255)                                                                              gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCa                                                                        <  1:408666/151‑1 (MQ=255)                                                                                     ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTgg                                                                 <  1:310646/151‑1 (MQ=255)                                                                                     ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTg                                                                  <  2:1166540/150‑1 (MQ=255)                                                                                     ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTg                                                                  <  1:51259/150‑1 (MQ=255)                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                             <  1:1148807/150‑1 (MQ=255)                                                                                               ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGATGCCCCCCTTTTGCAATAAAGTTTTCGCCTCGCGCCGTAGGCCGCGCAGGCCCCCCCCCCCCCCCCCCGGTGCAACGGTGGGATCGcc                                                          >  1:289741/1‑148 (MQ=255)                                                                                                tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                       <  1:585867/150‑1 (MQ=255)                                                                                                tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                       <  2:496946/150‑1 (MQ=255)                                                                                                tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATGCCCCCTTTTTTTAATAATGTTTTCGCCTCCCGCGGTAACCAGCGCCCGCCCCCCCCCCCCCCCCGCGGTGCAACGTTGGGAACGcc                                                          >  1:364662/1‑147 (MQ=255)                                                                                                          cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                            <  2:67218/151‑1 (MQ=255)                                                                                                           acTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATGCCCCCCTTTTTTCATAATTTTTTCGCCCCCCCCCGTAGCCAGCGCCAGCCCCCCCCCCCCCCCCGCCGTGCAACGGTGGTATCGCCCGTGATCTccccc                                             >  1:369496/1‑146 (MQ=255)                                                                                                             tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAAATCCCCCCTTTTTCAATAAATTTTTCCCCTCGCCCCGGAGCCAGCCCAGGCCCCCCCCCCCCCCCCGCGGTGCAACggttggatgccc                                                          >  1:647951/1‑123 (MQ=255)                                                                                                              cccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTa                                         <  1:1031686/150‑1 (MQ=255)                                                                                                                    cGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCaa                                  <  1:65224/151‑1 (MQ=255)                                                                                                                     gAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGAATGCCCCCTTTTTTCAATAAAGTTTTCGCCTCCCCCTGTAGCCACCGCAGGCCCCCCCCCCCCCCCCCCGGTGCCACGTTGGTATCGCCGGTGATAGCCAACAGTAGGAg                                     >  2:310801/1‑147 (MQ=255)                                                                                                                      aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAgg                                <  2:721045/151‑1 (MQ=255)                                                                                                                          gCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATc                            <  2:922619/151‑1 (MQ=255)                                                                                                                                 aTAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAt                     <  2:1061150/151‑1 (MQ=255)                                                                                                                                    aTCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGcc                  <  1:269090/151‑1 (MQ=255)                                                                                                                                       gagaTCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGt               <  1:326382/151‑1 (MQ=255)                                                                                                                                        agaTCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt              <  1:187432/151‑1 (MQ=255)                                                                                                                                              gc‑ggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt        <  2:809494/151‑1 (MQ=255)                                                                                                                                                 gggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      >  1:593937/1‑151 (MQ=255)                                                                                                                                                  ggggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa     <  2:1099515/151‑1 (MQ=255)                                                                                                                                                     gTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc  >  2:1157785/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646599‑646896

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