BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,114,335	T→C	A156A (GCA→GCG)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,114,335	0	T	C	99.6%	747.9 / ‑5.5	222	A156A (GCA→GCG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (107/114); total (108/114)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 235 total aligned reads displayed.

GTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCT  >  NC_000913/4114245‑4114424

gTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACgcgc                                                                                         >  1:6807956/1‑93 (MQ=255)

 ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCa                                                                                        >  1:8124104/1‑93 (MQ=255)

 ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCa                                                                                        >  2:3658317/1‑93 (MQ=255)

 ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCa                                                                                        >  1:23929/1‑93 (MQ=255)

   tgtgtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACg                                                                                      <  2:6950195/92‑1 (MQ=255)

    gtgtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGc                                                                                     <  1:8454713/93‑1 (MQ=255)

      gtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTc                                                                                   >  1:2146514/1‑93 (MQ=255)

       tAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTcg                                                                                  <  2:4917316/93‑1 (MQ=255)

         aTCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTcgcg                                                                                <  2:3072514/93‑1 (MQ=255)

          tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCga                                                                               >  1:1076315/1‑93 (MQ=255)

          tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCga                                                                               <  2:1172176/93‑1 (MQ=255)

          tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCga                                                                               <  2:4833415/93‑1 (MQ=255)

           cGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCgag                                                                              <  2:4510918/93‑1 (MQ=255)

            ggTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCgaga                                                                             >  2:8536836/1‑93 (MQ=255)

             gTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGagag                                                                            >  2:3771067/1‑93 (MQ=255)

                 caTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCtt                                                                        <  1:8694530/93‑1 (MQ=255)

                  aTGGACACGGCCCTGAGTCCTTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTc                                                                       >  2:7949724/1‑93 (MQ=255)

                    ggACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGCGCCTTCca                                                                     >  1:7717509/1‑93 (MQ=255)

                     gACACGGCCCTGAGTCATTTTCCAGGTAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCcac                                                                    >  1:16970/1‑93 (MQ=255)

                     gACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCcac                                                                    <  1:2403538/93‑1 (MQ=255)

                      gcacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCcaca                                                                   <  1:1884022/92‑1 (MQ=255)

                      acacGGCCCTGAGTCATTTTCCAGATAATCCACGTATCAACCGTACAAAACAGCAATCCACCACGACGCGCACGCTCGCGAGAGCCTTCcaca                                                                   <  2:5512690/93‑1 (MQ=255)

                      acacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCcaca                                                                   <  1:3959939/93‑1 (MQ=255)

                       cacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACAt                                                                  <  1:3660062/93‑1 (MQ=255)

                       cacGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACAt                                                                  <  2:246334/93‑1 (MQ=255)

                         cGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATgg                                                                >  1:1185187/1‑93 (MQ=255)

                         cGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATgg                                                                <  1:2955257/93‑1 (MQ=255)

                         cGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATgg                                                                <  1:4627146/93‑1 (MQ=255)

                          ggCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGt                                                               >  1:1640277/1‑93 (MQ=255)

                          ggCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGt                                                               <  2:7374245/93‑1 (MQ=255)

                           gCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTc                                                              <  1:8177550/93‑1 (MQ=255)

                            cccTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCg                                                             <  2:4727067/93‑1 (MQ=255)

                            cccTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCg                                                             >  2:4358902/1‑93 (MQ=255)

                            cccTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACACGGTCg                                                             <  1:4275191/93‑1 (MQ=255)

                             ccTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGa                                                            >  1:5504658/1‑93 (MQ=255)

                             ccTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGa                                                            <  2:8140612/93‑1 (MQ=255)

                             ccTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGa                                                            >  2:1626329/1‑93 (MQ=255)

                              cTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAg                                                           >  2:6318909/1‑93 (MQ=255)

                              cTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAg                                                           >  1:811255/1‑93 (MQ=255)

                              cTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAg                                                           >  1:7766303/1‑93 (MQ=255)

                               tGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCTCGAGAGCCTTCCACATGGTCGAgg                                                          >  1:2258056/1‑93 (MQ=255)

                               tGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAgg                                                          <  2:8704724/93‑1 (MQ=255)

                                 aGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGAt                                                        <  1:3194202/93‑1 (MQ=255)

                                 aGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGAt                                                        <  1:5206590/93‑1 (MQ=255)

                                 aGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGAt                                                        <  2:5655677/93‑1 (MQ=255)

                                 aGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGAt                                                        <  2:7561183/93‑1 (MQ=255)

                                  gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc                                                       >  2:3967811/1‑93 (MQ=255)

                                   tCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATcc                                                      <  1:2845466/93‑1 (MQ=255)

                                     aTTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                    >  2:3879535/1‑93 (MQ=255)

                                     aTTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                    >  1:5680493/1‑93 (MQ=255)

                                     aTTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                    >  2:1602282/1‑93 (MQ=255)

                                     aTTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcac                                                    <  2:793835/93‑1 (MQ=255)

                                      ttttCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcact                                                   <  2:2548391/93‑1 (MQ=255)

                                       tttCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcactt                                                  >  2:5206150/1‑93 (MQ=255)

                                            aGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcactt                                             >  1:4712154/1‑93 (MQ=255)

                                            aGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcactt                                             <  2:5377703/93‑1 (MQ=255)

                                            aGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcactt                                             <  1:4409763/93‑1 (MQ=255)

                                            aGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGCATCcacttcactt                                             >  1:7496992/1‑93 (MQ=255)

                                             gATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACttt                                            >  2:5885541/1‑93 (MQ=255)

                                             gATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACttt                                            <  1:1431532/93‑1 (MQ=255)

                                              aTAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTg                                           >  2:8183123/1‑93 (MQ=255)

                                                  gCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGc                                       >  1:3605605/1‑93 (MQ=255)

                                                  gCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGc                                       >  1:6411020/1‑93 (MQ=255)

                                                   ccACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGcc                                      >  1:6015069/1‑93 (MQ=255)

                                                    cACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCa                                     <  2:1535225/93‑1 (MQ=255)

                                                       gTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaa                                  >  1:7433462/1‑93 (MQ=255)

                                                        tATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaaa                                 >  2:8056621/1‑93 (MQ=255)

                                                        tATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaaa                                 <  1:7612287/93‑1 (MQ=255)

                                                        tATCAACCGTACCAAACAGCAATTCAACCCGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaaa                                 >  1:6993236/1‑93 (MQ=255)

                                                         aTCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaaaa                                >  2:7225886/1‑93 (MQ=255)

                                                         aTCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaaaa                                <  2:6037153/93‑1 (MQ=255)

                                                          tCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAg                               >  2:8422527/1‑93 (MQ=255)

                                                            aaCCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTa                             >  1:7209706/1‑93 (MQ=255)

                                                             aCCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTAc                            <  1:6257328/93‑1 (MQ=255)

                                                             aCCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTAc                            >  2:6850615/1‑93 (MQ=255)

                                                             aCCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACCTCACTTTGGTGCCAGAAAAGTAc                            <  1:917874/93‑1 (MQ=255)

                                                              ccGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACg                           <  2:1381576/93‑1 (MQ=255)

                                                               cGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                          <  2:340816/93‑1 (MQ=255)

                                                               cGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACgg                          <  2:7792512/93‑1 (MQ=255)

                                                                gTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACggg                         <  2:6935336/93‑1 (MQ=255)

                                                                gTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACggg                         <  2:8702664/93‑1 (MQ=255)

                                                                 tACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCTAGTAAAGTACgggc                        >  1:1748458/1‑92 (MQ=255)

                                                                 tACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGt                        >  2:829034/1‑93 (MQ=255)

                                                                  aCCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTc                       <  2:1956974/93‑1 (MQ=255)

                                                                     aaaCAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAAt                    >  1:4892446/1‑93 (MQ=255)

                                                                      aaCAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATc                   <  1:4318294/93‑1 (MQ=255)

                                                                      aaCAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATc                   >  2:3562145/1‑93 (MQ=255)

                                                                        cAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCAc                 <  1:5837379/93‑1 (MQ=255)

                                                                          gCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCa               <  2:2164705/93‑1 (MQ=255)

                                                                             aTTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAc            <  2:6744819/93‑1 (MQ=255)

                                                                             aTTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAc            <  1:6109261/93‑1 (MQ=255)

                                                                              ttCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGAcc           <  2:202859/93‑1 (MQ=255)

                                                                                caccacGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCgg         <  1:3439428/93‑1 (MQ=255)

                                                                                   cacGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCGCTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTAt      >  1:8601014/1‑93 (MQ=255)

                                                                                    acgacgCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTAtt     >  2:5177350/1‑93 (MQ=255)

                                                                                    acgacgCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTAtt     >  1:704904/1‑93 (MQ=255)

                                                                                    acgacgCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTAtt     <  2:5733274/93‑1 (MQ=255)

                                                                                      gacgCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgc   <  2:5049322/93‑1 (MQ=255)

                                                                                      gacgCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgc   <  1:38211/93‑1 (MQ=255)

                                                                                       acgCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgct  >  1:2545327/1‑93 (MQ=255)

GTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCT  >  NC_000913/4114245‑4114424

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