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	99.6%	890.2 / NA	242	D473N (GAC→AAC)	mopA	GroEL
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (116/125); total (116/126)
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 267 total aligned reads displayed.

GCACTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCT  >  NC_000913/4369929‑4370108

gCACTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAAc                                                                                         <  1:768515/93‑1 (MQ=255)

gCACTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGATGCTAACACCGTTAAAGGCGGCAAc                                                                                         <  1:1962569/93‑1 (MQ=255)

   cTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGc                                                                                      >  1:834484/1‑93 (MQ=255)

   cTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGc                                                                                      <  1:2643515/93‑1 (MQ=255)

   cTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGc                                                                                      <  2:5229686/93‑1 (MQ=255)

       gTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACt                                                                                  >  1:2139795/1‑93 (MQ=255)

        tGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTa                                                                                 >  2:3722993/1‑93 (MQ=255)

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

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

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

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

               gAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTAc                                                                          >  1:2337852/1‑93 (MQ=255)

                aaGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACa                                                                         >  2:2639541/1‑93 (MQ=255)

                 aGCTCCGCTGCGTCAGAGCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACaa                                                                        <  2:1324399/93‑1 (MQ=255)

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

                    tCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACgc                                                                     >  1:2878894/1‑93 (MQ=255)

                      cGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACgcag                                                                   >  1:1033783/1‑93 (MQ=255)

                       gCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACgcagc                                                                  <  1:2564226/93‑1 (MQ=255)

                       gCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACgcagc                                                                  >  1:2586229/1‑93 (MQ=255)

                       gCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACgcagc                                                                  >  2:724749/1‑93 (MQ=255)

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

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

                              cAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaa                                                           <  2:3173815/93‑1 (MQ=255)

                                gATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGCTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaaga                                                         >  1:5198110/1‑93 (MQ=255)

                                 aTCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaagaa                                                        <  2:2809184/93‑1 (MQ=255)

                                 aTCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCgaagaa                                                        >  1:489401/1‑93 (MQ=255)

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

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

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

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

                                    gTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATAc                                                     <  1:2708409/93‑1 (MQ=255)

                                    gTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATAc                                                     <  2:3981636/93‑1 (MQ=255)

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

                                    gTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATAc                                                     <  2:5792726/93‑1 (MQ=255)

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

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

                                       ttGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGc                                                  >  1:2995137/1‑93 (MQ=255)

                                       ttGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGc                                                  >  1:2318646/1‑93 (MQ=255)

                                         gAACTGCGGCGAAGAACCGTCTGTTCTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCaa                                                <  1:6086755/93‑1 (MQ=255)

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

                                          aaCTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAAc                                               >  1:5930542/1‑93 (MQ=255)

                                             tGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATg                                            <  2:1213381/93‑1 (MQ=255)

                                              gcggcgAAGAAGCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGa                                           <  2:1004268/93‑1 (MQ=255)

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

                                                ggcgAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATc                                         <  2:4337198/93‑1 (MQ=255)

                                                  cgAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGa                                       <  1:5767569/93‑1 (MQ=255)

                                                  cgAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGa                                       <  2:3010074/93‑1 (MQ=255)

                                                   gaagaaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGAc                                      <  2:3027754/93‑1 (MQ=255)

                                                    aagaaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACa                                     <  2:6193317/93‑1 (MQ=255)

                                                    aagaaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACa                                     <  2:5439376/93‑1 (MQ=255)

                                                    aagaaCCGTCTGTTGTTGCTAACACCGCTAAAGGGGGCAAAGGCCACTCCGGTTACAACGCAGCAACCGAAGAATACGGCAAGCTGATCGACa                                     >  2:774219/1‑93 (MQ=255)

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

                                                      gaaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATg                                   <  1:588404/93‑1 (MQ=255)

                                                      gaaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATg                                   <  2:5873472/93‑1 (MQ=255)

                                                       aaCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATgg                                  >  1:1406794/1‑93 (MQ=255)

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

                                                        aCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATggg                                 >  1:5277707/1‑93 (MQ=255)

                                                        aCCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATggg                                 <  1:409841/93‑1 (MQ=255)

                                                         ccGTCTGTTGTTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGt                                >  1:1732352/1‑93 (MQ=255)

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

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

                                                              tgttgtTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCt                           >  1:5079706/1‑93 (MQ=255)

                                                              tgttgtTGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCt                           >  2:2743983/1‑93 (MQ=255)

                                                               gttgttGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTg                          <  2:2489686/93‑1 (MQ=255)

                                                               gttgttGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTg                          <  1:873134/93‑1 (MQ=255)

                                                               gttgttGCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTg                          <  2:4686482/93‑1 (MQ=255)

                                                                ttgttgCTAACACCGTTAAAGGCGGCAAGGGCACCTAGGGTTCAAACACAGCAACCGAAGAAAACGGCACCAGGTTCGCCATGGGAATACTgg                         >  2:3824414/1‑93 (MQ=255)

                                                                     gCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATcca                    <  1:4296662/93‑1 (MQ=255)

                                                                     gCTAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATcca                    <  1:5542908/93‑1 (MQ=255)

                                                                       tAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATccaac                  <  2:2053714/93‑1 (MQ=255)

                                                                       tAACACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATccaac                  <  2:198989/93‑1 (MQ=255)

                                                                        aaCACCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATccaacc                 <  1:3765545/93‑1 (MQ=255)

                                                                         acacCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATccaacca                >  1:1656248/1‑93 (MQ=255)

                                                                           acCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCaaa              <  1:5110966/93‑1 (MQ=255)

                                                                           acCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCaaa              <  2:5559835/93‑1 (MQ=255)

                                                                           acCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCaaa              >  2:1736139/1‑93 (MQ=255)

                                                                           acCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCaaa              <  1:4730966/93‑1 (MQ=255)

                                                                           acCGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCaaa              <  1:250537/93‑1 (MQ=255)

                                                                            ccGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAg             >  2:3589853/1‑93 (MQ=255)

                                                                            ccGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAg             >  1:1134951/1‑93 (MQ=255)

                                                                            ccGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATAGACATGGGTATCCTGGATCGAACCAAAg             >  1:1230135/1‑93 (MQ=255)

                                                                             cGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAga            >  1:1297975/1‑92 (MQ=255)

                                                                             cGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGt            <  2:4633897/93‑1 (MQ=255)

                                                                             cGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGt            <  2:5705899/93‑1 (MQ=255)

                                                                             cGTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGt            <  1:754797/93‑1 (MQ=255)

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

                                                                              gTTAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTa           <  2:6264670/93‑1 (MQ=255)

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

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

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

                                                                                tAAAGGCGGCAACGGCAACTGCGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAAc         <  2:3472875/93‑1 (MQ=255)

                                                                                tAAAGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAAc         <  2:3413589/93‑1 (MQ=255)

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

                                                                                  aaGGCGGCAACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAGCATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTc       <  2:5902855/93‑1 (MQ=255)

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

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

                                                                                     gcggcaACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGtt    <  1:4861281/93‑1 (MQ=255)

                                                                                     gcggcaACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGtt    <  1:70826/93‑1 (MQ=255)

                                                                                     gcggcaACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGtt    <  1:1995480/93‑1 (MQ=255)

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

GCACTGCGTGCAATGGAAGCTCCGCTGCGTCAGATCGTATTGAACTGCGGCGAAGAACCGTCTGTTGTTGCTAACACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCT  >  NC_000913/4369929‑4370108

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