BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,749,776	T→C	A693A (GCA→GCG)	b1673 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,749,776	0	T	C	99.2%	418.9 / ‑4.5	132	A693A (GCA→GCG)	b1673	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (60/71); total (60/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

Only 100 of 136 total aligned reads displayed.

CGCCTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCTGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTCTTT  >  NC_000913/1749686‑1749865

cGCCTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGtgccgc                                                                                         >  2:2662356/1‑93 (MQ=255)

   cTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAg                                                                                      >  2:3064083/1‑93 (MQ=255)

   cTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAg                                                                                      >  1:625602/1‑93 (MQ=255)

    tCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCCGa                                                                                     >  2:4663819/1‑93 (MQ=255)

    tCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGa                                                                                     <  1:789625/93‑1 (MQ=255)

    tCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGa                                                                                     >  2:1162616/1‑93 (MQ=255)

    tCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGa                                                                                     >  2:167184/1‑93 (MQ=255)

    tCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGa                                                                                     >  2:453516/1‑93 (MQ=255)

      gTCAACTTCGCTGCCTTTCGCTAACGTCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATc                                                                                   <  1:4429716/93‑1 (MQ=255)

      gTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATc                                                                                   <  2:410937/93‑1 (MQ=255)

       tCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcg                                                                                  <  1:4146333/93‑1 (MQ=255)

       tCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcg                                                                                  <  1:4146334/93‑1 (MQ=255)

       tCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcg                                                                                  >  1:4629125/1‑93 (MQ=255)

        cAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcgg                                                                                 <  2:3782953/93‑1 (MQ=255)

         aaCTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcggc                                                                                <  1:3570822/93‑1 (MQ=255)

           cTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATcggcgg                                                                              >  1:2818662/1‑93 (MQ=255)

            ttCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCggcggc                                                                             >  1:3696850/1‑93 (MQ=255)

             tCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCa                                                                            >  1:1740191/1‑93 (MQ=255)

              cgcgGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCaa                                                                           <  1:1038740/93‑1 (MQ=255)

              cGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCaa                                                                           <  2:775499/93‑1 (MQ=255)

               gCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAat                                                                          >  1:4388578/1‑93 (MQ=255)

                 tGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAatat                                                                        >  1:1573936/1‑93 (MQ=255)

                 tGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAatat                                                                        >  2:8632/1‑93 (MQ=255)

                  gACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATc                                                                       <  2:4540353/93‑1 (MQ=255)

                   aCTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCt                                                                      <  2:96571/93‑1 (MQ=255)

                   aCTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCt                                                                      <  2:96573/93‑1 (MQ=255)

                     tttCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTc                                                                    <  2:2783750/93‑1 (MQ=255)

                     tttCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTc                                                                    >  1:1053156/1‑93 (MQ=255)

                      ttCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTcc                                                                   <  2:2718364/93‑1 (MQ=255)

                      ttCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTcc                                                                   <  1:4016950/93‑1 (MQ=255)

                      ttCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTcc                                                                   <  1:3270018/93‑1 (MQ=255)

                      ttCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTcc                                                                   <  1:4519348/93‑1 (MQ=255)

                       tCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCa                                                                  <  1:3670158/93‑1 (MQ=255)

                       tCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCa                                                                  >  1:2928440/1‑93 (MQ=255)

                        cGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAg                                                                 <  1:2309531/93‑1 (MQ=255)

                        cGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAg                                                                 >  1:684806/1‑93 (MQ=255)

                         gttAACGGCAATTCATCCTGATGGTTCATCTTACGCATTCCTTACGCAGGCAGTAGATTGTGTGCCTCGAGATCGGCGGCAATATCTTCCAGc                                                                <  2:1800579/91‑1 (MQ=255)

                           tAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGccc                                                              >  1:2927476/1‑93 (MQ=255)

                            aaCGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGcccc                                                             <  1:2021547/93‑1 (MQ=255)

                             aCGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGccccc                                                            >  1:1327151/1‑92 (MQ=255)

                              cGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAg                                                           <  1:1384437/93‑1 (MQ=255)

                                gCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCTGCCAGATCGGCGGCAATATCTTCCAGCCCCAGAc                                                         <  2:1587986/93‑1 (MQ=255)

                                  cattcattcGGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGc                                                       <  1:1583164/92‑1 (MQ=255)

                                  aaTTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGc                                                       <  1:1583163/93‑1 (MQ=255)

                                  aaTTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAAATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGc                                                       >  1:4290969/1‑93 (MQ=255)

                                    ttCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTg                                                     <  1:3736093/93‑1 (MQ=255)

                                    ttCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTg                                                     <  2:4156915/93‑1 (MQ=255)

                                     tCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGc                                                    <  1:2725618/93‑1 (MQ=255)

                                     tCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGc                                                    >  1:1823491/1‑93 (MQ=255)

                                     tCATCCTGATGGTTCATCTTACCCATTCCTCACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGc                                                    <  2:3920171/93‑1 (MQ=255)

                                      cATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCa                                                   >  1:2228868/1‑93 (MQ=255)

                                       aTCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCaa                                                  >  1:4407453/1‑93 (MQ=255)

                                       aTCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCaa                                                  <  2:2014760/93‑1 (MQ=255)

                                        tCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGGTGCAAt                                                 >  2:3259899/1‑93 (MQ=255)

                                          cTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGt                                               <  1:1318382/93‑1 (MQ=255)

                                          cTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGt                                               >  2:3988586/1‑93 (MQ=255)

                                          cTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGt                                               <  2:3559390/93‑1 (MQ=255)

                                          cTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCAAGATCGGCGGCAATAAATTCCAGCCCCAGACGCTGCCATGt                                               >  2:860380/1‑93 (MQ=255)

                                           tGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATggt                                              >  2:2640669/1‑91 (MQ=255)

                                           tGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTACAGCCTCAAAAGCTGCAATCtt                                              >  1:1930196/1‑93 (MQ=255)

                                            gATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCTCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGttt                                             <  2:3340292/93‑1 (MQ=255)

                                            gATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGttt                                             >  2:135856/1‑93 (MQ=255)

                                            gATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGttt                                             >  2:135855/1‑93 (MQ=255)

                                             aTGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTc                                            <  2:295423/93‑1 (MQ=255)

                                             aTGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTc                                            <  2:4620362/93‑1 (MQ=255)

                                              tGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTcg                                           <  1:2242183/93‑1 (MQ=255)

                                               ggTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTcgc                                          <  1:3843816/93‑1 (MQ=255)

                                               ggTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTcgc                                          <  2:2069127/93‑1 (MQ=255)

                                                  tCATCTTACCCATGCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCggg                                       <  1:98270/93‑1 (MQ=255)

                                                   cATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGt                                      <  2:2817/93‑1 (MQ=255)

                                                     tCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTg                                    >  2:3487508/1‑93 (MQ=255)

                                                      cTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTgg                                   <  2:4522568/93‑1 (MQ=255)

                                                       ttccccATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGa                                  <  1:2851706/90‑1 (MQ=255)

                                                       ttACCCATTCCTTACGCAGGCAGTAGATTGTGTCCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGa                                  <  1:2449529/93‑1 (MQ=255)

                                                         aCCCATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGAGGGTTGGACa                                >  1:1127417/1‑93 (MQ=255)

                                                           ccATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAAc                              >  1:3915165/1‑93 (MQ=255)

                                                            cATTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAAcc                             <  1:1845066/93‑1 (MQ=255)

                                                             aTTCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCa                            >  2:4378415/1‑93 (MQ=255)

                                                               tCCTTACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAg                          >  2:3161055/1‑93 (MQ=255)

                                                                   tACGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCtgt                      <  2:1159652/93‑1 (MQ=255)

                                                                    aCGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCtgtg                     <  2:3449752/93‑1 (MQ=255)

                                                                     cGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTgg                    <  2:4062801/93‑1 (MQ=255)

                                                                     cGCAGGCAGTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTgg                    <  2:43174/93‑1 (MQ=255)

                                                                      gcagccagTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTggg                   <  2:1631809/93‑1 (MQ=255)

                                                                       caggcagTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGt                  <  1:2083833/93‑1 (MQ=255)

                                                                        aggcagTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAATCTGTGGGTc                 <  1:955448/93‑1 (MQ=255)

                                                                         ggcagTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTcc                >  1:516414/1‑93 (MQ=255)

                                                                         ggcagGAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTcc                <  2:513955/93‑1 (MQ=255)

                                                                          gcagTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTccc               <  2:1498411/93‑1 (MQ=255)

                                                                           cagTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCa              <  1:3545396/93‑1 (MQ=255)

                                                                             gTAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGGCGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGc            <  2:2006341/93‑1 (MQ=255)

                                                                              tAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGcc           >  2:2356789/1‑93 (MQ=255)

                                                                              tAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGcc           <  1:519003/93‑1 (MQ=255)

                                                                              tAGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGcc           >  1:3891285/1‑93 (MQ=255)

                                                                               aGATTGTGTGCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGccc          >  1:2464471/1‑93 (MQ=255)

                                                                                    gtgtgCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTc     >  2:2060696/1‑93 (MQ=255)

                                                                                    gtgtgCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTc     <  2:1772859/93‑1 (MQ=255)

                                                                                     tgtgCCGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTCt    <  2:1660675/93‑1 (MQ=255)

                                                                                       tgccgcCAGATCGGCTGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTCttt  >  2:1185369/1‑93 (MQ=255)

                                                                                       tgccgcCAGATCGGCGGCAATATCTTCCAGCCCCTGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTCttt  <  2:2142199/93‑1 (MQ=255)

CGCCTCGTCAACTTCGCTGACTTTCGCTAACGGCAATTCATCCTGATGGTTCATCTTACCCATTCCTTACGCAGGCAGTAGATTGTGTGCTGCCAGATCGGCGGCAATATCTTCCAGCCCCAGACGCTGCAATGTTTCGCGGGTTGGACAACCAAGCTGTGGGTCCCAGCCCATTTCTTT  >  NC_000913/1749686‑1749865

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