Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,393	C→G	A159G (GCT→GGT)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,393	0	C	G	100.0%	258.6 / NA	75	A159G (GCT→GGT)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (31/44);  total (31/44)

GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GC‑TGGCG‑‑‑‑A‑‑‑‑ACT‑GG‑CTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458                                                                                       |                                                                                 gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a                                                                       <  2:1404079/93‑1 (MQ=255)  cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a‑‑‑‑a                                                                  <  1:477398/93‑1 (MQ=255)  cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a‑‑‑‑a                                                                  <  2:5007933/93‑1 (MQ=255)  cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a‑‑‑‑a                                                                  <  1:2722700/93‑1 (MQ=255)   gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a‑‑‑‑ag                                                                 >  1:1750405/1‑92 (MQ=255)    aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGTCATGCGCTTGAAG‑GCGTGGCgcaga                                                                        >  1:2318867/1‑89 (MQ=255)     tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a‑‑‑‑agaa                                                               >  2:3486276/1‑90 (MQ=255)      cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a‑‑‑‑agaaa                                                              >  1:1260781/1‑89 (MQ=255)                                      tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtt                                 >  1:1148257/1‑93 (MQ=255)                                      tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:5720128/93‑1 (MQ=255)                                      tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:3052296/93‑1 (MQ=255)                                      tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:2423527/93‑1 (MQ=255)                                        tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               <  2:4945205/93‑3 (MQ=255)                                        tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               <  2:2720365/93‑3 (MQ=255)                                        tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               <  2:2612051/93‑3 (MQ=255)                                        tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:5203346/93‑3 (MQ=255)                                        tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:6351818/93‑3 (MQ=255)                                        tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               >  1:5279004/1‑91 (MQ=255)                                         ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttgg                              >  2:2368692/1‑90 (MQ=255)                                         ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttgg                              <  1:1390840/93‑4 (MQ=255)                                          gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑AAGAGACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttggg                             >  2:3799169/1‑89 (MQ=255)                                          gAAGTGATGGCGGGCGGCGTCTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:2138944/93‑5 (MQ=255)                                                aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  2:3643840/93‑1 (MQ=255)                                                aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       >  1:4895227/1‑93 (MQ=255)                                                aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  1:4539782/93‑1 (MQ=255)                                                 tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  2:2652227/93‑1 (MQ=255)                                                  ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:1612030/93‑1 (MQ=255)                                                  ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:1677885/93‑1 (MQ=255)                                                  ggcgcgcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACT‑GGCCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:819012/93‑1 (MQ=255)                                                   gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg                    <  2:4096718/93‑1 (MQ=255)                                                    cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  1:3726700/93‑1 (MQ=255)                                                     gggcggCGTTTCAAAAGGGCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  2:3318022/93‑1 (MQ=255)                                                     gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  1:3501601/93‑1 (MQ=255)                                                     gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  >  2:201803/1‑93 (MQ=255)                                                      ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  1:4668012/93‑1 (MQ=255)                                                      ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  2:4169501/1‑93 (MQ=255)                                                      ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  2:2472633/93‑1 (MQ=255)                                                      ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  2:2412012/1‑93 (MQ=255)                                                      ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  2:1145412/1‑93 (MQ=255)                                                       gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                <  2:733206/93‑1 (MQ=255)                                                       gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                <  1:3946254/93‑1 (MQ=255)                                                        cggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:5029312/93‑1 (MQ=255)                                                        cggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:3061014/93‑1 (MQ=255)                                                        cggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:3733356/93‑1 (MQ=255)                                                         ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  1:838731/1‑93 (MQ=255)                                                         ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  2:6002290/93‑1 (MQ=255)                                                         ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  2:506903/1‑93 (MQ=255)                                                         ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  2:6349720/93‑1 (MQ=255)                                                         ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  1:1387788/93‑1 (MQ=255)                                                         ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  1:2279083/93‑1 (MQ=255)                                                          gcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  1:1938896/1‑93 (MQ=255)                                                          gcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             <  1:1704662/93‑1 (MQ=255)                                                          gcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  1:550341/1‑93 (MQ=255)                                                           cgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  2:2284238/93‑1 (MQ=255)                                                            gTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:1050943/1‑93 (MQ=255)                                                            gTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:5795906/1‑93 (MQ=255)                                                            gTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  2:3718784/93‑1 (MQ=255)                                                            gTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:3054417/1‑93 (MQ=255)                                                                caaaaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:2486652/91‑1 (MQ=255)                                                                 aaaaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:269946/2‑93 (MQ=255)                                                                 aaaaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:2970238/2‑93 (MQ=255)                                                                  aaaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  2:1611781/1‑93 (MQ=255)                                                                   aaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  2:4329317/93‑1 (MQ=255)                                                                   aaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  1:501772/1‑93 (MQ=255)                                                                   aaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  1:1121834/1‑93 (MQ=255)                                                                   aaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  1:5171996/1‑93 (MQ=255)                                                                    aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:6392472/1‑93 (MQ=255)                                                                    aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:871705/93‑1 (MQ=255)                                                                    aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:3051662/1‑93 (MQ=255)                                                                    aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:2466575/1‑93 (MQ=255)                                                                    aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:4138537/1‑93 (MQ=255)                                                                    aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:616170/1‑93 (MQ=255)                                                                    aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:4308666/1‑93 (MQ=255)                                                                    aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:5530256/1‑93 (MQ=255)                                                                    aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:4406304/1‑93 (MQ=255)                                                                     ggCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:3128863/93‑1 (MQ=255)                                                                     ggCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  >  1:317445/1‑93 (MQ=255)                                                                     ggCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:3417397/93‑1 (MQ=255)                                                                     ggCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:1872397/93‑1 (MQ=255)                                                                                       |                                                                                 GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GC‑TGGCG‑‑‑‑A‑‑‑‑ACT‑GG‑CTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458

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