BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	631,938	C→T	A255T (GCA→ACA)	ybdH ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	631,938	0	C	T	99.6%	842.8 / ‑5.7	236	A255T (GCA→ACA)	ybdH	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (100/135); total (100/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

Only 100 of 248 total aligned reads displayed.

CACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTGCGTTAAC  >  NC_000913/631849‑632029

cACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTgtcgc                                                                                          >  1:7094017/1‑93 (MQ=255)

  ccAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTgtcgccc                                                                                        >  2:2458850/1‑92 (MQ=255)

   cAGAATTCCGTAGGCGTCTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCATACCGTTATGCACGGCATGAGCTGTCGCcac                                                                                       <  1:7229445/93‑1 (MQ=255)

    aGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCcaca                                                                                      <  1:1032272/93‑1 (MQ=255)

    aGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCcaca                                                                                      <  2:6925736/93‑1 (MQ=255)

     gAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCacac                                                                                     <  2:6149196/93‑1 (MQ=255)

      aaTTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACAcg                                                                                    <  2:4748589/93‑1 (MQ=255)

           cggaGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTaa                                                                               <  1:6255296/90‑1 (MQ=255)

           cGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTaa                                                                               <  1:2045552/93‑1 (MQ=255)

            gTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTaaa                                                                              <  1:6284947/93‑1 (MQ=255)

             tAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTaaaa                                                                             >  1:2360298/1‑93 (MQ=255)

             tAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTaaaa                                                                             >  2:2943820/1‑93 (MQ=255)

             tAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTaaaa                                                                             <  1:3132073/93‑1 (MQ=255)

              aGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAAc                                                                            >  2:116034/1‑93 (MQ=255)

              aGGCAACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAAc                                                                            <  1:1721230/93‑1 (MQ=255)

               ggCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACg                                                                           <  1:1000002/93‑1 (MQ=255)

                gCGACTTTGGTGCCGTGGAGAAAGTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGa                                                                          <  2:5652491/93‑1 (MQ=255)

                gCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGa                                                                          <  2:2135739/93‑1 (MQ=255)

                 ggACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGAt                                                                         <  1:2806504/92‑1 (MQ=255)

                 ggACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGAt                                                                         <  1:5582034/92‑1 (MQ=255)

                 cGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGAt                                                                         <  1:6355421/93‑1 (MQ=255)

                 cGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGAt                                                                         <  1:689066/93‑1 (MQ=255)

                  gACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATc                                                                        <  1:7753014/93‑1 (MQ=255)

                  gACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATc                                                                        <  1:7928145/93‑1 (MQ=255)

                  gACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATc                                                                        >  2:374737/1‑93 (MQ=255)

                  gACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATc                                                                        <  2:69889/93‑1 (MQ=255)

                   aCTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCg                                                                       >  1:5487506/1‑93 (MQ=255)

                   aCTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCg                                                                       >  2:3778376/1‑93 (MQ=255)

                   aCTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCg                                                                       >  2:299884/1‑93 (MQ=255)

                     tttGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGcc                                                                     >  1:547743/1‑93 (MQ=255)

                     tttGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGcc                                                                     <  1:4450897/93‑1 (MQ=255)

                     tttGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGcc                                                                     >  2:1699849/1‑93 (MQ=255)

                      ttGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGccc                                                                    <  2:3721186/93‑1 (MQ=255)

                      ttGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGccc                                                                    <  2:5872480/93‑1 (MQ=255)

                       tGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCa                                                                   <  1:2568274/93‑1 (MQ=255)

                          tGCCGTGGAGAATCTACTCAGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGac                                                                <  2:7040146/93‑1 (MQ=255)

                           gCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATTAGCTGTCGCCACACGCGTAAAACGATCGCCCAGacc                                                               >  1:4188028/1‑93 (MQ=255)

                             cGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGaccac                                                             >  2:7407964/1‑93 (MQ=255)

                             cGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGaccac                                                             >  1:900004/1‑93 (MQ=255)

                              gTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGaccacc                                                            <  2:5367648/93‑1 (MQ=255)

                              gTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGaccacc                                                            <  2:2202431/93‑1 (MQ=255)

                              gTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGaccacc                                                            >  2:5228759/1‑93 (MQ=255)

                              gTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGaccacc                                                            >  2:751106/1‑93 (MQ=255)

                               tGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGAccacca                                                           >  1:1399368/1‑93 (MQ=255)

                               tGGAGAAACTTCTCGGTTTGCGGCAGAACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGCTCGCCCAGAccacca                                                           >  1:6330726/1‑93 (MQ=255)

                                 gagaAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCaccaac                                                         >  2:7921032/1‑93 (MQ=255)

                                 gagaAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCaccaac                                                         <  1:5698250/93‑1 (MQ=255)

                                   gaAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCaccaacca                                                       >  2:1877054/1‑93 (MQ=255)

                                   gaAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCaccaacca                                                       <  2:3845483/93‑1 (MQ=255)

                                     aaCTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATc                                                     <  1:4820705/93‑1 (MQ=255)

                                      aCTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATcc                                                    >  1:5099153/1‑93 (MQ=255)

                                      aCTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATcc                                                    <  2:3445242/93‑1 (MQ=255)

                                      aCTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATcc                                                    <  2:3148726/93‑1 (MQ=255)

                                        ttCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGGCACACGCGTAAAACGATCGCCCAGACCACCAACCATCcca                                                  <  1:2968470/93‑1 (MQ=255)

                                         tctcGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCccac                                                 <  1:6170500/93‑1 (MQ=255)

                                         tctcGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCccac                                                 >  1:5273194/1‑93 (MQ=255)

                                         tctcGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCccac                                                 >  2:5479729/1‑93 (MQ=255)

                                              gTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCccaccacca                                            >  1:3597060/1‑93 (MQ=255)

                                                ttGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGc                                          <  2:562050/93‑1 (MQ=255)

                                                    ggCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCaata                                      <  1:3205758/93‑1 (MQ=255)

                                                     gcagcaCGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCaataa                                     >  1:5036432/1‑93 (MQ=255)

                                                     gcagcaCGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCaataa                                     <  2:4736606/93‑1 (MQ=255)

                                                      cagcaCGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCaataat                                    <  2:7077996/93‑1 (MQ=255)

                                                       agcaCGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAataata                                   <  2:2848650/93‑1 (MQ=255)

                                                        gcaCGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAg                                  >  2:1643938/1‑93 (MQ=255)

                                                        gcaCGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAg                                  >  2:7062713/1‑93 (MQ=255)

                                                         caCGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGc                                 >  2:183086/1‑93 (MQ=255)

                                                           cGGTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCAt                               >  1:7523704/1‑93 (MQ=255)

                                                             gTCAGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAATCATCCCACCACCAGCAATAATAGCATcc                             >  1:1855427/1‑93 (MQ=255)

                                                              tCAGACCGTGAGGCACGGCATGAGCTGTTGCCACACGCGGAATACGAGCGCTCAGACCACCAACCATCCCACCATCATCAACAATAGCAAcca                            >  1:7024061/1‑93 (MQ=255)

                                                                aGACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATccacc                          >  2:1550849/1‑93 (MQ=255)

                                                                 gACCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATccacca                         <  2:5540537/93‑1 (MQ=255)

                                                                  gccGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCcaccac                        <  2:626480/92‑1 (MQ=255)

                                                                  aCCGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCcaccac                        >  1:660165/1‑93 (MQ=255)

                                                                   ccGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACcaca                       <  1:6942165/93‑1 (MQ=255)

                                                                   ccGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACcaca                       <  1:4046561/93‑1 (MQ=255)

                                                                    cGTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACAt                      >  2:5367697/1‑93 (MQ=255)

                                                                     gTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAATGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATc                     <  1:3532924/93‑1 (MQ=255)

                                                                     gTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATc                     <  1:1870641/93‑1 (MQ=255)

                                                                     gTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATc                     <  1:7275538/93‑1 (MQ=255)

                                                                     gTTATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATc                     <  2:5025519/93‑1 (MQ=255)

                                                                       tATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGc                   <  2:423086/93‑1 (MQ=255)

                                                                       tATGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGc                   >  1:4132444/1‑93 (MQ=255)

                                                                         tGCACGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCaa                 >  2:5195258/1‑93 (MQ=255)

                                                                            aCGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAAt              >  2:1946756/1‑93 (MQ=255)

                                                                            aCGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAAt              <  1:646798/93‑1 (MQ=255)

                                                                             cGGCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATg             >  1:2898998/1‑93 (MQ=255)

                                                                              ggCATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGa            <  2:2804765/93‑1 (MQ=255)

                                                                                cATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGAtt          <  2:6705415/93‑1 (MQ=255)

                                                                                cATGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGAtt          <  1:3981290/93‑1 (MQ=255)

                                                                                 aTGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTg         <  2:647787/93‑1 (MQ=255)

                                                                                 aTGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTg         <  2:1379081/93‑1 (MQ=255)

                                                                                 aTGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTg         <  1:1155691/93‑1 (MQ=255)

                                                                                  tGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTGc        >  2:6375924/1‑93 (MQ=255)

                                                                                  tGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTGc        >  2:3351268/1‑93 (MQ=255)

                                                                                  tGAGCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTGc        <  2:2548836/93‑1 (MQ=255)

                                                                                  tGAGCTGCCGCCACACGCGTAAACCGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTGc        <  2:4299180/93‑1 (MQ=255)

                                                                                     gCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTGCGtt     >  1:6092457/1‑93 (MQ=255)

                                                                                     gCTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTGCGtt     <  1:759988/93‑1 (MQ=255)

                                                                                      cTGTCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTGCGTTa    <  2:321400/93‑1 (MQ=255)

CACCAGAATTCCGTAGGCGACTTTGGTGCCGTGGAGAAACTTCTCGGTTTGCGGCAGCACGGTCAGACCGTTATGCACGGCATGAGCTGCCGCCACACGCGTAAAACGATCGCCCAGACCACCAACCATCCCACCACCAGCAATAATAGCATCCACCACATCGCAAAATGATTGCGTTAAC  >  NC_000913/631849‑632029

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