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

TCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAA  >  NC_000913/646605‑646894                                                                                                                                                |                                                                                                                                                          tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGATATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑g                                                                                                                                                 <  2:851709/150‑1 (MQ=255) tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAt                                                                                                                                                   <  2:1035600/149‑1 (MQ=255)    aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑cccc                                                                                                                                            <  1:732506/151‑1 (MQ=255)      aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCt                                                                                                                                          <  1:355631/151‑1 (MQ=255)          aTCAAGTGCTTAAATAATTAAATCTGCGCTAACAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACCCTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGTGGTTC‑G‑AAC‑G‑GCCCCcttt                                                                                                                                       >  2:816673/1‑150 (MQ=255)            cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTG‑G‑CAT‑G‑ACCCCctttt                                                                                                                                      >  2:403812/1‑149 (MQ=255)             aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCa                                                                                                                                   <  1:458800/151‑1 (MQ=255)                  cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑AAT‑G‑CCCCCCTTTTTCaataa                                                                                                                               >  1:90517/1‑150 (MQ=255)                   ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑GAT‑G‑CCCCCTTTTTTCaaaa                                                                                                                                >  2:9098/1‑146 (MQ=255)                    tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑AAT‑G‑CCCCCcttttttta                                                                                                                                  >  1:853489/1‑141 (MQ=255)                     aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGt                                                                                                                            <  2:370076/150‑1 (MQ=255)                      aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGttt                                                                                                                          <  2:706737/151‑1 (MQ=255)                       ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGC‑G‑AAT‑G‑CCCCCCTTTTTCcataaatt                                                                                                                            >  2:710939/1‑145 (MQ=255)                         aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑GAT‑G‑CCCCCCTTTTTCTATAATGttt                                                                                                                          >  1:11938/1‑148 (MQ=255)                           ttAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑C‑AAT‑T‑CCCCCCCTTTTCCATAATATTTTCGcc                                                                                                                     >  1:891575/1‑151 (MQ=255)                            tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑AAT‑G‑TTCCCCTTTTTTAATAAAATTTtcgc                                                                                                                      >  2:1042328/1‑149 (MQ=255)                             aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCt                                                                                                                    <  1:168477/150‑1 (MQ=255)                              aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑CGT‑G‑GCCCCCCTTTTCAAAAAAGTTTGGGcccc                                                                                                                   >  2:975032/1‑148 (MQ=255)                                  tgtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTG‑C‑AAT‑G‑CCCCCctttttctaaaaaattt                                                                                                                          >  2:94383/1‑127 (MQ=255)                                    tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑AAT‑G‑CCCCCCTTTTTTAATAATATTTTTGCCTccc                                                                                                                 >  1:56614/1‑144 (MQ=255)                                     gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑AAT‑G‑CCCCCCTTTTTTAATAATttttttgccccccgc                                                                                                               >  1:111528/1‑136 (MQ=255)                                       tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑AAA‑G‑CCCCCCTTTTTCCATAATGTTTTCGCCCCG‑CGCTGTAg                                                                                                         >  1:606800/1‑150 (MQ=255)                                        aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAg                                                                                                         <  1:470740/149‑1 (MQ=255)                                        aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGC‑G‑CGT‑GCCCCCCTTTTTGCAA‑AATATTTGTGccccccccggttg                                                                                                          >  2:114333/1‑138 (MQ=255)                                         aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCa                                                                                                      <  2:1172082/151‑1 (MQ=255)                                         aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCa                                                                                                      <  1:640015/151‑1 (MQ=255)                                           aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑GAT‑G‑GCCCCCTTTTTCCATAATATTTTTGCCCCG‑CGCCGTAGcccg                                                                                                     >  2:152582/1‑148 (MQ=255)                                             ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑GAT‑G‑GCCCCCTTTTTCAATAAAATTTTTGCCTCG‑CGCTGTGGCCAgcg                                                                                                   >  2:8494/1‑150 (MQ=255)                                               aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGG‑G‑ATT‑G‑CCCCCTTTTT‑TAATAATAtttttggcccgcgccgt                                                                                                            >  2:497026/1‑126 (MQ=255)                                                ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑gg                                                                                              <  1:493710/151‑1 (MQ=255)                                                   aaGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGcc                                                                                            <  1:407217/150‑1 (MQ=255)                                                    aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑GAA‑G‑CCCCCCTTTTTTAATAATATTATCGCCTCC‑CGCCGTAGCCACCGCA‑GGgccc                                                                                          >  1:662690/1‑151 (MQ=255)                                                      gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCCCTTTTTAATAAAATTTTTGCCCCG‑CGCCGTTGCCAGCGCC‑GGcccccc                                                                                        >  2:1102126/1‑151 (MQ=255)                                                         cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGCC‑G‑GAT‑G‑CCCCCCTTTTTCAAGAATGTTTTCGCCCCG‑CGCTGTGGCCCGCCCACGGcccccccc                                                                                      >  2:242877/1‑149 (MQ=255)                                                           gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑CGT‑GACCCCCTTTTTGCAAAAATTTTTCCGCC‑CG‑CGCGGGTGCCCCCGCA‑GGccccccccccc                                                                                   >  2:20407/1‑146 (MQ=255)                                                              ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCAcc                                                                                <  2:56613/151‑1 (MQ=255)                                                                tCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑T‑AAT‑G‑TCCCCCTTTTGTAAAAAATATTTCGCCTCC‑CCCGGGAGCCCGCGCA‑GGgcccccccctcaccac                                                                              >  1:406960/1‑141 (MQ=255)                                                                   gCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑GAT‑G‑CCCCCCTTTTTTTATAATATTTTCGCCTCC‑CGCCGTAGCCCGCGCA‑GGCCCCCCCCCCCCCcccgc                                                                            >  1:917866/1‑150 (MQ=255)                                                                     aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTT‑G‑AAA‑G‑CCCCCCTTTTTTAATAATGTTTTCGCCCCG‑CGCTGGGGCCCGCGCA‑GGGCCCCCCCCCCCCcccgcg                                                                           >  1:188270/1‑149 (MQ=255)                                                                     aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑G‑GATGG‑CCCCCTTTTTG‑AAAAATATTTTTGCCCCC‑CCCGCTAGGCAGCCCA‑AGCCCCCCCCGCccccc                                                                               >  1:475943/1‑144 (MQ=255)                                                                        gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑G‑AAT‑G‑CCCCCCTTTTTCAATAATGTTTTCGCCTCC‑CGCCGTAGCCAGCGCA‑GGCCCCCCCCCCCCCCCCCCGATg                                                                        >  1:882654/1‑149 (MQ=255)                                                                           gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACg                                                                   <  2:551206/151‑1 (MQ=255)                                                                            ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGG‑GCAAT‑G‑CCCCCTTTTTCCGA‑AATTTTTGTGCCTCA‑CGCTGTAGCCAGGCCC‑GGCCCCCCCCCCCACCCCCCGGTCCAAcac                                                                  >  2:744247/1‑149 (MQ=255)                                                                               ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTgg                                                               <  1:1064504/151‑1 (MQ=255)                                                                                           cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                    <  1:474023/150‑1 (MQ=255)                                                                                            gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑AAT‑C‑CCCCCTTTTT‑CAATAATGTTTTCGCCCCC‑CGCTGTAGCCCGCGCC‑GGCCCCCCCCCCCCCCCCCCGATGCAACGTTGGTATCGCCCGTGAt                                                  >  1:955017/1‑151 (MQ=255)                                                                                              tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                                 <  2:11938/150‑1 (MQ=255)                                                                                                        cccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTa                                       <  1:83938/150‑1 (MQ=255)                                                                                                          cTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑C‑AAT‑T‑CCCCCCTTTTTTTATAATGTTTTCGCCCCGCCGCTTTGGCC‑GCGGA‑GGCCCCCCCCCCCCCCCCCCGCAGCCACGCTGGGATTGCCGCTGATCGCCCTCCGTAg                                      >  1:708197/1‑149 (MQ=255)                                                                                                                aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAgg                              <  1:1079006/151‑1 (MQ=255)                                                                                                                    gCCGGAGATAATCGAGATCGGC‑GTGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCC‑GGCTCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATc                          <  2:743698/151‑1 (MQ=255)                                                                                                                    gCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTa                            <  1:952486/149‑1 (MQ=255)                                                                                                                        gagaTAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                      <  1:121285/151‑1 (MQ=255)                                                                                                                           aTAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAt                   <  1:352497/151‑1 (MQ=255)                                                                                                                                    aTCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa          <  2:946030/151‑1 (MQ=255)                                                                                                                                       ggc‑ggGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGa       <  1:170614/151‑1 (MQ=255)                                                                                                                                         c‑ggGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt      >  2:1024040/1‑150 (MQ=255)                                                                                                                                            ggggTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTCGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt    >  2:798404/1‑150 (MQ=255)                                                                                                                                             gggTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTaa  <  1:99535/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                          TCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑G‑AAT‑G‑CCCCCTTTTTGCAATAATGTTTGCGCCTCG‑CGCTGTAGCCAGCGCA‑GGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAA  >  NC_000913/646605‑646894

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