BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	318,413	G→A	R280R (CGC→CGT)	ykgC ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	318,413	0	G	A	95.9%	654.6 / 4.3	193	R280R (CGC→CGT)	ykgC	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base G (3/5); new base A (107/78); total (110/83)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 203 total aligned reads displayed.

GTAAGTAAATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATCCACCGCCA  >  NC_000913/318323‑318504

gtaagtaaATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCcacg                                                                                           >  2:830853/1‑93 (MQ=255)

 taagtaaATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCacgc                                                                                          >  1:1423396/1‑93 (MQ=255)

  aagtaaATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCt                                                                                         >  2:238205/1‑93 (MQ=255)

  aagtaaATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCt                                                                                         >  2:1954636/1‑93 (MQ=255)

     taaATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGt                                                                                      <  2:590198/93‑1 (MQ=255)

      aaaTTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGtt                                                                                     <  2:1708368/93‑1 (MQ=255)

      aaaTTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGtt                                                                                     <  1:2065841/93‑1 (MQ=255)

       aaTTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGttt                                                                                    >  1:1057953/1‑93 (MQ=255)

       aaTTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGttt                                                                                    <  1:1609413/93‑1 (MQ=255)

       aaTTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATAGCTTGTCAACGACAATTGCCCCACGCTCGttt                                                                                    >  1:2520815/1‑93 (MQ=255)

         ttGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTAc                                                                                  >  1:770815/1‑93 (MQ=255)

          tGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACt                                                                                 >  1:2457780/1‑93 (MQ=255)

            cAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGc                                                                               <  2:1569321/93‑1 (MQ=255)

            cAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGc                                                                               <  2:4643285/93‑1 (MQ=255)

            cAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGc                                                                               <  2:4528069/93‑1 (MQ=255)

             aGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCg                                                                              <  2:2543661/93‑1 (MQ=255)

              gCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGa                                                                             >  1:2989642/1‑93 (MQ=255)

              gCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGa                                                                             >  2:500989/1‑93 (MQ=255)

              gCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGa                                                                             <  1:2343958/93‑1 (MQ=255)

               cccgcccgTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGAt                                                                            <  1:216961/93‑1 (MQ=255)

               cccGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGAt                                                                            >  1:177222/1‑93 (MQ=255)

               cccGCCAGTAACATCTCCCATCGCCCGAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACACTTGCCCCACGCTCGTTTACTGCGAt                                                                            <  2:4217692/93‑1 (MQ=255)

                  gccgGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATAcc                                                                         <  1:2258785/93‑1 (MQ=255)

                  gccgGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATAcc                                                                         >  2:284740/1‑93 (MQ=255)

                    cgGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCgg                                                                       >  1:3126250/1‑93 (MQ=255)

                     ggTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGc                                                                      >  2:2003533/1‑93 (MQ=255)

                     ggTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGc                                                                      >  1:1152603/1‑93 (MQ=255)

                       tAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCAt                                                                    >  2:2302938/1‑93 (MQ=255)

                        taCATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGGAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCAtt                                                                   <  1:1591219/92‑1 (MQ=255)

                        aaCATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCAtt                                                                   >  1:3056530/1‑93 (MQ=255)

                         aCATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCAttt                                                                  >  1:3606919/1‑93 (MQ=255)

                         aCATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCAttt                                                                  >  1:1192421/1‑93 (MQ=255)

                         aCATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCAttt                                                                  >  2:1482103/1‑93 (MQ=255)

                         aCATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCAttt                                                                  >  1:1106815/1‑93 (MQ=255)

                           aTCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCTTCATTTTc                                                                >  2:4355958/1‑93 (MQ=255)

                             ctcCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTg                                                              <  1:3992281/93‑1 (MQ=255)

                             ctcCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTg                                                              >  2:3839337/1‑93 (MQ=255)

                               cccATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGa                                                            >  2:2638362/1‑93 (MQ=255)

                               cccATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGa                                                            <  1:2535273/93‑1 (MQ=255)

                                ccATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGAt                                                           >  1:1794096/1‑93 (MQ=255)

                                ccATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGAt                                                           >  2:811565/1‑93 (MQ=255)

                                ccATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGAt                                                           <  1:529248/93‑1 (MQ=255)

                                  aTCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATgg                                                         >  2:2689549/1‑92 (MQ=255)

                                  aTCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATga                                                         >  1:1025881/1‑92 (MQ=255)

                                  aTCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGt                                                         <  2:3926795/93‑1 (MQ=255)

                                  aTCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGt                                                         >  2:3769353/1‑93 (MQ=255)

                                  aTCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGt                                                         <  1:1923045/93‑1 (MQ=255)

                                   tCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTa                                                        >  1:916043/1‑93 (MQ=255)

                                   tCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTa                                                        <  1:3803556/93‑1 (MQ=255)

                                     gCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAAc                                                      >  2:2188043/1‑93 (MQ=255)

                                        cAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGaa                                                   <  2:3441333/93‑1 (MQ=255)

                                        cAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGaa                                                   >  2:3536688/1‑93 (MQ=255)

                                        cAAATATTGTCTGCGGTGGGATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGaa                                                   <  2:1541296/93‑1 (MQ=255)

                                         aaaTATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAg                                                  >  1:3832281/1‑93 (MQ=255)

                                         aaaTATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAg                                                  <  1:3969340/93‑1 (MQ=255)

                                          aaTATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGc                                                 <  2:4525180/93‑1 (MQ=255)

                                          aaTATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGc                                                 >  2:1333214/1‑93 (MQ=255)

                                            tatTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCgg                                               <  2:3258165/93‑1 (MQ=255)

                                            tatTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCgg                                               >  1:402405/1‑93 (MQ=255)

                                             atTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGt                                              >  2:4423329/1‑93 (MQ=255)

                                              ttGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTa                                             <  2:4020096/93‑1 (MQ=255)

                                              ttGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTa                                             >  2:1610251/1‑93 (MQ=255)

                                               tGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAg                                            <  1:264455/93‑1 (MQ=255)

                                                 tCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGcc                                          >  1:3761856/1‑93 (MQ=255)

                                                  cTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCg                                         <  2:3785000/93‑1 (MQ=255)

                                                    gCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGt                                       <  1:3419777/93‑1 (MQ=255)

                                                       gtggtATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTga                                    >  2:3820246/1‑93 (MQ=255)

                                                       gtggtATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTga                                    <  2:1385129/93‑1 (MQ=255)

                                                       gtggtATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTga                                    >  2:3091020/1‑93 (MQ=255)

                                                        tggtATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTgac                                   >  1:3728560/1‑93 (MQ=255)

                                                          gtATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTATCTGGATGTAACGAAGCGGTAGCCGGTTgacga                                 >  1:4588354/1‑93 (MQ=255)

                                                           tATGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTgacgac                                >  1:4184645/1‑93 (MQ=255)

                                                             tGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCg                              >  1:1873353/1‑93 (MQ=255)

                                                             tGTAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCg                              >  1:1357340/1‑93 (MQ=255)

                                                               tAATCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGa                            <  1:1369885/93‑1 (MQ=255)

                                                                aaTCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGaa                           <  2:4304163/93‑1 (MQ=255)

                                                                 aTCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAg                          >  2:1241295/1‑93 (MQ=255)

                                                                  tCGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGc                         >  2:322450/1‑93 (MQ=255)

                                                                   cGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCt                        <  2:482198/93‑1 (MQ=255)

                                                                   cGCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCt                        <  1:4043481/93‑1 (MQ=255)

                                                                    gCTTGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTa                       <  2:414712/93‑1 (MQ=255)

                                                                      ttGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTAtt                     <  1:3923987/93‑1 (MQ=255)

                                                                       tGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTa                    >  1:1269694/1‑93 (MQ=255)

                                                                       tGTCAACGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTa                    <  1:4200537/93‑1 (MQ=255)

                                                                          cAACGACAATTGCCCCACGCTCGTTTACTGCGATACCCGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACTGGACGCTATTACCa                 >  1:604391/1‑93 (MQ=255)

                                                                           aaCGACAGTTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAg                <  1:4572978/93‑1 (MQ=255)

                                                                           aaCGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAg                <  2:3159053/93‑1 (MQ=255)

                                                                           aaCGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAg                >  2:3619249/1‑93 (MQ=255)

                                                                            aCGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGt               <  1:819324/93‑1 (MQ=255)

                                                                             cGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTg              >  2:3492076/1‑93 (MQ=255)

                                                                             cGACAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTg              >  1:4672318/1‑93 (MQ=255)

                                                                               aCAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCa            >  1:4276378/1‑93 (MQ=255)

                                                                               aCAATTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCa            >  1:3658251/1‑93 (MQ=255)

                                                                                 aaTTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATc          >  2:3354739/1‑93 (MQ=255)

                                                                                  aTTGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATcc         >  1:1545289/1‑93 (MQ=255)

                                                                                  aTTGCCACACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATcc         >  1:2699331/1‑93 (MQ=255)

                                                                                    tGCCCCACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATCCAc       <  2:1367674/93‑1 (MQ=255)

                                                                                       cccACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGACGCTATTAACAGTGCATCCACCGc    >  2:3387922/1‑93 (MQ=255)

                                                                                       cccACGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATCCACCGc    <  1:2687366/93‑1 (MQ=255)

                                                                                         cgcgCTCGTTTACTGCTATACCGGCAGTATCTGGATGTAACGAAGCGGTAGTGGGTTGACGACCGGAAGCTATTAACAGTGCATCCACCGCCa  <  1:1660944/93‑1 (MQ=255)

GTAAGTAAATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATCCACCGCCA  >  NC_000913/318323‑318504

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