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	81.2%	617.9 / 146.6	304	A156A (GCA→GCG)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base T (22/35); new base C (108/139); total (130/174)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 308 total aligned reads displayed.

GCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCG  >  NC_000913/4114243‑4114427

  gTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACATATCAACCGTACCAAACAGCAATTCACCACGACgcgc                                                                                            >  1:6250986/1‑93 (MQ=255)

   ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCa                                                                                           <  2:400302/93‑1 (MQ=255)

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

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

   ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCa                                                                                           <  1:3903319/93‑1 (MQ=255)

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

   ttGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCa                                                                                           <  1:79523/93‑1 (MQ=255)

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

        gtAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTc                                                                                      <  2:1659378/93‑1 (MQ=255)

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

         tAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTcg                                                                                     <  2:5434717/93‑1 (MQ=255)

          aaTCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTcgc                                                                                    >  1:5145147/1‑93 (MQ=255)

          aaTCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTcgc                                                                                    >  2:3039174/1‑93 (MQ=255)

           aTCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTcgcg                                                                                   >  2:1765900/1‑93 (MQ=255)

            tCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCga                                                                                  <  1:382875/93‑1 (MQ=255)

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

            tCGGTCACATGGACACGGCCATGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCga                                                                                  >  2:2849776/1‑93 (MQ=255)

               gTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGagag                                                                               <  2:2240138/93‑1 (MQ=255)

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

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

                tCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACATATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGc                                                                              <  2:3365233/93‑1 (MQ=255)

                 cacaTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGcc                                                                             <  1:2147973/93‑1 (MQ=255)

                  acaTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCt                                                                            <  1:2336303/93‑1 (MQ=255)

                  acaTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCt                                                                            >  1:4568834/1‑93 (MQ=255)

                  acaTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCt                                                                            >  2:2326137/1‑93 (MQ=255)

                   caTGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCtt                                                                           <  1:6442903/93‑1 (MQ=255)

                     tGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTcc                                                                         <  1:2025020/93‑1 (MQ=255)

                     tGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTcc                                                                         >  1:1489642/1‑93 (MQ=255)

                      ggACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCca                                                                        <  1:33069/93‑1 (MQ=255)

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

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

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

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

                              cccTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCg                                                                >  2:1457548/1‑93 (MQ=255)

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

                               ccTGAGTCATTTTCCAGATAAGCCACATATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGa                                                               >  1:4144787/1‑93 (MQ=255)

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

                                 tGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAgg                                                             >  1:1784787/1‑93 (MQ=255)

                                 tGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAgg                                                             >  2:6514770/1‑93 (MQ=255)

                                  gAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGa                                                            >  1:2255229/1‑93 (MQ=255)

                                  gAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGa                                                            >  1:995197/1‑93 (MQ=255)

                                  gAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGa                                                            <  2:4394925/93‑1 (MQ=255)

                                    gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc                                                          <  1:6664505/93‑1 (MQ=255)

                                    gTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATc                                                          <  2:2690476/93‑1 (MQ=255)

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

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

                                     tCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATcc                                                         <  2:6286283/93‑1 (MQ=255)

                                     tCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATcc                                                         >  1:1914797/1‑93 (MQ=255)

                                      cATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCca                                                        <  2:6658225/93‑1 (MQ=255)

                                      cATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCca                                                        <  2:6114479/93‑1 (MQ=255)

                                      cATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCca                                                        <  2:7158323/93‑1 (MQ=255)

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

                                         tttCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcactt                                                     <  2:6048385/93‑1 (MQ=255)

                                          ttCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttc                                                    <  2:2473278/93‑1 (MQ=255)

                                          ttCCAGATAAGCCACGTATCAACCGCACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttc                                                    <  1:5558673/93‑1 (MQ=255)

                                            ccAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcac                                                  <  2:377884/93‑1 (MQ=255)

                                            ccAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcac                                                  >  1:3666853/1‑93 (MQ=255)

                                            ccAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCcacttcac                                                  <  2:7376603/93‑1 (MQ=255)

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

                                                aTAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTg                                              <  2:3080438/93‑1 (MQ=255)

                                                 tAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTgg                                             >  2:514416/1‑93 (MQ=255)

                                                 tAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTgg                                             >  2:6705328/1‑93 (MQ=255)

                                                 tAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTgg                                             >  2:495953/1‑93 (MQ=255)

                                                 tAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTgg                                             <  2:1724550/93‑1 (MQ=255)

                                                   aGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTg                                           >  1:7223514/1‑93 (MQ=255)

                                                   aGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTg                                           >  2:6920034/1‑93 (MQ=255)

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

                                                      cacaTATCAACCGGACCAAACAGCAATTCACCACGACTCGCACGCTCGGGAGAGAATTCCACATGTTCGAGGTTCCACTTCACTTTGGGGCCa                                        >  1:1638111/1‑93 (MQ=255)

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

                                                       aCGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAg                                       >  1:5987224/1‑93 (MQ=255)

                                                         gTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGaa                                     <  1:2891502/93‑1 (MQ=255)

                                                         gTATCAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCTAGGATCCACTTCACTTTGGTGCCAGaa                                     >  2:70086/1‑93 (MQ=255)

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

                                                             cAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGt                                 >  1:2553793/1‑93 (MQ=255)

                                                             cAACCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGt                                 <  1:5033568/93‑1 (MQ=255)

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

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

                                                               aCCGTACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGCGCCAGAAAAGTAc                               >  1:7182848/1‑93 (MQ=255)

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

                                                                ccGTACCAAACAGCAATCCCCCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACg                              <  1:1748644/93‑1 (MQ=255)

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

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

                                                                   tACCAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGt                           <  2:796533/93‑1 (MQ=255)

                                                                    aCCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTc                          <  1:2391592/93‑1 (MQ=255)

                                                                      cAAACAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCaa                        <  2:5536988/93‑1 (MQ=255)

                                                                       aaaCAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAAt                       <  1:698920/93‑1 (MQ=255)

                                                                         aCAGCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCa                     <  2:6263126/93‑1 (MQ=255)

                                                                           aGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCAcc                   >  2:4254288/1‑93 (MQ=255)

                                                                           aGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCAcc                   <  1:2621214/93‑1 (MQ=255)

                                                                            gCAATTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCa                  >  1:752067/1‑93 (MQ=255)

                                                                              aaTTCACCACGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGa                <  2:5942628/93‑1 (MQ=255)

                                                                                  caccacGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAAAAAAGTACGGGTCAATCACCAGACCgg            <  1:5305313/93‑1 (MQ=255)

                                                                                     cacGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTAt         <  2:7427721/93‑1 (MQ=255)

                                                                                     cacGACGCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTAt         <  2:1384744/93‑1 (MQ=255)

                                                                                      acgacgTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTAtt        <  1:2458579/93‑1 (MQ=255)

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

                                                                                        gacgCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATACACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgc      <  2:203888/93‑1 (MQ=255)

                                                                                         acgTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgct     <  2:5071238/93‑1 (MQ=255)

                                                                                         acgCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgct     <  2:6606194/93‑1 (MQ=255)

                                                                                         acgCGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATtgct     <  2:5789978/93‑1 (MQ=255)

GCGTTGGTGTAATCGGTCACATGGACACGGCCCTGAGTCATTTTCCAGATAAGCCACGTATCAACCGTACCAAACAGCAATTCACCACGACGTGCACGCTCGCGAGAGCCTTCCACATGGTCGAGGATCCACTTCACTTTGGTGCCAGAAAAGTACGGGTCAATCACCAGACCGGTATTGCTGCG  >  NC_000913/4114243‑4114427

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