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%	260.9 / NA	83	K158N (AAG→AAC)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (34/49);  total (34/49)

CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAA‑G‑GCTGGCG‑‑‑‑AA‑‑‑‑‑CT‑GGCTACAG‑CTGAAGGATTG‑ATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458                                                                                               |                                                                                    cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑g                                                                                      <  1:5707178/93‑1 (MQ=255)  tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAAGCCCTGCGCTGGAA‑gc                                                                                     >  1:773125/1‑92 (MQ=255)   cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑gcg                                                                                    <  2:4574996/93‑3 (MQ=255)   cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑gcg                                                                                    <  1:3195704/93‑3 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGGTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                                   >  1:4813315/1‑90 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                                   <  1:2019961/93‑4 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                                   <  1:4623304/93‑4 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                                   <  2:3083575/93‑4 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                                   <  2:3097571/93‑4 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑gcgg                                                                                   <  2:3282794/93‑4 (MQ=255)    gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAAGCCATGCGCTGGAA‑gcgg                                                                                   >  1:4041485/1‑90 (MQ=255)     ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑gcggt                                                                                  <  2:3738063/93‑5 (MQ=255)          gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑a                                                                         <  2:722109/93‑1 (MQ=255)          gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑a                                                                         >  1:882407/1‑93 (MQ=255)          gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑a                                                                         >  2:2003960/1‑93 (MQ=255)          gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGGCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑a                                                                         >  2:4689360/1‑93 (MQ=255)           cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aa                                                                        <  2:2103880/93‑1 (MQ=255)           cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aa                                                                        <  2:4357872/93‑1 (MQ=255)             aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                                      <  2:4328059/93‑3 (MQ=255)             aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                                      <  2:219070/93‑3 (MQ=255)             aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                                      <  2:5143571/93‑3 (MQ=255)             aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                                      >  1:4242722/1‑91 (MQ=255)             aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                                      <  2:3925526/93‑3 (MQ=255)             aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aaga                                                                      >  1:3268770/1‑91 (MQ=255)             aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCACTGGCTGGAAGG‑GGTGGCG‑‑‑‑aaaa                                                                      >  1:1410130/1‑91 (MQ=255)              tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aagaa                                                                     >  1:2474435/1‑90 (MQ=255)              tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aagaa                                                                     >  2:2923090/1‑90 (MQ=255)               cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aagaaa                                                                    <  2:2759649/93‑5 (MQ=255)               cgcgTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑aagaaa                                                                    <  1:3617388/93‑5 (MQ=255)                                               tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                 >  2:5367497/1‑93 (MQ=255)                                               tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tt                                 <  1:3067422/93‑1 (MQ=255)                                                cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                                >  2:2690426/1‑92 (MQ=255)                                                cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                                >  2:2607369/1‑92 (MQ=255)                                                cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑ttt                                <  1:3423220/93‑2 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               >  2:4708075/1‑91 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               <  1:1051706/93‑3 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               >  2:1156237/1‑91 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               >  1:5995372/1‑91 (MQ=255)                                                 tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttg                               >  1:5951006/1‑91 (MQ=255)                                                  ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑tttgg                              >  2:630739/1‑90 (MQ=255)                                                          tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                      <  2:349093/93‑1 (MQ=255)                                                          tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                      <  2:5970655/93‑1 (MQ=255)                                                          tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                      <  1:100493/93‑1 (MQ=255)                                                          tGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGa                      >  2:6045798/1‑93 (MQ=255)                                                           ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt                     <  2:2148939/93‑1 (MQ=255)                                                           ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCG‑‑‑‑AAGAAACCG‑GGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGAt                     <  2:1686766/93‑1 (MQ=255)                                                            gcgggcggCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg                    <  2:2236559/93‑1 (MQ=255)                                                            gcgggcggCGTGTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATg                    <  2:790947/93‑1 (MQ=255)                                                             cgggcggCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGa                   <  1:5478423/93‑1 (MQ=255)                                                              gggcggCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                  <  1:2635131/93‑1 (MQ=255)                                                              gggcggCGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGaa                  <  2:3445771/93‑1 (MQ=255)                                                               ggcggcGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                 >  1:3993000/1‑93 (MQ=255)                                                               ggcggcGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                 <  2:2772791/93‑1 (MQ=255)                                                               ggcggcGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                 <  1:4957056/93‑1 (MQ=255)                                                               ggcggcGTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAac                 <  2:331868/93‑1 (MQ=255)                                                                gcggcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg                <  1:2957114/93‑1 (MQ=255)                                                                gcggcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacg                >  1:4814006/1‑93 (MQ=255)                                                                 cggcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacga               <  1:2546223/93‑1 (MQ=255)                                                                 cggcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacga               <  2:2166732/93‑1 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCGTTTT‑‑GTGACGGGATGAacgac              >  2:5773184/1‑93 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac              <  2:4279944/93‑1 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac              >  2:3119511/1‑93 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac              <  2:4498313/93‑1 (MQ=255)                                                                  ggcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgac              >  2:5281397/1‑93 (MQ=255)                                                                   gcgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAacgacg             >  2:4300379/1‑93 (MQ=255)                                                                    cgTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGc            >  2:5938791/1‑93 (MQ=255)                                                                     gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGcc           <  1:4418546/93‑1 (MQ=255)                                                                     gTTTCAAAAGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGcc           <  2:4808827/93‑1 (MQ=255)                                                                       ttcaaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa         <  2:2580522/89‑1 (MQ=255)                                                                       ttcaaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa         >  2:5932837/5‑93 (MQ=255)                                                                       ttcaaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGa         >  1:6270373/5‑93 (MQ=255)                                                                        tcaaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaa        >  2:476989/4‑93 (MQ=255)                                                                         caaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  2:723500/91‑1 (MQ=255)                                                                         caaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  2:808001/91‑1 (MQ=255)                                                                         caaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  1:3434171/91‑1 (MQ=255)                                                                         caaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  1:995345/91‑1 (MQ=255)                                                                         caaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  1:5766444/91‑1 (MQ=255)                                                                         caaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       <  1:4553544/91‑1 (MQ=255)                                                                         caaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       >  1:4825533/3‑93 (MQ=255)                                                                         caaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGaaa       >  1:4825534/3‑93 (MQ=255)                                                                          aaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      >  2:4130459/2‑93 (MQ=255)                                                                          aaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      >  1:5825289/2‑93 (MQ=255)                                                                          aaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      >  2:3678743/2‑93 (MQ=255)                                                                          aaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      <  2:2946514/92‑1 (MQ=255)                                                                          aaaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAAt      >  2:1526869/2‑93 (MQ=255)                                                                           aaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGCAATg     >  2:2085646/1‑93 (MQ=255)                                                                           aaaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATg     >  1:290367/1‑93 (MQ=255)                                                                            aaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:6192030/93‑1 (MQ=255)                                                                            aaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc    <  1:2935716/93‑1 (MQ=255)                                                                            aaGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGc    >  2:3848605/1‑93 (MQ=255)                                                                             aGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:1197056/93‑1 (MQ=255)                                                                             aGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt   >  2:6045723/1‑93 (MQ=255)                                                                             aGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:193059/93‑1 (MQ=255)                                                                             aGGCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCt   <  1:3735094/93‑1 (MQ=255)                                                                              ggCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  2:2269624/93‑1 (MQ=255)                                                                              ggCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:5385678/93‑1 (MQ=255)                                                                              ggCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:241840/93‑1 (MQ=255)                                                                              ggCCATGCGCTGGAA‑GCGGTGGCGAAGAAA‑‑‑‑‑CTGGGCTACAGCCTGAAGGATTGTATTGCG‑‑TTTGGTGACGGGATGAACGACGCCGAAATGCTg  <  1:4092752/93‑1 (MQ=255)                                                                                               |                                                                                    CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAA‑G‑GCTGGCG‑‑‑‑AA‑‑‑‑‑CT‑GGCTACAG‑CTGAAGGATTG‑ATTGCG‑‑TT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458

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