BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	262,642	T→C	E84E (GAA→GAG)	b0245 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	262,642	0	T	C	88.3%	716.9 / 79.8	282	E84E (GAA→GAG)	b0245	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (17/16); new base C (122/127); total (139/143)
Fisher's exact test for biased strand distribution p-value = 8.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 283 total aligned reads displayed.

TTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTGATACCG  >  NC_000913/262551‑262729

   aTCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCt                                                                                     >  1:2651092/1‑93 (MQ=255)

    tCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTg                                                                                    >  1:5624507/1‑93 (MQ=255)

    tCGTACTACGTTGTGACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTg                                                                                    >  2:5182805/1‑93 (MQ=255)

      gTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGcc                                                                                  >  2:1462336/1‑93 (MQ=255)

      gTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGcc                                                                                  >  2:6842298/1‑93 (MQ=255)

      gTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGcc                                                                                  >  1:2366602/1‑93 (MQ=255)

        actacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAg                                                                                >  2:2023800/1‑93 (MQ=255)

         ctacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGc                                                                               >  1:6701400/1‑93 (MQ=255)

          tacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCt                                                                              >  1:4064602/1‑93 (MQ=255)

          tacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCt                                                                              >  1:500850/1‑93 (MQ=255)

           acGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTa                                                                             <  1:2247380/93‑1 (MQ=255)

            cGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTaa                                                                            <  2:4949375/93‑1 (MQ=255)

            cGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTaa                                                                            <  1:19177/93‑1 (MQ=255)

             gttgttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTaaa                                                                           <  1:4764845/93‑1 (MQ=255)

             gttgttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTaaa                                                                           <  2:3101376/93‑1 (MQ=255)

             gttgttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTaaa                                                                           <  1:1487026/93‑1 (MQ=255)

              ttgttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAAt                                                                          >  1:1681535/1‑93 (MQ=255)

              ttgttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAAt                                                                          <  1:2889033/93‑1 (MQ=255)

              ttgttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAAt                                                                          >  1:422718/1‑93 (MQ=255)

              ttgttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAAt                                                                          >  2:3294230/1‑93 (MQ=255)

                gttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATcc                                                                        <  1:2004836/93‑1 (MQ=255)

                 ttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATccc                                                                       >  2:6906516/1‑93 (MQ=255)

                    cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTc                                                                    <  2:6305953/93‑1 (MQ=255)

                    cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTc                                                                    >  2:356753/1‑93 (MQ=255)

                    cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTc                                                                    <  2:2207662/93‑1 (MQ=255)

                     ggCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTcc                                                                   <  2:672192/93‑1 (MQ=255)

                      gCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCt                                                                  >  2:2891285/1‑93 (MQ=255)

                      gCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCt                                                                  <  2:2324285/93‑1 (MQ=255)

                      gCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCt                                                                  >  2:5624382/1‑93 (MQ=255)

                       cTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTg                                                                 <  1:3985375/93‑1 (MQ=255)

                        ttggCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGt                                                                <  1:901107/90‑1 (MQ=255)

                         ttGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTc                                                               <  2:6497130/93‑1 (MQ=255)

                           gCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGa                                                             <  2:5157501/93‑1 (MQ=255)

                           gCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGa                                                             >  1:390289/1‑93 (MQ=255)

                            ccGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGAt                                                            <  1:2138137/93‑1 (MQ=255)

                              gCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATAc                                                          >  1:1089650/1‑93 (MQ=255)

                                   aaGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAAcc                                                     <  1:6751497/93‑1 (MQ=255)

                                     gCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAg                                                   <  1:5786208/93‑1 (MQ=255)

                                     gCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAg                                                   >  2:3818640/1‑93 (MQ=255)

                                     gCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAg                                                   >  1:4206096/1‑93 (MQ=255)

                                     gCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAg                                                   >  2:1666841/1‑93 (MQ=255)

                                     gCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAg                                                   <  1:4618979/93‑1 (MQ=255)

                                        aGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTc                                                <  2:5972626/93‑1 (MQ=255)

                                        aGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTc                                                <  2:7189800/93‑1 (MQ=255)

                                        aGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTc                                                >  2:565208/1‑93 (MQ=255)

                                         gTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCg                                               >  2:2834803/1‑93 (MQ=255)

                                         gTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCg                                               >  2:7511666/1‑93 (MQ=255)

                                           tgcctgccGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTa                                             <  2:4708944/93‑1 (MQ=255)

                                           tgcctgccGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTa                                             <  1:6994130/93‑1 (MQ=255)

                                            gcctgccGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTAt                                            >  2:4191495/1‑93 (MQ=255)

                                             cctgccGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTAtt                                           <  2:2728687/93‑1 (MQ=255)

                                               tgccGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTAtttt                                         >  2:333925/1‑93 (MQ=255)

                                               tgccGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTAtttt                                         <  1:5488340/93‑1 (MQ=255)

                                                  cgcgCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCt                                      >  1:469368/1‑93 (MQ=255)

                                                   gcgcTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTa                                     <  2:1163919/93‑1 (MQ=255)

                                                    cgcTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTAc                                    >  2:369536/1‑93 (MQ=255)

                                                      cTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCa                                  <  1:4307854/93‑1 (MQ=255)

                                                        cGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGa                                <  1:1876567/93‑1 (MQ=255)

                                                         gCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGaa                               >  2:941980/1‑93 (MQ=255)

                                                         gCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGaa                               >  1:3470127/1‑93 (MQ=255)

                                                          cAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGaaa                              >  2:4236792/1‑93 (MQ=255)

                                                          cAGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGaaa                              >  2:3911871/1‑93 (MQ=255)

                                                           aGAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAACTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGaaaa                             >  2:6328079/1‑93 (MQ=255)

                                                            gAATGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAAt                            <  1:4809901/93‑1 (MQ=255)

                                                             aaTGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAAtt                           >  2:4576300/1‑93 (MQ=255)

                                                              aTGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTc                          >  1:7340119/1‑93 (MQ=255)

                                                              aTGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTc                          >  1:37827/1‑93 (MQ=255)

                                                              aTGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTc                          >  1:5262688/1‑93 (MQ=255)

                                                               tGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCa                         >  1:6478827/1‑93 (MQ=255)

                                                               tGTCTGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCa                         <  2:5521970/93‑1 (MQ=255)

                                                                gTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCAc                        <  2:1949340/93‑1 (MQ=255)

                                                                   tGCAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGc                     >  2:5270010/1‑93 (MQ=255)

                                                                     cAGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCAt                   <  2:4481452/93‑1 (MQ=255)

                                                                      aGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATc                  <  1:7455567/93‑1 (MQ=255)

                                                                      aGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATc                  <  1:3086312/93‑1 (MQ=255)

                                                                      aGCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATc                  <  1:4490464/93‑1 (MQ=255)

                                                                       gCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCg                 >  2:4303493/1‑93 (MQ=255)

                                                                       gCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCg                 <  1:5970740/93‑1 (MQ=255)

                                                                       gCCCGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCg                 >  2:5043624/1‑93 (MQ=255)

                                                                         ccGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTTTACCAGATAATTCACGGCATCGGc               <  2:7359817/93‑1 (MQ=255)

                                                                         ccGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGc               >  1:1000135/1‑93 (MQ=255)

                                                                         ccGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGc               >  1:1201798/1‑93 (MQ=255)

                                                                         ccGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGc               <  1:6238288/93‑1 (MQ=255)

                                                                         ccGGAGATAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGc               <  2:1480032/93‑1 (MQ=255)

                                                                         ccGGAGATAAGGAGATTGCTACTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGc               >  1:2973353/1‑93 (MQ=255)

                                                                           ggAGATAAAGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTa             >  1:5090936/1‑93 (MQ=255)

                                                                            gagaTAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAg            <  2:5874381/93‑1 (MQ=255)

                                                                              gaTAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAggg          <  2:275777/93‑1 (MQ=255)

                                                                              gaTAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAggg          >  2:253961/1‑93 (MQ=255)

                                                                               aTAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGt         >  2:4847317/1‑93 (MQ=255)

                                                                                tAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTg        >  1:4362722/1‑93 (MQ=255)

                                                                                tAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTg        >  1:3321896/1‑93 (MQ=255)

                                                                                tAAGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTg        >  1:2669844/1‑93 (MQ=255)

                                                                                  aGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTGAt      <  2:4407833/93‑1 (MQ=255)

                                                                                  aGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTGAt      <  2:2180204/93‑1 (MQ=255)

                                                                                  aGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTGAt      <  1:3389673/93‑1 (MQ=255)

                                                                                  aGGAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGGGAt      >  2:5402215/1‑93 (MQ=255)

                                                                                   ggAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTGATa     <  2:6012245/93‑1 (MQ=255)

                                                                                   ggAGATTGCTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTGATa     <  2:3447233/93‑1 (MQ=255)

                                                                                      gaTTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTGATACCg  >  2:1029982/1‑93 (MQ=255)

TTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTGATACCG  >  NC_000913/262551‑262729

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