Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,765	0	G	T	42.4%	77.3 / 37.4	60	Q249K (CAA→AAA)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base G (4/30);  new base T (25/0);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 7.84e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAG  >  NC_000913/646616‑646906                                                                                                                                                      |                                                                                                                                              cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTg                                                                                                                                               <  1:105956/150‑1 (MQ=255) cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGc                                                                                                                                              <  2:125613/151‑1 (MQ=255)  aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCa                                                                                                                                             <  2:998936/151‑1 (MQ=255)        ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCtataa                                                                                                                                         >  1:349117/1‑149 (MQ=255)            ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGATGGCCCCCTTTTTCAATAAtttt                                                                                                                                     >  2:88364/1‑149 (MQ=255)                ttAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTCAATAATGTTTTCGcc                                                                                                                               >  1:396313/1‑151 (MQ=255)                    aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTTGCCCcgc                                                                                                                           >  1:303911/1‑151 (MQ=255)                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTg                                                                                                                       <  1:695699/151‑1 (MQ=255)                        gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGGAGTCCCCTTTTTGTTATAATATTTTCGCCTcg                                                                                                                            >  1:517700/1‑146 (MQ=255)                         tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTGCCATAATGtttttgccccccgc                                                                                                                         >  2:720054/1‑137 (MQ=255)                          gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGt                                                                                                                      <  2:301216/150‑1 (MQ=255)                          gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATTTTTTTGCCTCGCGCTg                                                                                                                       >  1:979082/1‑149 (MQ=255)                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                                  <  2:91540/151‑1 (MQ=255)                                aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgc                                                                                                               <  2:484594/151‑1 (MQ=255)                                     ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAAAGTTTTCGCCCCGCGCCGTAGCCAgcg                                                                                                              >  1:946920/1‑147 (MQ=255)                                     ggTAAGGATCAGTAGGTCAGCACTGCCGACTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTTGCTGTCCCCGTTTTTTAAAAACGTTTTCTCGGCCCCCTGTAGCCAGCCCAg                                                                                                           >  1:558427/1‑150 (MQ=255)                                      gTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGc                                                                                                         <  2:10037/151‑1 (MQ=255)                                           gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCTTTTTTGAATAATTTTTTCGCCCCGCGCTGGAGCCAGCGCAGGgccccc                                                                                                    >  1:139090/1‑151 (MQ=255)                                            aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCAAGCCCCCCTTTTTCCATAATGTTTTCGCCCCCCGCCGTAGCCAGCGCCGGgcccccc                                                                                                   >  1:127110/1‑150 (MQ=255)                                            aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATGCCCCCCTTTTTTAATAAAGTTTTTGCCCCCCCCGGTAGCCCCCGCAGGccccccc                                                                                                   >  2:831836/1‑150 (MQ=255)                                              cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                                  <  2:927407/150‑1 (MQ=255)                                              cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATTCCCCCCTTTTTCTATAAATTTTTCGCCCCCCGCCTTAGCCAGCGCAGGcccccccc                                                                                                  >  1:934759/1‑148 (MQ=255)                                                 tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAATTTTTCGCCCCGCGCGGTTGCCAGCGCCGGcccccccccccc                                                                                              >  2:848115/1‑145 (MQ=255)                                                    gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcg                                                                                           <  1:177926/151‑1 (MQ=255)                                                    gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcg                                                                                           <  1:365629/151‑1 (MQ=255)                                                     tCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCTTTTTTCAATAAAATTTTTCCCTCGCCCTGTGGGCAACGCAGGcccccccccccccccc                                                                                          >  1:607568/1‑141 (MQ=255)                                                      cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTTAATAAAATTTTTGCCCCCCCCTGTTGCCCCCGCACGcccccccccccccccc                                                                                          >  2:831223/1‑140 (MQ=255)                                                        gCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGCTGCCCCCCTTTTTGAATAATGTTTTTCCCCCGCGGAGTTGCAAGCGCCGGCCCCCCCCGCcccc                                                                                            >  2:263864/1‑146 (MQ=255)                                                         cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcgcg                                                                                       <  1:125719/150‑1 (MQ=255)                                                           cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCGATAAAGTTTTCGCCTCGCCCTGTAGCCCGCGCAGGCCCCCCCCGCCACCCCGCGGTg                                                                                    >  1:104593/1‑151 (MQ=255)                                                           cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTCAATAAATTTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCCCGCCCCCCCGCGGTg                                                                                    >  1:613157/1‑151 (MQ=255)                                                               cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATTACCCCCTTTTTTAATAATTTTTTTGCCCCGCGCCGTAGGCAGCCCCGGCCCCCCTCcccccccccc                                                                                        >  2:729287/1‑139 (MQ=255)                                                                gccTGGACTGAGATTTCGTTCGATACTCTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACg                                                                               <  2:418777/151‑1 (MQ=255)                                                                   tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAGTGCCCCCTTTTTTCAATGATGTTTGTCGGGCGCGCTGTAGGCGGCGCCCGCCCCCCCCCCCCCCCCGCGGGGCAACggt                                                                             >  2:429498/1‑148 (MQ=255)                                                                         gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                                       <  2:793565/150‑1 (MQ=255)                                                                         gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCg                                                                      <  1:21850/151‑1 (MQ=255)                                                                              ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCAAAGTTGGTATCGcc                                                                    >  1:80131/1‑148 (MQ=255)                                                                                cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                                <  2:853609/150‑1 (MQ=255)                                                                                  ttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAATTTTTTTGCCTCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGGGCAACGGTGGTATc                                                                       >  2:577743/1‑141 (MQ=255)                                                                                   tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTTAATAATATTTTTGCCCCGCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGATGGAACGGTTGGATTGCCGGTg                                                                >  1:491061/1‑147 (MQ=255)                                                                                   tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCCATAATTTTTTCGCCCCGCGCGGTAGGCAGCGCAGGCCCCCCCCGCCACCCCGCGAAGCAAAGGTGGGATTGCCGGTggt                                                              >  1:784547/1‑147 (MQ=255)                                                                                       ttcgcgcccGGTCGAACGGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCGCCGCGCGATGCAACGTTGGTGTCGCCGTTGATCg                                                            <  2:183612/141‑1 (MQ=255)                                                                                       tACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCGATAATATTTTTGCCCCGCGCCGGAGCCCGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCAACGCTGGGGTCGCCCTTGTTCGc                                                           >  2:9329/1‑148 (MQ=255)                                                                                           cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAg                                                     <  2:261861/150‑1 (MQ=255)                                                                                            tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCTGTAGCCAGCGCAGGGCCCCCTCGCCACCCCGCGCGGCGACGTTGGGATCGCCGTTGAtag                                                            >  1:865459/1‑140 (MQ=255)                                                                                              ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAgg                                                 <  2:270840/151‑1 (MQ=255)                                                                                                 gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                               <  2:740859/150‑1 (MQ=255)                                                                                                      aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGt                                         <  1:332179/151‑1 (MQ=255)                                                                                                             gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                                  <  2:968605/151‑1 (MQ=255)                                                                                                                aTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAt                               <  2:127099/151‑1 (MQ=255)                                                                                                                           cggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAgg                    <  2:330954/151‑1 (MQ=255)                                                                                                                            ggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGa                   <  2:214427/151‑1 (MQ=255)                                                                                                                            ggcgaGGGTTCGGATGACCCCTTTTTCCAATAATGGTTGCGCCTCGCGCTGCAGCCAGCGGCGGCCCCCCACGCCACCGCGCGATGCAACGCTTGTATCCCCGATGATCGCCATCAGTAGGAGCTAGGTATCGAGCCATGCCAGt                         >  2:790351/1‑145 (MQ=255)                                                                                                                             gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt                  <  2:772332/151‑1 (MQ=255)                                                                                                                              cggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc                 <  1:88378/151‑1 (MQ=255)                                                                                                                              cggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc                 <  2:971583/151‑1 (MQ=255)                                                                                                                                  ggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAc             >  1:69254/1‑151 (MQ=255)                                                                                                                                         aTGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGAt      >  1:749951/1‑151 (MQ=255)                                                                                                                                            cccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCca   <  1:218443/151‑1 (MQ=255)                                                                                                                                            cccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCca   <  2:45233/151‑1 (MQ=255)                                                                                                                                             ccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcag  <  1:197526/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                              CAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAG  >  NC_000913/646616‑646906

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