BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,264,376	T→C	62.3%	I221I (ATT→ATC)	hemA →	glutamyl tRNA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,264,376	0	T	C	62.3%	103.1 / 112.6	117	I221I (ATT→ATC)	hemA	glutamyl tRNA reductase
Reads supporting (aligned to +/- strand): ref base T (25/18); new base C (36/36); total (61/56)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 117 total aligned reads displayed.

GCGCGGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGCCAGCC  >  NC_000913/1264279‑1264470

gcgcGGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCt                                                                                              <  1:1961143/100‑1 (MQ=255)

  gcgGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTgg                                                                                            <  2:3587016/100‑1 (MQ=255)

  gcgGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTgg                                                                                            <  2:1884742/100‑1 (MQ=255)

  gcgGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTgg                                                                                            <  2:486989/100‑1 (MQ=255)

   cgGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGc                                                                                           <  2:1949507/100‑1 (MQ=255)

   cgGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGc                                                                                           <  2:2496405/100‑1 (MQ=255)

   cgGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGc                                                                                           <  1:2201027/100‑1 (MQ=255)

    gggCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCa                                                                                          >  1:3416748/1‑100 (MQ=255)

    gggCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCa                                                                                          <  1:944881/100‑1 (MQ=255)

    gggCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCa                                                                                          <  2:160653/100‑1 (MQ=255)

     ggCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAg                                                                                         >  1:1610096/1‑100 (MQ=255)

       ccaaaCTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTACCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGa                                                                                        <  2:3578235/97‑1 (MQ=255)

       cGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGa                                                                                        >  2:3182606/1‑99 (MQ=255)

        gAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATg                                                                                      <  1:964705/100‑1 (MQ=255)

         aaaCTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGa                                                                                     >  1:2957581/1‑100 (MQ=255)

              aTCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTc                                                                                 >  2:426237/1‑99 (MQ=255)

              aTCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCg                                                                                >  1:44651/1‑100 (MQ=255)

              aTCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCg                                                                                >  1:106432/1‑100 (MQ=255)

              aTCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCg                                                                                >  1:1077189/1‑100 (MQ=255)

                cGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGgc                                                                              >  2:1138645/1‑100 (MQ=255)

                cGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGgc                                                                              >  1:2537741/1‑100 (MQ=255)

                 gAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGgcg                                                                             >  1:2693056/1‑100 (MQ=255)

                 gAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGgcg                                                                             >  2:3721080/1‑100 (MQ=255)

                 gAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGgcg                                                                             >  2:1719286/1‑100 (MQ=255)

                 gAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGgcg                                                                             >  2:3775144/1‑100 (MQ=255)

                  aGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGgcgc                                                                            >  2:241451/1‑100 (MQ=255)

                  aGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGgcgc                                                                            >  2:593195/1‑100 (MQ=255)

                   gCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGcgcg                                                                           <  1:1526765/100‑1 (MQ=255)

                    cTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCgg                                                                          >  2:2135313/1‑100 (MQ=255)

                     tggtggCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATAGCCAACCGCACACGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCgg                                                                          >  2:2296618/1‑99 (MQ=255)

                       gtggCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAg                                                                       <  1:3020082/100‑1 (MQ=255)

                       gtggCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAg                                                                       <  2:2860787/100‑1 (MQ=255)

                       gtggCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAg                                                                       >  1:2246788/1‑100 (MQ=255)

                        tggCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGt                                                                      >  1:3122859/1‑100 (MQ=255)

                        tggCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGt                                                                      <  1:455854/100‑1 (MQ=255)

                           cgcgTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGAt                                                                   <  2:370505/100‑1 (MQ=255)

                             cgTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTg                                                                 >  1:1387031/1‑100 (MQ=255)

                                cATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGccc                                                              >  2:2134585/1‑100 (MQ=255)

                                 aTCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCt                                                             >  2:725988/1‑100 (MQ=255)

                                    tGCGCGAACACAAAGTACAGAAGATGATTATCGCCAAGCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAg                                                          >  2:2203020/1‑100 (MQ=255)

                                    tGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAGGGGATTGCCCTGAg                                                          >  2:1918745/1‑100 (MQ=255)

                                    tGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAg                                                          >  1:646301/1‑100 (MQ=255)

                                    tGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAg                                                          <  1:1693026/100‑1 (MQ=255)

                                     gtgcgaaCACAAAGTACAGAAGATGATTATCGCCAACCGCACGCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGT                                                         >  1:2799618‑M1/3‑100 (MQ=255)

                                      cgcgAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATGCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTG                                                        <  2:374123‑M1/100‑1 (MQ=255)

                                      cgcgAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTG                                                        >  2:4025458‑M1/1‑100 (MQ=255)

                                       gcgAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGA                                                       >  1:3991365‑M1/1‑100 (MQ=255)

                                       gcgAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGA                                                       <  1:593200‑M1/100‑1 (MQ=255)

                                        cgAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGAT                                                      >  1:3094875‑M1/1‑100 (MQ=255)

                                        cgAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGAT                                                      >  2:2730680‑M1/1‑100 (MQ=255)

                                        cgAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAACTGATTGCCCTGCGCGAT                                                      >  2:3955121‑M1/1‑100 (MQ=255)

                                         gAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATA                                                     <  2:3000010‑M1/100‑1 (MQ=255)

                                         gAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATA                                                     <  1:607636‑M1/100‑1 (MQ=255)

                                          aaCACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATAT                                                    <  1:846503‑M1/100‑1 (MQ=255)

                                           acacAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATC                                                   <  2:2537727‑M1/100‑1 (MQ=255)

                                            cacaAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCG                                                  >  1:3419379‑M1/1‑100 (MQ=255)

                                             acaAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATcga                                                 >  2:2998007/1‑100 (MQ=255)

                                              caAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATcgac                                                <  1:3245419/100‑1 (MQ=255)

                                              caAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATcgac                                                <  2:815331/100‑1 (MQ=255)

                                               aaaGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATcgacg                                               >  1:2860194/1‑100 (MQ=255)

                                                     cAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTc                                         >  2:3596287/1‑100 (MQ=255)

                                                      agaagaTGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCAGAAGAGATTGCCCTGAGTGATATCGACGAACGTc                                         >  2:1398954/1‑99 (MQ=255)

                                                       gaagaTGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTg                                       >  2:2335858/1‑100 (MQ=255)

                                                            tgatTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGa                                  >  2:1295346/1‑100 (MQ=255)

                                                                tATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCc                              >  1:4084423‑M2/1‑99 (MQ=255)

                                                                tATCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCc                              >  1:93467‑M2/1‑99 (MQ=255)

                                                                  tCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCcga                            <  1:3656350‑M2/100‑4 (MQ=255)

                                                                  tCGCCAACCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCcga                            >  1:1235223‑M2/1‑97 (MQ=255)

                                                                    gCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCcgata                          <  1:813273‑M2/100‑6 (MQ=255)

                                                                     ccAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCcgatat                         >  2:1049516‑M2/1‑94 (MQ=255)

                                                                     ccAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCcgatat                         <  1:3331194‑M2/100‑7 (MQ=255)

                                                                       aaCCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCcgatatca                       >  1:3340273‑M2/1‑92 (MQ=255)

                                                                       aaCCGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCcgatatca                       >  2:2380499‑M2/1‑92 (MQ=255)

                                                                         ccGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATatcatc                     >  1:343432/1‑100 (MQ=255)

                                                                         ccGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAATGATATCGACGAACGTCTGGGCGAAGCcgatagca                       <  2:4072431‑M2/98‑9 (MQ=255)

                                                                         ccGCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATatcatc                     <  2:1491027/100‑1 (MQ=255)

                                                                           gCACTCGCGAATGTGCCCAAATCCTGGCAGATGAAGTCGGCGAGGAAGTGATTGCCCTGAGGCATATCGACGAACGTCTGCGCGAATCCGATatcatcat                   <  2:2550677/100‑1 (MQ=255)

                                                                           gCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGgcttcgcgcagacgttcgtcgatatcactca                   <  2:3106871‑M1/100‑32 (MQ=255)

                                                                           gCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGgcttcgcgcagacgttcgtcgatatcactca                   <  1:824844‑M1/100‑32 (MQ=255)

                                                                           gCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATatcatcat                   >  1:1046003/1‑100 (MQ=255)

                                                                           gCACTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATatcat                      >  2:4023947/1‑97 (MQ=255)

                                                                             aCTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGgcttcgcgcagacgttcgtcgatatcactcagg                 <  1:1716908‑M1/100‑34 (MQ=255)

                                                                             aCTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATAtcatcatca                 <  2:2925837/100‑1 (MQ=255)

                                                                              cTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGgcttcgcgcagacgttcgtcgatatcactcaggg                <  1:3726723‑M1/100‑35 (MQ=255)

                                                                              cTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAg                >  2:663763/1‑100 (MQ=255)

                                                                              cTCGCGAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGgcttcgcgcagacgttcgtcgatatcactcaggg                <  1:391494‑M1/100‑35 (MQ=255)

                                                                                cgcgAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGtt              <  2:3410338/100‑1 (MQ=255)

                                                                                cgcgAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGtt              >  2:3046692/1‑100 (MQ=255)

                                                                                cgcgAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGtt              <  2:3377499/100‑1 (MQ=255)

                                                                                cgcgAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGtt              <  1:1886027/100‑1 (MQ=255)

                                                                                   gAACGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCa           <  2:3970784/100‑1 (MQ=255)

                                                                                     aCGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCAcc         <  2:3795066/100‑1 (MQ=255)

                                                                                     aCGTGCCCAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCAcc         <  1:3003253/100‑1 (MQ=255)

                                                                                      cGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCg        <  2:3614345/100‑1 (MQ=255)

                                                                                          cccAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATAGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGCCAg    >  1:202777/1‑100 (MQ=255)

                                                                                          cccAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGCCAg    >  1:464510/1‑100 (MQ=255)

                                                                                           ccAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGCCAGc   >  2:3695937/1‑100 (MQ=255)

                                                                                           ccAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGAGATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGCCAg    >  1:2520125/1‑99 (MQ=255)

                                                                                            cAAATCCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGTGCGAAGCCGATATCATCATCAGTTCCACCGCCAGcc  <  1:3710183/100‑1 (MQ=255)

                                                                                              aaTGCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCTATGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCAACAGTTCCACCGCCa     <  2:1621271/95‑1 (MQ=255)

GCGCGGGCGAAACTATCGAGCTGGTGGCGCGTCATCTGCGCGAACACAAAGTACAGAAGATGATTATCGCCAACCGCACTCGCGAACGTGCCCAAATTCTGGCAGATGAAGTCGGCGCGGAAGTGATTGCCCTGAGTGATATCGACGAACGTCTGCGCGAAGCCGATATCATCATCAGTTCCACCGCCAGCC  >  NC_000913/1264279‑1264470

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`