BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,845,852	T→C	T61A (ACC→GCC)	selD ←	selenophosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,845,852	0	T	C	99.5%	586.4 / NA	190	T61A (ACC→GCC)	selD	selenophosphate synthase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (95/94); minor base A (1/0); total (96/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

GCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAT  >  NC_000913/1845763‑1845941

 cTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCGTAAAGAAGTCGGTGGTACTGATAACGCTGGAGCCa                                                                                       >  1:1274891/1‑93 (MQ=255)

 cTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCa                                                                                       <  2:5958735/93‑1 (MQ=255)

 cTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCa                                                                                       >  1:2303813/1‑93 (MQ=255)

  tGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAt                                                                                      <  2:5072460/93‑1 (MQ=255)

   gATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAtt                                                                                     >  1:5709067/1‑93 (MQ=255)

    aTGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTg                                                                                    <  1:5990538/93‑1 (MQ=255)

       gCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGccc                                                                                 >  2:6697906/1‑93 (MQ=255)

        ggTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCa                                                                                <  2:3808760/92‑1 (MQ=255)

        cGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCa                                                                                <  1:113513/93‑1 (MQ=255)

          ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCCTTGGCCCGa                                                                              >  2:3215717/1‑93 (MQ=255)

          ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                              >  1:2068634/1‑93 (MQ=255)

           tAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGAt                                                                             <  2:7101124/93‑1 (MQ=255)

           tAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGAt                                                                             <  2:4545143/93‑1 (MQ=255)

            aGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATc                                                                            >  2:7732133/1‑93 (MQ=255)

             gTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCg                                                                           <  1:5866968/93‑1 (MQ=255)

              tCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGt                                                                          >  2:5161319/1‑93 (MQ=255)

                gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTac                                                                        >  1:3449271/1‑93 (MQ=255)

                  cgcAATGCGGCCAAAATCGGAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTacac                                                                      >  1:3988428/1‑93 (MQ=255)

                   gcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAcc                                                                     >  1:3297923/1‑93 (MQ=255)

                    cAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccg                                                                    >  2:2645870/1‑93 (MQ=255)

                    cAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccg                                                                    >  1:5060534/1‑93 (MQ=255)

                     caTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgc                                                                   <  1:44361/92‑1 (MQ=255)

                     caTGCCGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgc                                                                   <  1:44358/92‑1 (MQ=255)

                     aaTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgc                                                                   <  2:685210/93‑1 (MQ=255)

                     aaTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgc                                                                   >  1:3968444/1‑93 (MQ=255)

                      aTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgcc                                                                  >  2:1589507/1‑93 (MQ=255)

                       tGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgccg                                                                 >  2:802981/1‑93 (MQ=255)

                        gCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACcgccgc                                                                >  1:3921622/1‑93 (MQ=255)

                          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAt                                                              >  1:2443359/1‑93 (MQ=255)

                           gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtc                                                             >  1:2986020/1‑93 (MQ=255)

                           gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtc                                                             <  1:10813/93‑1 (MQ=255)

                           gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtc                                                             <  1:6899588/93‑1 (MQ=255)

                              aaaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtc                                                          >  2:7355216/1‑93 (MQ=255)

                               aaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtcg                                                         <  2:1323198/93‑1 (MQ=255)

                                aaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTcgc                                                        <  1:1666746/93‑1 (MQ=255)

                                 aTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTcgcg                                                       >  2:1199985/1‑93 (MQ=255)

                                  tCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTTGTACTGATAACGCTGGCGCCATTGCGCAGATCGTACACCGCCGCATCGTCGCgg                                                      <  2:4834423/93‑1 (MQ=255)

                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    >  1:6134868/1‑93 (MQ=255)

                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                    <  1:4241082/93‑1 (MQ=255)

                                     aaaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGtt                                                   >  2:2578080/1‑93 (MQ=255)

                                     aaaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGtt                                                   <  2:7134746/93‑1 (MQ=255)

                                      aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGttt                                                  <  2:2814293/93‑1 (MQ=255)

                                       aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                 <  1:1709829/93‑1 (MQ=255)

                                        ggcttATCAACGATCGGCATAAAGAAGTCGGTGGTACTGACAACGCTGGCGCCATTGCCCAGATTGTACACCGCCGCATCGTCGCGGGTTTCa                                                <  2:5271143/90‑1 (MQ=255)

                                        ggATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCa                                                >  2:1326259/1‑93 (MQ=255)

                                         gATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAt                                               >  1:5952360/1‑93 (MQ=255)

                                           ttATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTa                                             <  1:5506391/93‑1 (MQ=255)

                                            tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                            <  2:145106/93‑1 (MQ=255)

                                              tCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAccc                                          >  2:6226832/1‑93 (MQ=255)

                                              tCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAccc                                          <  2:2942533/93‑1 (MQ=255)

                                               cAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCca                                         <  1:4842100/93‑1 (MQ=255)

                                               cAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCca                                         >  2:2793258/1‑93 (MQ=255)

                                                aaCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCcac                                        >  2:1041586/1‑93 (MQ=255)

                                                aaCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCcac                                        >  1:3753885/1‑93 (MQ=255)

                                                 aCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCcaca                                       >  1:927992/1‑93 (MQ=255)

                                                  cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                      <  2:3988427/93‑1 (MQ=255)

                                                   gATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAg                                     >  2:3030637/1‑93 (MQ=255)

                                                     tCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCa                                   >  1:780268/1‑93 (MQ=255)

                                                      cGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCaa                                  <  1:57487/93‑1 (MQ=255)

                                                       ggCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCaaa                                 <  2:7889072/93‑1 (MQ=255)

                                                        gCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                                <  2:248248/93‑1 (MQ=255)

                                                        gCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                                >  1:6526241/1‑93 (MQ=255)

                                                         cATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAtt                               <  2:7291041/93‑1 (MQ=255)

                                                         cATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAtt                               <  1:2636312/93‑1 (MQ=255)

                                                          aTAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTc                              <  1:1530448/93‑1 (MQ=255)

                                                           tAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                             >  2:1149456/1‑93 (MQ=255)

                                                           tAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                             >  2:2355558/1‑93 (MQ=255)

                                                            aaaGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCgg                            >  1:1353715/1‑93 (MQ=255)

                                                            aaaGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCgg                            >  2:1808130/1‑93 (MQ=255)

                                                            aaaGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCgg                            >  1:6315782/1‑93 (MQ=255)

                                                             aagaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGa                           >  1:1618576/1‑93 (MQ=255)

                                                             aagaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGa                           <  2:6623438/93‑1 (MQ=255)

                                                                 agTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaa                       <  2:6404000/93‑1 (MQ=255)

                                                                  gTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaac                      <  2:2172837/93‑1 (MQ=255)

                                                                  gTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaac                      <  2:3653596/93‑1 (MQ=255)

                                                                   tCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaaca                     <  2:2598893/93‑1 (MQ=255)

                                                                     ggtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACaaa                   <  1:5415467/93‑1 (MQ=255)

                                                                     ggtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTCTCATTACCCACAAGCAAATTCGGATCAACaaa                   <  2:6586408/93‑1 (MQ=255)

                                                                      gtggtACTGATAACTCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAAc                  <  2:2037550/93‑1 (MQ=255)

                                                                      gtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAAc                  <  2:873767/93‑1 (MQ=255)

                                                                        ggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACtt                >  2:5469786/1‑93 (MQ=255)

                                                                         gtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc               >  2:7477962/1‑93 (MQ=255)

                                                                         gtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc               <  1:520827/93‑1 (MQ=255)

                                                                            cTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGcc            <  2:196934/93‑1 (MQ=255)

                                                                             tGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCt           >  2:6768719/1‑93 (MQ=255)

                                                                              gATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTg          >  1:1898488/1‑93 (MQ=255)

                                                                              gATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGGTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTg          >  2:6923048/1‑93 (MQ=255)

                                                                               aTAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGc         <  2:7587495/93‑1 (MQ=255)

                                                                               aTAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGc         <  1:7437348/93‑1 (MQ=255)

                                                                               aTAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGc         <  1:7670968/93‑1 (MQ=255)

                                                                                 aaCGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTc       <  1:5774812/93‑1 (MQ=255)

                                                                                 aaCGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTc       <  2:5885525/93‑1 (MQ=255)

                                                                                   cGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc     <  2:725129/93‑1 (MQ=255)

                                                                                   cGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc     >  2:3731167/1‑93 (MQ=255)

                                                                                   cGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc     <  1:6844982/93‑1 (MQ=255)

                                                                                     cTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa   <  1:1699206/93‑1 (MQ=255)

                                                                                     cTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa   <  2:2851018/93‑1 (MQ=255)

                                                                                     cTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa   >  1:7207974/1‑93 (MQ=255)

                                                                                      tGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAATTTCGCCTGCTCACTAt  >  1:7034281/1‑93 (MQ=255)

                                                                                      tGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAt  >  1:2796738/1‑93 (MQ=255)

GCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAT  >  NC_000913/1845763‑1845941

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: `‑`