Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,391	+C	coding (474/516 nt)	b4555 →	hypothetical protein

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

CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GC‑TGGCG‑‑‑‑A‑‑‑‑ACT‑GG‑CTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458                                                                                              |                                                                                   cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAg                                                                                     <  2:3237564/93‑1 (MQ=255)  tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                                    <  1:665885/93‑2 (MQ=255)  tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                                    <  1:4450303/93‑2 (MQ=255)  tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                                    <  2:5312058/93‑2 (MQ=255)  tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                                    <  1:1938513/93‑2 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                                  >  1:2846952/1‑90 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                                  <  1:2544476/93‑4 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                                  <  2:2866793/93‑4 (MQ=255)     ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                                 <  2:2095072/93‑5 (MQ=255)     ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                                 <  1:5572441/93‑5 (MQ=255)     ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                                 <  2:2844229/93‑5 (MQ=255)     ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggg                                                                                 >  1:975712/1‑89 (MQ=255)     ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggg                                                                                 >  2:1580045/1‑89 (MQ=255)     ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAgcggg                                                                                 >  2:2826432/1‑89 (MQ=255)          gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a                                                                       <  2:1404079/93‑1 (MQ=255)           cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a‑‑‑‑a                                                                  <  1:2722700/93‑1 (MQ=255)           cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a‑‑‑‑a                                                                  <  2:5007933/93‑1 (MQ=255)           cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCG‑‑‑‑a‑‑‑‑a                                                                  <  1:477398/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:2612051/93‑3 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               <  2:2720365/93‑3 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:6351818/93‑3 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               <  2:4945205/93‑3 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               >  1:5279004/1‑91 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:5203346/93‑3 (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                       >  1:4895227/1‑93 (MQ=255)                                                         aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg                       <  2:3643840/93‑1 (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                     <  2:1677885/93‑1 (MQ=255)                                                           ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:1612030/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                 <  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)                                                               ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  1:4668012/93‑1 (MQ=255)                                                               ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 >  2:4169501/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:3733356/93‑1 (MQ=255)                                                                 cggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:3061014/93‑1 (MQ=255)                                                                 cggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:5029312/93‑1 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  2:6349720/93‑1 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  2:6002290/93‑1 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  1:838731/1‑93 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              >  2:506903/1‑93 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  1:2279083/93‑1 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac              <  1:1387788/93‑1 (MQ=255)                                                                   gcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  1:1938896/1‑93 (MQ=255)                                                                   gcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  1:550341/1‑93 (MQ=255)                                                                   gcgTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             <  1:1704662/93‑1 (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:3054417/1‑93 (MQ=255)                                                                     gTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  2:3718784/93‑1 (MQ=255)                                                                     gTTTCAAAAGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           >  2:5795906/1‑93 (MQ=255)                                                                         caaaaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:2486652/91‑1 (MQ=255)                                                                          aaaaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:2970238/2‑93 (MQ=255)                                                                          aaaaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:269946/2‑93 (MQ=255)                                                                           aaaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  2:1611781/1‑93 (MQ=255)                                                                            aaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  1:1121834/1‑93 (MQ=255)                                                                            aaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  1:5171996/1‑93 (MQ=255)                                                                            aaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    >  1:501772/1‑93 (MQ=255)                                                                            aaGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc    <  2:4329317/93‑1 (MQ=255)                                                                             aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:6392472/1‑93 (MQ=255)                                                                             aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:2466575/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   >  1:4138537/1‑93 (MQ=255)                                                                             aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:4308666/1‑93 (MQ=255)                                                                             aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:4406304/1‑93 (MQ=255)                                                                             aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:5530256/1‑93 (MQ=255)                                                                             aGGCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   >  1:616170/1‑93 (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:3128863/93‑1 (MQ=255)                                                                              ggCCATGCGCTGGAAGCGG‑TGGCGAAGAA‑‑‑‑ACTGGG‑CTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:1872397/93‑1 (MQ=255)                                                                                              |                                                                                   CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GC‑TGGCG‑‑‑‑A‑‑‑‑ACT‑GG‑CTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458

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