Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,131,440	0	T	C	41.4%	94.9 / 14.5	58	G323G (GGA→GGG)	yghQ	putative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base T (29/5);  new base C (0/24);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.85e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCG  >  NC_000913/3131291‑3131578                                                                                                                                                          |                                                                                                                                           gaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTc                                                                                                                                           >  1:514789/1‑151 (MQ=255) gaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTc                                                                                                                                           <  2:302165/151‑1 (MQ=255)  aaatcaCTGAGTCAGCCGAGAAGAATTTTCCCGCTTATTCGCTCCTTCCCTAAATCAGCTCATACGCTTCGAGATACTTACCGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTcc                                                                                                                                          <  2:379003/151‑1 (MQ=255)    atgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTccgc                                                                                                                                        <  1:815295/151‑1 (MQ=255)      gaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTccgccc                                                                                                                                      >  2:960854/1‑151 (MQ=255)       aCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTccgcccg                                                                                                                                     <  1:891622/151‑1 (MQ=255)         tGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCcgcccgcc                                                                                                                                   >  2:354460/1‑151 (MQ=255)          gAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCcgcccgcc                                                                                                                                   >  2:399187/1‑150 (MQ=255)             tCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGt                                                                                                                                >  2:435503/1‑150 (MQ=255)                    agaagaATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCgg                                                                                                                         >  2:1158066/1‑150 (MQ=255)                         atttcccccGCTTATCCCCCCCTCCCCGAAATCGGGTAATCCGCTGGAAGATCCTAACCGCA‑AAC‑CCCCAGCAAAAGGAGCGGCTTGCCGAAGATCACCCCCGCGGC‑CGCCACA‑AAATTTCCGCTCCCCCCCGCCAGTAAACCGGACTTca                                                                                                                   <  1:403594/146‑1 (MQ=255)                            ttCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCacac                                                                                                                >  2:675598/1‑151 (MQ=255)                              ccccGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCt                                                                                                              >  1:46283/1‑151 (MQ=255)                               ccctttttgtccccccTCCCATAAACCAGGTCATCCGCTCAGAATTACTTACCGCC‑AAC‑CCCCGGAAAAAGGAGGGCCTTGCCGACGACCACCACCGCGGG‑CGCCACC‑AGAACACCAATCCCGCCCGCCAGTAAACCGGACTTCACACCTa                                                                                                             <  1:859089/137‑1 (MQ=255)                               cccGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCt                                                                                                              >  1:32337/1‑150 (MQ=255)                                ccttttttcccccccTCCCCAAAATGGGGTCATACTCTTCGAAATTCTTACCCCC‑AAA‑CCCCAGCAAAATGTGGGGCGTGCCGAGGCTCAACACCGCGAG‑CCCAACA‑AGATTCCCGTTCCGGCCCGCCAGTAAACCGGACTTCACAcc                                                                                                               <  1:1131410/135‑1 (MQ=255)                                   cTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTaaca                                                                                                          >  2:929473/1‑150 (MQ=255)                                    ttatcccccccTCCCCAAAATCGGTCCATACGCTTCGAGATACTAACCGCC‑AAA‑CCCAAGAAAAAGGGGCGGCTGCCCGACGTCCAACCCCGCGAG‑CCCCCCA‑AAAATACCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTaacaa                                                                                                         <  1:860702/141‑1 (MQ=255)                                      aTTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAAcc                                                                                                       >  1:850729/1‑150 (MQ=255)                                       tcccccccTTCCCAAAACCAGCTCATCCGCTTGGAAAT‑TTTAACCCCAAAA‑CCCAAGCAAAAGGGGCGGCTTGCCAACTACAAACCCCGCGCG‑C‑CCACCAAAAATCCCGTTCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCa                                                                                                      <  1:806914/144‑1 (MQ=255)                                         cgcccctcccaaaaaTCAGGTAATCCGCTTCGAGATACTAACCGCA‑AAA‑CCCCAGAAAAAGGGGCGGCTTGCCGACAACCACCCCCGCGCG‑CGCCCCA‑AAAATCCCGATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACgg                                                                                                   <  1:240697/139‑1 (MQ=255)                                         cgcccctcccaaaaaCCAGGTAATACCCTTCGAAATCCTAACCCCA‑AAA‑CCCAAGCAAAAGGGGCGGCTGGCCGACGACAACCCCCGCGAG‑CCCCCCA‑AAAAACCCGTTCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACg                                                                                                    <  1:77137/138‑1 (MQ=255)                                           cACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTc                                                                                                 >  1:721293/1‑151 (MQ=255)                                            aCCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCt                                                                                                >  1:981552/1‑151 (MQ=255)                                             ccccccaaaaaCCAGGTCATCCGCTTAGAAATACTAAACACC‑AAC‑CCCCGGGGAAATGGGCGGCTCACCGCCAATCACCCCCGCGAG‑CCCCCCA‑AAAACACCGATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCtg                                                                                               <  2:211897/143‑1 (MQ=255)                                             ccTCCCCAAAAACAGGTAATGCTCTTAGAGATACTAAACCCA‑AAA‑CCCCAGCGAAATGAGCGGCTTCCCAACGACCACCCCCGCGGG‑CCCCACC‑AAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCt                                                                                                <  1:315177/150‑1 (MQ=255)                                              cTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGCAAACCGGACTTCACACCTAACAACCACGGTCtgg                                                                                              >  1:1093783/1‑151 (MQ=255)                                                     ataTCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGCATACCGATTCCGCCCGCCAGTCAACCGGCCTTCACACCTAACAACCACGGTCTGGTGGTtcg                                                                                       >  2:504894/3‑151 (MQ=255)                                                         tagtccccccGCTTCGACACAATCAACACA‑AAA‑CACCAGAAAATGGGGCGGCTTCCCGAC‑AACTACCCCGCGCGCCCCCACC‑AAATCCCCGCGCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCgg                                                                                    <  1:808500/141‑1 (MQ=255)                                                          aGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGAt                                                                                  >  2:603564/1‑151 (MQ=255)                                                               accccctgggaaaTACTAAACGCA‑AAC‑CCCCAGAGAAGGGAGCGGTTGCCCGGAAAACAACCCCGCGCG‑CCCCACA‑AGAACACCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGa                                                                                   <  1:328864/133‑1 (MQ=255)                                                                tACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAAc                                                                            >  1:538175/1‑151 (MQ=255)                                                                     ttCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATc                                                                       >  2:129728/1‑151 (MQ=255)                                                                       ggagaTATTCAACGAC‑AAC‑CACCAATAAAGTGGGCGGCTTCCCCACGCCCAACCCGGCGAG‑CCCACCC‑AGAAGGCCTTTTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCAc                                                                     <  2:1069544/150‑1 (MQ=255)                                                                       cGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCAc                                                                     >  2:555240/1‑151 (MQ=255)                                                                        gagaTACTTAACGCC‑AAA‑CAGCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGCTGCGCGGATCTAAACGCATCAcc                                                                    >  2:26445/1‑151 (MQ=255)                                                                         gagttactAACCCC‑AAACCCCCAGCAAAAGGGGCGGCTTGCCGCCCACCACCCCCGCGGG‑CCCCCCA‑AAAACCCTCATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCAc                                                                     <  1:675807/143‑1 (MQ=255)                                                                         agaTACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACct                                                                   >  1:744532/1‑151 (MQ=255)                                                                           aTACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCACCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACctct                                                                 >  2:1059040/1‑151 (MQ=255)                                                                            tACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTg                                                                >  2:894326/1‑151 (MQ=255)                                                                             actaaccgac‑AAA‑CACCAGCAAAACGAGCGGCTTCCCGACAATCAACCCCGCGCG‑CCCCACC‑AAAATACCGTTCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACctct                                                                 <  1:355086/140‑1 (MQ=255)                                                                               ttAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGt                                                             >  2:16418/1‑151 (MQ=255)                                                                                 accGCC‑AAA‑CCCCAGAGAAATGGGCGGCTTGCCGACGATCACCCCCGCGAG‑CGCCACC‑AAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGt                                                             <  1:283159/147‑1 (MQ=255)                                                                                     cc‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAAc                                                       >  1:996452/1‑151 (MQ=255)                                                                                        aaa‑CCCCAGAAAAAGGGGCGGCTTCCCGACGACCACCCCCGCGAG‑CCCCACA‑AAAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAAc                                                       <  1:156471/149‑1 (MQ=255)                                                                                                        gggcgcccccccaacgaTCACCCCCGCGAG‑CCCCACA‑AAAACCCCGTCCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAAcc                                       <  2:168459/137‑1 (MQ=255)                                                                                                            ggCTTGCCGACAATCAACCCCGCGCG‑CCCCACC‑AGACTACCGATCCCCCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTg                                  <  1:1156653/151‑1 (MQ=255)                                                                                                               tcgCCGACAAACAACCCCGCGAG‑CCCCCCC‑AAAACCCCGTTCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTg                                  <  1:444141/146‑1 (MQ=255)                                                                                                                     acgaTCAACCCCGCGAC‑CCCCCCC‑AGAAAACCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCt                                   <  2:939087/141‑1 (MQ=255)                                                                                                                     acgaTCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcgg                         >  1:216201/1‑151 (MQ=255)                                                                                                                      caaccacccccGCGAG‑CCCCACA‑AGAATCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggc                        <  1:435518/142‑1 (MQ=255)                                                                                                                       gaTCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggcg                       >  1:520438/1‑151 (MQ=255)                                                                                                                               ccGCGAG‑CCCCACA‑AGAATCCCCACCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTc                  <  2:197289/148‑1 (MQ=255)                                                                                                                               ccGCGAC‑CGCCACA‑AAAATCCCGCCCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcggcg                       <  2:983444/143‑1 (MQ=255)                                                                                                                                   gAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAAt           >  2:599634/1‑151 (MQ=255)                                                                                                                                    gacccccacc‑aAAACCCCGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaag              <  2:865187/143‑1 (MQ=255)                                                                                                                                               aaaaccccGATCCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGc   <  1:504592/143‑1 (MQ=255)                                                                                                                                               aGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCg  >  1:14384/1‑150 (MQ=255)                                                                                                                                                          |                                                                                                                                           GAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCC‑AAA‑CACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAG‑CGCCACC‑AGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCG  >  NC_000913/3131291‑3131578

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