BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,370,019	G→A	D473N (GAC→AAC)	mopA →	GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,370,019	0	G	A	93.9%	1433.5 / 37.5	443	D473N (GAC→AAC)	mopA	GroEL
Reads supporting (aligned to +/- strand): ref base G (9/18); new base A (186/230); total (195/248)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 477 total aligned reads displayed.

GCACTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCT  >  NC_000913/4369929‑4370108

 cACTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACg                                                                                        >  2:9703170/1‑93 (MQ=255)

 cACTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACg                                                                                        <  2:2159129/93‑1 (MQ=255)

    tGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGGAACGGCa                                                                                     >  1:7297337/1‑93 (MQ=255)

    tGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCa                                                                                     >  2:4506398/1‑93 (MQ=255)

    tGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCa                                                                                     >  2:3212364/1‑93 (MQ=255)

         gCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTAc                                                                                <  2:80588/93‑1 (MQ=255)

         gCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTAc                                                                                >  2:6658082/1‑93 (MQ=255)

         gCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTAc                                                                                <  1:7602889/93‑1 (MQ=255)

         gCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTAc                                                                                <  2:9096255/93‑1 (MQ=255)

          cAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACg                                                                               <  2:3792631/93‑1 (MQ=255)

           aaTGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACgg                                                                              >  1:2319312/1‑93 (MQ=255)

           aaTGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACgg                                                                              <  2:2948370/93‑1 (MQ=255)

           aaTGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGGTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACgg                                                                              <  1:391007/93‑1 (MQ=255)

            aTGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGt                                                                             <  1:9553397/93‑1 (MQ=255)

            aTGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGt                                                                             >  1:4867531/1‑93 (MQ=255)

               gAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTaa                                                                          >  1:6405515/1‑92 (MQ=255)

                aaGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGGTACa                                                                         >  2:7252984/1‑93 (MQ=255)

                 aGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACaa                                                                        <  1:2961203/93‑1 (MQ=255)

                 aGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACaa                                                                        <  1:8744271/93‑1 (MQ=255)

                 aGCGCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCGACGGCAACTACGGTTACaa                                                                        <  1:6224008/93‑1 (MQ=255)

                  gCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAAc                                                                       <  1:4638564/93‑1 (MQ=255)

                    tCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACgc                                                                     <  2:9451983/93‑1 (MQ=255)

                     ccGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACgca                                                                    >  2:6103788/1‑93 (MQ=255)

                      cGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAA‑GCCGCGACGGCAACTTCGGTTACAAcgccgc                                                                  >  2:4921662/1‑93 (MQ=255)

                      cGCTGCGTCAGAGCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTGACAACgcag                                                                   >  2:8979000/1‑93 (MQ=255)

                       gCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAATgcagc                                                                  >  2:2001352/1‑93 (MQ=255)

                         tGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCaa                                                                <  1:4371207/93‑1 (MQ=255)

                          gCGTCAGATCGTATTGAACTGCGGCGAAGACCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAAc                                                               <  2:3421277/93‑1 (MQ=255)

                          gCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAAc                                                               >  1:9927086/1‑93 (MQ=255)

                             tCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCga                                                            >  2:3042562/1‑93 (MQ=255)

                               aGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaag                                                          >  2:6377117/1‑93 (MQ=255)

                               aGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaag                                                          >  1:6634248/1‑93 (MQ=255)

                                gATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaaga                                                         >  1:5817554/1‑93 (MQ=255)

                                gATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaaga                                                         <  2:8173567/93‑1 (MQ=255)

                                gATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaaga                                                         <  1:5059309/93‑1 (MQ=255)

                                gATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaaga                                                         >  2:6522479/1‑93 (MQ=255)

                                  tCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGGAt                                                       >  2:1458393/1‑93 (MQ=255)

                                  tCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAAt                                                       <  2:8716718/93‑1 (MQ=255)

                                  tCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAAt                                                       <  2:5843389/93‑1 (MQ=255)

                                   cGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATa                                                      <  2:1785544/93‑1 (MQ=255)

                                   cGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATa                                                      <  1:4014684/93‑1 (MQ=255)

                                   cGTATTGAACTGCAGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATa                                                      <  1:7125413/93‑1 (MQ=255)

                                   cGTACTGACCTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATa                                                      <  2:3840/93‑1 (MQ=255)

                                    gTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATAc                                                     >  1:5150956/1‑93 (MQ=255)

                                    gTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCAGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATAc                                                     <  1:1980099/93‑1 (MQ=255)

                                     tATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCACCCGACGAATACg                                                    >  1:2703823/1‑93 (MQ=255)

                                     tATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACg                                                    >  1:93502/1‑93 (MQ=255)

                                      aTTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACgg                                                   >  1:6948334/1‑93 (MQ=255)

                                      aTTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACgg                                                   <  1:7877252/93‑1 (MQ=255)

                                       ttGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGc                                                  <  2:6250340/93‑1 (MQ=255)

                                        tGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCa                                                 >  2:4171114/1‑93 (MQ=255)

                                          aaCTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAAc                                               <  2:8850153/93‑1 (MQ=255)

                                          aaCTGCGGCGAAGAACCGTCTGTTGTTGCTAACAACGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAAc                                               <  2:7165397/93‑1 (MQ=255)

                                           aCTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAGCGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACa                                              <  2:2774548/93‑1 (MQ=255)

                                           aCTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACa                                              >  2:9188946/1‑93 (MQ=255)

                                              gcggcgAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGa                                           <  1:7498217/93‑1 (MQ=255)

                                                 gcgAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCg                                        >  1:10008092/1‑93 (MQ=255)

                                                   gaagaaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGAc                                      >  2:2582061/1‑93 (MQ=255)

                                                     agaaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACAt                                    <  2:4048085/93‑1 (MQ=255)

                                                      gaaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATg                                   >  2:10029769/1‑93 (MQ=255)

                                                       aaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATgg                                  <  2:9598131/93‑1 (MQ=255)

                                                         ccGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGCCATGGGt                                >  2:8531056/1‑93 (MQ=255)

                                                           gTCTGTTGTTGCTAACACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTAt                              >  1:2052180/1‑93 (MQ=255)

                                                            tCTGTTGTTGCTAACACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATc                             <  2:469992/93‑1 (MQ=255)

                                                            tCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATc                             >  1:3473348/1‑93 (MQ=255)

                                                            tCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACGACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATc                             <  1:8563901/93‑1 (MQ=255)

                                                             cTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATcc                            <  1:8852506/93‑1 (MQ=255)

                                                              tgttgtTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCt                           <  2:7493223/93‑1 (MQ=255)

                                                              tgttgtTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCt                           <  2:70642/93‑1 (MQ=255)

                                                              tgttgtTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATAGACATGGGTATccc                           >  1:4258393/1‑92 (MQ=255)

                                                                ttgttgCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTgg                         <  2:8017060/93‑1 (MQ=255)

                                                                    tgCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATcc                     <  1:9715240/93‑1 (MQ=255)

                                                                      cTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATccaa                   <  2:826488/93‑1 (MQ=255)

                                                                      cTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATccaa                   <  1:3377038/93‑1 (MQ=255)

                                                                      cTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATAAGGCAACATGATCGACATGGGCATCCTGGATccaa                   >  1:2174579/1‑93 (MQ=255)

                                                                        aaCGCCGTTAAAGGCGGCAACGGCAACGACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATccaacc                 <  1:1887036/93‑1 (MQ=255)

                                                                        aaCACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCAAAGAATACGGCAACATGATCGACATGGGTATCCTGGATccaacc                 >  2:644481/1‑93 (MQ=255)

                                                                          cacCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATccaaccaa               <  1:1524878/93‑1 (MQ=255)

                                                                            ccGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAg             <  1:6768605/93‑1 (MQ=255)

                                                                            ccGTGAAAGGCGGCAAGTGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAg             <  1:8683473/93‑1 (MQ=255)

                                                                             cGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGt            <  2:8456352/93‑1 (MQ=255)

                                                                              gTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTa           >  1:5529036/1‑93 (MQ=255)

                                                                              gTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTa           >  1:9009925/1‑93 (MQ=255)

                                                                               ttAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTaa          <  1:8406965/93‑1 (MQ=255)

                                                                               ttAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGaaa          >  1:3254673/1‑93 (MQ=255)

                                                                               ttAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTaa          >  2:1196954/1‑93 (MQ=255)

                                                                               ttAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTaa          >  1:5728611/1‑93 (MQ=255)

                                                                               ttAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTaa          <  1:1095636/93‑1 (MQ=255)

                                                                                tAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTAGATCCAACCAAAGTAAc         >  2:5314433/1‑93 (MQ=255)

                                                                                 aaaGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATTGGTATCCTGGATCCAACCAAAGTAACt        <  1:5969350/93‑1 (MQ=255)

                                                                                 aaaGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACt        <  2:6750101/93‑1 (MQ=255)

                                                                                  aaGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTc       <  1:6596112/93‑1 (MQ=255)

                                                                                  aaGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTc       <  1:2545210/93‑1 (MQ=255)

                                                                                   aGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATGCTGGATCCAACCACAGTAACTCg      >  2:5534799/1‑93 (MQ=255)

                                                                                   aGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACAAGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCg      <  1:5220329/93‑1 (MQ=255)

                                                                                    ggcggcAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGt     >  2:4276235/1‑93 (MQ=255)

                                                                                    ggcggcAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGt     <  2:2233915/93‑1 (MQ=255)

                                                                                    ggcggcAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGt     >  1:1206274/1‑93 (MQ=255)

                                                                                     gcggcaACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGtt    >  2:7841909/1‑93 (MQ=255)

                                                                                       ggcaacggcaacTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTtct  <  2:5701599/93‑1 (MQ=255)

GCACTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCT  >  NC_000913/4369929‑4370108

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 17 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`