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

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATC  >  NC_000913/646599‑646891                                                                                                                                                     |                                                                                                                                                    aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtt                                                                                                                                                    <  2:747287/150‑1 (MQ=255) aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTc                                                                                                                                                   <  2:296875/151‑1 (MQ=255)      ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                               <  1:421788/150‑1 (MQ=255)         cAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGccc                                                                                                                                           >  1:321745/1‑151 (MQ=255)          aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGcccc                                                                                                                                          <  2:273673/151‑1 (MQ=255)          aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGcccc                                                                                                                                          >  2:371006/1‑151 (MQ=255)           aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGccccc                                                                                                                                         >  2:288618/1‑151 (MQ=255)           aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCATGccccc                                                                                                                                         >  2:362234/1‑151 (MQ=255)           aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGccccc                                                                                                                                         >  2:798819/1‑151 (MQ=255)              ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCctt                                                                                                                                      >  1:67844/1‑151 (MQ=255)               tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCttt                                                                                                                                      <  2:53696/150‑1 (MQ=255)                  cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTg                                                                                                                                   <  2:410880/150‑1 (MQ=255)                     gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCaa                                                                                                                                >  2:367586/1‑150 (MQ=255)                       gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCaataa                                                                                                                             <  1:876097/151‑1 (MQ=255)                       gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGATGCCCCCctttttcta                                                                                                                                >  1:837100/1‑144 (MQ=255)                           aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGCTTGGATGCCCCCCTTTTGTactaaat                                                                                                                           >  2:284214/1‑147 (MQ=255)                            aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTCAATAATGttt                                                                                                                        >  1:322159/1‑151 (MQ=255)                              taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGttt                                                                                                                        >  1:749538/1‑149 (MQ=255)                                  tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGcc                                                                                                                   <  1:868455/150‑1 (MQ=255)                                   aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCTTTTTTCAATAATGTTTTTGCCTc                                                                                                                 >  1:844497/1‑151 (MQ=255)                                    aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAAAGCCCCCCTTTTTCTATAATGTTTtcgtc                                                                                                                   >  1:520000/1‑146 (MQ=255)                                        tgtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCTATAATGTTTTCGCCTccc                                                                                                               >  1:299947/1‑146 (MQ=255)                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCt                                                                                                            <  1:231373/150‑1 (MQ=255)                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAg                                                                                                        <  1:146284/149‑1 (MQ=255)                                               aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAAGCCCCCCTTTTTTTATAATGTTTTTGCCCCCCGCTGGAGcc                                                                                                      >  2:815851/1‑150 (MQ=255)                                                aaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAg                                                                                                    <  2:168509/151‑1 (MQ=255)                                                     aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTAGcc                                                                                                      >  1:464756/1‑144 (MQ=255)                                                        tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAAGCCCCCCTTTTTGAATAATGTTTTGTCCCCCCGCCGGAGCCAgcgc                                                                                                 >  1:278529/1‑146 (MQ=255)                                                         aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑cc                                                                                          <  1:416750/150‑1 (MQ=255)                                                            gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGATGACCCCCTTTTTTAAAAATGTTTTCGCCCCGCGCTGTAGACACCGG‑CGG‑gccccc                                                                                      >  1:365118/1‑151 (MQ=255)                                                             aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGGATGCCCCCCTTTTTCACTAATGTTTGTGCCCCGCGCTGGAGCCAGCGC‑AGG‑ccccccc                                                                                     >  2:429304/1‑150 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTc                                                                                    <  2:120529/150‑1 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCg                                                                                   <  2:636150/151‑1 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCg                                                                                   <  2:321719/151‑1 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGCCATTCCCCCCTTTTTTGATAATTTTTTGTCCCCCCGCCGTTGGCACCGC‑AGG‑ccccccccc                                                                                   >  1:477069/1‑148 (MQ=255)                                                                   aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCa                                                                                <  2:112017/150‑1 (MQ=255)                                                                   aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCAATAATTTTTTCGCCTCCCCCTGTAGCCCGCGC‑AGG‑ccccccccccccc                                                                               >  1:104115/1‑144 (MQ=255)                                                                        aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACcgcgc                                                                          <  2:178842/151‑1 (MQ=255)                                                                          cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATTCCCCCCTTTTTTTAAAATGTTTTCGCCCCCCGCCGTAGCCCGCGC‑AGG‑cccccccccccccccccc                                                                          >  2:274643/1‑137 (MQ=255)                                                                            cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTgcg                                                                                                                     <  2:144320/106‑1 (MQ=255)                                                                             tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGTGGGTTGTATGCCCCCTCTTTGCGATAATGTTTGCGAACCACGCCGTTTCCCGCGCGCGG‑CCCCccccgcccccgccccagg                                                                      >  1:795769/1‑145 (MQ=255)                                                                                 gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGTATGGCCCCCTTTTGTAATAATGTTTGCGCCCCTCGCCGCAGCCAGCGC‑GAG‑CCCCCCCCGCCCCctcgccgagcaa                                                                   >  1:717246/1‑142 (MQ=255)                                                                                  ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGt                                                               <  2:722332/151‑1 (MQ=255)                                                                                     ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGTATGCCCCCCTTTTGCAATAATTTTTTTGCCCCGCCCTGTAGCCAGCGC‑CGG‑GCCCCCCCCCCCCCCCCCGGTG‑CCACGTTg                                                             >  1:184647/1‑150 (MQ=255)                                                                                      gACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGt                                                           <  2:211327/151‑1 (MQ=255)                                                                                            gaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGcc                                                     <  1:22701/151‑1 (MQ=255)                                                                                                 cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTg                                                 <  1:171987/150‑1 (MQ=255)                                                                                                  gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGAt                                               <  1:366386/151‑1 (MQ=255)                                                                                                    tccaTACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGTCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CGCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATc                                              <  2:674642/147‑1 (MQ=255)                                                                                                     cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCg                                             <  2:844430/150‑1 (MQ=255)                                                                                                     cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATGCCCCCCTTTTTCCATAATGTTTTTGCCCCCCGCCGGAGCCCGCCC‑CGG‑GCCCCCCCCCCCCCCCCCGCTGCCACCGTGGGATTCGCCGTTgtcg                                              >  2:55120/1‑147 (MQ=255)                                                                                                      gATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGcc                                           <  1:419196/151‑1 (MQ=255)                                                                                                        tACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCTGTAGCCAGCGC‑AGG‑CCCCCCCCCCCCCCCCCCGATG‑CAACGGTTGGATTGCCGGTGATCGcc                                           >  1:41491/1‑149 (MQ=255)                                                                                                               ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAgg                                  <  1:91125/151‑1 (MQ=255)                                                                                                                        cTGCCGGAGATAATCGAGATCGGCGGGGGGGTGGATGGCCCCCTTTTGCAAAAAATGTTGTGGCCCGCGGGGGAGGCAGCGA‑AGG‑CCCCCCCCCCCCCCCCGCGCAG‑CAAAGCGGGGATCACCGTGGGTCGCCCTCAGAAGGGGGAAggg                          >  2:258169/1‑149 (MQ=255)                                                                                                                            cGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                     <  2:11847/151‑1 (MQ=255)                                                                                                                               agaTAATCGAGATCGGCGGGGGGTTGAATGGCCCCTTTTTTCAATAAATTTTTTCCCCCCCCCCGGAGCCCGCGC‑AGGCCCCCCCCCCCCCCCCCCCGTTG‑CAA‑GGTGGGATCGCCGCTGGTCGCCCTcccgaaggggaaggg                          >  2:200610/1‑129 (MQ=255)                                                                                                                                      cGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCa            <  1:453542/150‑1 (MQ=255)                                                                                                                                      cGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCa            <  2:354256/150‑1 (MQ=255)                                                                                                                                         gagcggcggGGGTTCGAATTCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTATTGCCAGCGC‑AGG‑CTTGTCTCGCCTGGGGGCGATG‑CAGCGTTGGTATCGCCGTTGAGTGCCTGCAGTAGGTGCAAGGTAACGTGCAATGTTAGTTc        <  2:878307/148‑1 (MQ=255)                                                                                                                                            cggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAg      <  1:275692/150‑1 (MQ=255)                                                                                                                                            cggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCCCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCCCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAgg     >  1:108148/1‑151 (MQ=255)                                                                                                                                               cggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc  >  1:700391/1‑151 (MQ=255)                                                                                                                                               cggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc  <  1:367612/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGC‑AGG‑CCCCCCTCGCCACCGCGCGATG‑CAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATC  >  NC_000913/646599‑646891

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