BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,935,303	A→G	Q128Q (CAA→CAG)	rbsK →	ribokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,935,303	0	A	G	84.7%	725.5 / 115.0	314	Q128Q (CAA→CAG)	rbsK	ribokinase
Reads supporting (aligned to +/- strand): ref base A (20/28); new base G (144/122); total (164/150)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 315 total aligned reads displayed.

GCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAA  >  NC_000913/3935211‑3935393

 cgcTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAAc                                                                                           <  2:1523473/93‑1 (MQ=255)

  gcTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCg                                                                                          >  1:6863144/1‑93 (MQ=255)

  gcTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCg                                                                                          >  1:3978160/1‑93 (MQ=255)

  gcTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCg                                                                                          >  2:3032157/1‑93 (MQ=255)

   cTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGt                                                                                         >  2:2493981/1‑93 (MQ=255)

   cTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGt                                                                                         >  1:14248/1‑93 (MQ=255)

    tGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTg                                                                                        >  2:3420228/1‑93 (MQ=255)

     gATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGa                                                                                       <  1:5592533/93‑1 (MQ=255)

     gATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGa                                                                                       <  1:6439277/93‑1 (MQ=255)

     gATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGa                                                                                       <  2:4075112/93‑1 (MQ=255)

       tttttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGc                                                                                     >  2:2942273/1‑93 (MQ=255)

       tttttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGc                                                                                     >  2:2230240/1‑93 (MQ=255)

          ttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTa                                                                                  <  1:2806537/93‑1 (MQ=255)

          ttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTa                                                                                  <  2:2204198/93‑1 (MQ=255)

           tGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTAt                                                                                 <  2:7512102/93‑1 (MQ=255)

           tGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTAt                                                                                 <  2:170703/93‑1 (MQ=255)

              tAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGc                                                                              >  2:4370441/1‑93 (MQ=255)

               aaTGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGcc                                                                             <  2:3679005/93‑1 (MQ=255)

                aTGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCCGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCa                                                                            <  1:5844470/93‑1 (MQ=255)

                 tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCaa                                                                           <  1:2332528/93‑1 (MQ=255)

                 tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCaa                                                                           <  1:4649946/93‑1 (MQ=255)

                 tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCaa                                                                           <  1:3478219/93‑1 (MQ=255)

                 tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCaa                                                                           >  2:6837488/1‑93 (MQ=255)

                   gCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAAcg                                                                         >  2:5627870/1‑93 (MQ=255)

                    cGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAAcgc                                                                        >  2:5755915/1‑93 (MQ=255)

                      aaGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGt                                                                      >  2:2729261/1‑93 (MQ=255)

                      aaGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGt                                                                      <  2:4276524/93‑1 (MQ=255)

                      aaGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGt                                                                      >  2:6990533/1‑93 (MQ=255)

                       aGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTc                                                                     <  2:932183/93‑1 (MQ=255)

                        ggTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCa                                                                    <  2:662163/93‑1 (MQ=255)

                        ggTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCa                                                                    <  1:4501066/93‑1 (MQ=255)

                        ggTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCa                                                                    <  1:1755569/93‑1 (MQ=255)

                         gTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAg                                                                   <  1:2020019/93‑1 (MQ=255)

                             gaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCatt                                                               >  2:2913886/1‑93 (MQ=255)

                              aaTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCatta                                                              <  1:6318226/93‑1 (MQ=255)

                               aTGTCATCGGTATTCATGCCGGCGCTAGTGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCattat                                                             >  2:6421972/1‑93 (MQ=255)

                                tGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCattatt                                                            <  2:2613406/93‑1 (MQ=255)

                                tGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCattatt                                                            <  1:2710515/93‑1 (MQ=255)

                                   cATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAAt                                                         >  2:6081745/1‑93 (MQ=255)

                                   cATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAAt                                                         >  1:7158672/1‑93 (MQ=255)

                                    aTCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATg                                                        >  2:4682640/1‑93 (MQ=255)

                                    aTCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATg                                                        >  2:2314911/1‑93 (MQ=255)

                                     tCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGc                                                       >  2:3470457/1‑93 (MQ=255)

                                     tCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGc                                                       >  1:2143760/1‑93 (MQ=255)

                                     tCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGc                                                       <  1:5931928/93‑1 (MQ=255)

                                     tCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGc                                                       >  2:2068900/1‑93 (MQ=255)

                                     tCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGc                                                       <  1:5737606/93‑1 (MQ=255)

                                      cGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCa                                                      <  1:6978876/93‑1 (MQ=255)

                                       ggTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAg                                                     >  2:533163/1‑93 (MQ=255)

                                        gTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGc                                                    >  2:1976175/1‑93 (MQ=255)

                                        gTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGc                                                    <  2:6725149/93‑1 (MQ=255)

                                        gTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGc                                                    >  1:7316343/1‑93 (MQ=255)

                                        gTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGc                                                    <  1:4085251/93‑1 (MQ=255)

                                         tATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCt                                                   >  1:6749723/1‑93 (MQ=255)

                                          aTTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTg                                                  >  1:1036771/1‑93 (MQ=255)

                                           ttCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTgg                                                 <  1:2399454/93‑1 (MQ=255)

                                           ttCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTgg                                                 >  1:3534299/1‑93 (MQ=255)

                                           ttCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTgg                                                 >  2:7165474/1‑93 (MQ=255)

                                            tCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGa                                                >  2:5880690/1‑93 (MQ=255)

                                              aTGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAAt                                              >  1:6400030/1‑93 (MQ=255)

                                               tGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                             <  1:6598570/93‑1 (MQ=255)

                                               tGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTCCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                             >  1:5794563/1‑93 (MQ=255)

                                                 ccGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcac                                           >  1:4893714/1‑93 (MQ=255)

                                                 ccGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcac                                           <  2:3312882/93‑1 (MQ=255)

                                                    gcgcTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcaccac                                        <  1:2497612/93‑1 (MQ=255)

                                                    gcgcTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcaccac                                        <  1:4514600/93‑1 (MQ=255)

                                                    gcgcTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcaccac                                        <  1:1890181/93‑1 (MQ=255)

                                                     cgcTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACt                                       >  1:433111/1‑93 (MQ=255)

                                                      gcTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTc                                      <  2:6384596/93‑1 (MQ=255)

                                                      gcTAATGATGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTc                                      >  1:369196/1‑93 (MQ=255)

                                                        tAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGa                                    >  1:1353626/1‑93 (MQ=255)

                                                         aaTGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGaa                                   <  2:3997237/93‑1 (MQ=255)

                                                           tgctgcCCTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAg                                 <  2:4553474/93‑1 (MQ=255)

                                                                cccTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGa                            <  2:5895504/93‑1 (MQ=255)

                                                                 ccTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGAt                           <  2:1241139/93‑1 (MQ=255)

                                                                  cTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATg                          >  2:635889/1‑93 (MQ=255)

                                                                  cTTTCCCCGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATg                          >  2:70885/1‑93 (MQ=255)

                                                                     tCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCa                       <  2:191806/93‑1 (MQ=255)

                                                                      ccccGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAg                      >  1:5045775/1‑93 (MQ=255)

                                                                      ccccGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAg                      <  2:1044747/93‑1 (MQ=255)

                                                                      ccccGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGAGTGATGGCAg                      <  2:1477170/93‑1 (MQ=255)

                                                                       cccGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgc                     <  2:1794347/93‑1 (MQ=255)

                                                                       cccGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgc                     >  1:6138810/1‑93 (MQ=255)

                                                                        ccGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcg                    <  1:2893583/93‑1 (MQ=255)

                                                                        ccGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcg                    >  2:1128705/1‑93 (MQ=255)

                                                                        ccGGCGCTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcg                    >  2:6889650/1‑93 (MQ=255)

                                                                            cgcTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGa                >  2:5044104/1‑93 (MQ=255)

                                                                            cgcTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGa                <  2:355667/93‑1 (MQ=255)

                                                                             gcTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaa               >  2:7375050/1‑93 (MQ=255)

                                                                              cTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaa              >  1:1282182/1‑93 (MQ=255)

                                                                              cTGGTGGAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaa              >  1:745168/1‑93 (MQ=255)

                                                                               tggtggAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaaa             >  1:3104705/1‑93 (MQ=255)

                                                                               tggtggAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaaa             <  2:5181865/93‑1 (MQ=255)

                                                                                 gtggAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATc           >  1:6013875/1‑93 (MQ=255)

                                                                                  tggAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCg          >  1:3212832/1‑93 (MQ=255)

                                                                                   ggAAGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGc         <  1:828412/93‑1 (MQ=255)

                                                                                     aaGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGccc       <  1:7074253/93‑1 (MQ=255)

                                                                                     aaGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGccc       <  2:2438448/93‑1 (MQ=255)

                                                                                      aGCGCAGCGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCa      >  2:4957581/1‑93 (MQ=255)

                                                                                       gcgcAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCAt     >  1:546167/1‑93 (MQ=255)

GCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAA  >  NC_000913/3935211‑3935393

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