BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,387,576	T→G	K205Q (AAA→CAA)	aaeB ←	p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,387,576	0	T	G	98.6%	506.2 / ‑4.1	145	K205Q (AAA→CAA)	aaeB	p‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand): ref base T (1/1); new base G (71/72); total (72/73)
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 146 total aligned reads displayed.

CGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTT  >  NC_000913/3387479‑3387671

cGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCGGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGAt                                                                                               >  2:135652/1‑100 (MQ=255)

cGTCGATTGGCCCGCGCCCAGCGGGAAGATTACATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCCGGCTTGAt                                                                                               >  1:2127584/1‑100 (MQ=255)

 gTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATc                                                                                              >  2:12272/1‑100 (MQ=255)

 gTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATc                                                                                              <  1:1081346/100‑1 (MQ=255)

 gTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATc                                                                                              <  1:1360427/100‑1 (MQ=255)

 gTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATc                                                                                              <  1:1850053/100‑1 (MQ=255)

  tCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCg                                                                                             <  2:190401/100‑1 (MQ=255)

     aTTGGCCCGCGCCCAGCGGGAAGATTCCCTATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACa                                                                                          <  2:1547912/100‑1 (MQ=255)

     aTTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCATCAGGTCGCCCCAGGCTTGATCGACa                                                                                          <  1:3649540/100‑1 (MQ=255)

     aTTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACa                                                                                          <  1:1421272/100‑1 (MQ=255)

     aTTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTGGTAGCGCCGTGGTGCGTCGCAACAGGTCGCCCCAGGCTTGATCGACa                                                                                          <  1:3646568/100‑1 (MQ=255)

        ggCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACt                                                                                       >  1:73944/1‑100 (MQ=255)

        ggCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACt                                                                                       <  1:3578421/100‑1 (MQ=255)

        ggCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACt                                                                                       >  1:753681/1‑100 (MQ=255)

        ggCCCGCCCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGCCAACt                                                                                       >  2:45666/1‑100 (MQ=255)

         gCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACtt                                                                                      >  2:3424110/1‑100 (MQ=255)

           ccGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCa                                                                                    <  2:1395654/100‑1 (MQ=255)

              cgcCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCa                                                                                 >  1:2007122/1‑100 (MQ=255)

                cccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATc                                                                               >  1:4051265/1‑100 (MQ=255)

                 ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTTGCCCCAGGCTTGATCGACAACTTCACCATCg                                                                              >  2:1356202/1‑100 (MQ=255)

                  cAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGc                                                                             >  1:3314305/1‑100 (MQ=255)

                    gCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCa                                                                           >  2:1100377/1‑100 (MQ=255)

                    gCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCa                                                                           >  2:2404931/1‑100 (MQ=255)

                    gCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCa                                                                           <  1:2320966/100‑1 (MQ=255)

                    gCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCa                                                                           >  1:1867100/1‑100 (MQ=255)

                     cGGGAAGATTCAATATACAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCAt                                                                          >  2:3878190/1‑100 (MQ=255)

                      gggAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATg                                                                         <  2:598875/100‑1 (MQ=255)

                       ggAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGc                                                                        >  1:2092661/1‑100 (MQ=255)

                        gAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCt                                                                       <  2:2717417/100‑1 (MQ=255)

                         aaGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCtt                                                                      <  2:2969613/100‑1 (MQ=255)

                          aGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTg                                                                     <  1:2940182/100‑1 (MQ=255)

                            aTTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGAt                                                                   <  1:2378804/100‑1 (MQ=255)

                              tCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATAc                                                                 >  1:1124854/1‑100 (MQ=255)

                               ccATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACa                                                                >  1:154942/1‑100 (MQ=255)

                                cATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGCCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACag                                                               >  2:922558/1‑100 (MQ=255)

                                      cAGGTTGCTGCGCATGCCTTGTAGCGTCGTGGTGCGTCGAACCAGCTCGCCCCAGGCTTGATCGACAACTACACCGTCGCCATGCTTGATACAGAGTTGc                                                         <  1:3811601/100‑1 (MQ=255)

                                      cAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGc                                                         >  1:3142233/1‑100 (MQ=255)

                                       aGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCATCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCa                                                        <  2:3961133/100‑1 (MQ=255)

                                        ggTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCAt                                                       >  1:2200941/1‑100 (MQ=255)

                                         gTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCAtt                                                      >  1:3683175/1‑100 (MQ=255)

                                          ttGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTa                                                     <  1:2194553/100‑1 (MQ=255)

                                           tgctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTaa                                                    >  2:3560278/1‑100 (MQ=255)

                                            gctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAt                                                   <  2:786778/100‑1 (MQ=255)

                                             ctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAtt                                                  >  1:2674476/1‑100 (MQ=255)

                                             ctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAtt                                                  <  2:1288917/100‑1 (MQ=255)

                                             ctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAtt                                                  <  1:2796183/100‑1 (MQ=255)

                                             ctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAtt                                                  <  1:2796324/100‑1 (MQ=255)

                                             ctgcGCATGCCTTGAAGCGCCGTGGTGCGTCGAACCAGGTCGCACCACGCTTGATCGACAACTTCACCATCGCCATGCTTGATACGGAGTTGCATTAAtt                                                  <  2:2761429/100‑1 (MQ=255)

                                              tgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGCTACAGAGTTGCATTAATTg                                                 >  2:2311662/1‑100 (MQ=255)

                                              tgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTg                                                 <  2:911287/100‑1 (MQ=255)

                                                   aTGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATAtt                                            >  2:2152991/1‑100 (MQ=255)

                                                     gCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgc                                          >  1:858245/1‑100 (MQ=255)

                                                      ccgtGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGTCTTGATCGACAACTTTACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgcg                                         <  2:2245238/97‑1 (MQ=255)

                                                      ccTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgcg                                         >  1:2874475/1‑100 (MQ=255)

                                                      ccTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgcg                                         >  1:1897750/1‑100 (MQ=255)

                                                      ccTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCGTGATCGAAAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATGTTgcg                                         >  1:3872734/1‑100 (MQ=255)

                                                       cTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgcgc                                        <  1:2849884/100‑1 (MQ=255)

                                                       cTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgcgc                                        >  2:254643/1‑100 (MQ=255)

                                                        ttGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGcgcg                                       <  1:3452719/100‑1 (MQ=255)

                                                          ntAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGAc                                     >  1:3260467/2‑100 (MQ=255)

                                                             gcgcCGTGGTGAGTCGCACCAGGTCGCCCCAGGCTTGCTAGATAACTTCACCAACGACATGCTTGACACAGAGTTGTTTTTATTGATATTGCGCGACCAg                                  <  2:3023632/100‑1 (MQ=255)

                                                              cgcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGc                                 <  1:1521503/100‑1 (MQ=255)

                                                               gcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCa                                >  2:1771220/1‑100 (MQ=255)

                                                               gcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCa                                >  2:2051310/1‑100 (MQ=255)

                                                               gacctggtgCGACGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCa                                >  1:581604/5‑100 (MQ=255)

                                                               gacctggtgCGACGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATAGCCCTGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCa                                >  1:922409/5‑100 (MQ=255)

                                                                acctggtgCGACGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaa                               >  1:1909877/4‑100 (MQ=255)

                                                                  gtggtgCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAAc                             <  1:387942/100‑1 (MQ=255)

                                                                  gtggtgCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAAc                             <  1:3812546/100‑1 (MQ=255)

                                                                  gtggtgCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAAc                             <  2:1585940/100‑1 (MQ=255)

                                                                  gtggtgCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAAc                             <  2:1967994/100‑1 (MQ=255)

                                                                     gtgCGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACttt                          <  1:1501658/100‑1 (MQ=255)

                                                                      tgCGACGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTc                         >  1:3164235/1‑100 (MQ=255)

                                                                        cGTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCa                       <  1:1247699/100‑1 (MQ=255)

                                                                        cGTCGCACCAGGTCGCCACAGGCTAGATCGACAACTTCACCATCGCCATGCTTGATAAAGAGTTGCATTAATTGATATTGCGAGACCAGCAAACTTTCCa                       >  2:1763357/1‑100 (MQ=255)

                                                                         gTCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAg                      >  1:3074317/1‑100 (MQ=255)

                                                                          tCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGc                     <  2:3092703/100‑1 (MQ=255)

                                                                          tCGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGAAATTGCGCGACCAGCAAACTTTCCAGc                     <  2:2951653/100‑1 (MQ=255)

                                                                           cGCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGct                    >  1:2777452/1‑100 (MQ=255)

                                                                            gCACCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGctc                   >  2:1074587/1‑100 (MQ=255)

                                                                              aCCAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCtctc                 <  2:583242/100‑1 (MQ=255)

                                                                                cAGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGa               >  1:641792/1‑100 (MQ=255)

                                                                                 aGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGAt              <  1:2469912/100‑1 (MQ=255)

                                                                                 aGGTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGAGACAGAGTTGCATTAATTGATATTGCGCGACCAGCGAACTTTCCAGCTCTCGAt              <  2:374068/100‑1 (MQ=255)

                                                                                  ggTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATc             <  1:98229/100‑1 (MQ=255)

                                                                                  ggTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGGTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATc             >  1:1627644/1‑100 (MQ=255)

                                                                                   ggcgacccAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGATTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATcc            <  2:1917858/95‑1 (MQ=255)

                                                                                   gTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATcc            <  2:3664422/100‑1 (MQ=255)

                                                                                   gTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATcc            <  1:4070842/100‑1 (MQ=255)

                                                                                   gTCGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATcc            >  1:2760367/1‑100 (MQ=255)

                                                                                     cGCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCAc          <  1:717462/100‑1 (MQ=255)

                                                                                      gCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCAct         <  2:3702356/100‑1 (MQ=255)

                                                                                      gCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCAct         >  2:3915968/1‑100 (MQ=255)

                                                                                      gCCCCAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCAct         >  2:257339/1‑100 (MQ=255)

                                                                                       ccccAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCActt        >  1:2555508/1‑100 (MQ=255)

                                                                                       ccccAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTACGCGACCAGCAAACTTTCCAGCTCTCGATCCActt        <  2:2129442/100‑1 (MQ=255)

                                                                                          cAGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGTTCTCGATCCActtctt     >  1:3781341/1‑100 (MQ=255)

                                                                                           aGGCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCActtcttn    <  1:746162/100‑2 (MQ=255)

                                                                                            ggCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGt   >  1:3481621/1‑100 (MQ=255)

                                                                                             gCTTGATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGtt  <  2:656351/100‑1 (MQ=255)

CGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTT  >  NC_000913/3387479‑3387671

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