BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,604,832	A→G	R146R (CGT→CGC)	yjjT ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,832	0	A	G	99.3%	443.1 / ‑4.5	138	R146R (CGT→CGC)	yjjT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (58/79); total (58/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

Only 100 of 144 total aligned reads displayed.

ACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCC  >  NC_000913/4604744‑4604921

acacGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAggcgac                                                                                       <  2:171295/93‑1 (MQ=255)

acacGCCAGGCAGCGTTTTGACCGTCAGGCCATAGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAggcgac                                                                                       >  2:2675343/1‑93 (MQ=255)

 cacGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGAcc                                                                                      >  1:4434160/1‑93 (MQ=255)

 cacGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGAcc                                                                                      >  2:4476609/1‑93 (MQ=255)

  acGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCa                                                                                     >  2:669083/1‑93 (MQ=255)

  acGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCa                                                                                     >  2:3097193/1‑93 (MQ=255)

     ccagccagCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCaaaa                                                                                  <  2:4702221/93‑1 (MQ=255)

     ccAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCaaaa                                                                                  <  2:4702220/93‑1 (MQ=255)

     ccAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCaaaa                                                                                  <  1:60492/93‑1 (MQ=255)

        ggCAGCGTTTTGACCGTTAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATag                                                                               <  1:3338763/93‑1 (MQ=255)

        ggCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATag                                                                               <  1:1916249/93‑1 (MQ=255)

        ggCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATag                                                                               <  1:1756262/93‑1 (MQ=255)

        ggCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATag                                                                               <  1:3873384/93‑1 (MQ=255)

         gCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATaga                                                                              >  2:4339510/1‑93 (MQ=255)

         gCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATaga                                                                              <  1:3524562/93‑1 (MQ=255)

         gCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATaga                                                                              <  2:2423595/93‑1 (MQ=255)

         gCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATaga                                                                              >  2:2632780/1‑93 (MQ=255)

         gCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATaga                                                                              <  1:3317188/93‑1 (MQ=255)

          gaGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATagag                                                                             <  1:3732373/92‑1 (MQ=255)

          cAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATagag                                                                             <  1:4136693/93‑1 (MQ=255)

            gCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGc                                                                           <  2:881108/93‑1 (MQ=255)

            gCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGc                                                                           >  2:2056465/1‑93 (MQ=255)

              gTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCca                                                                         >  1:579505/1‑93 (MQ=255)

               ttttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcac                                                                        >  2:523578/1‑93 (MQ=255)

                tttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcaca                                                                       <  1:1758626/93‑1 (MQ=255)

                tttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcaca                                                                       <  2:2846630/93‑1 (MQ=255)

                tttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcaca                                                                       <  1:138168/93‑1 (MQ=255)

                    aCCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcga                                                                   <  1:4318640/93‑1 (MQ=255)

                       gTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgacga                                                                <  2:1412398/93‑1 (MQ=255)

                        tCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCCAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgacgcg                                                               >  2:4101339/1‑91 (MQ=255)

                        tCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAg                                                               <  1:2198434/93‑1 (MQ=255)

                           ggCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgcgc                                                            <  1:1470626/93‑1 (MQ=255)

                           ggCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgcgc                                                            >  1:663335/1‑93 (MQ=255)

                            gCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCt                                                           <  2:1921068/93‑1 (MQ=255)

                              cATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGt                                                         <  2:3516819/93‑1 (MQ=255)

                               aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        >  1:2468520/1‑93 (MQ=255)

                               aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        >  1:673508/1‑93 (MQ=255)

                                  gACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGat                                                     >  1:906819/1‑92 (MQ=255)

                                  gACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGAc                                                     <  2:4377938/93‑1 (MQ=255)

                                   aCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACt                                                    <  2:485836/93‑1 (MQ=255)

                                   aCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACt                                                    <  2:542284/93‑1 (MQ=255)

                                   aCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACt                                                    >  2:1112744/1‑93 (MQ=255)

                                    cGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACtt                                                   <  1:3587531/93‑1 (MQ=255)

                                    cGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACtt                                                   >  2:4477886/1‑93 (MQ=255)

                                      cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 <  2:1058009/93‑1 (MQ=255)

                                      cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 >  1:1200340/1‑93 (MQ=255)

                                      cTGTATGCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 <  1:2354861/93‑1 (MQ=255)

                                       tGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAt                                                <  2:4012842/93‑1 (MQ=255)

                                        gTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAtt                                               <  1:2878518/93‑1 (MQ=255)

                                        gTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAtt                                               <  1:1645866/93‑1 (MQ=255)

                                         tATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTc                                              <  2:1258992/93‑1 (MQ=255)

                                          aTTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCa                                             >  1:1583724/1‑93 (MQ=255)

                                             cGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACg                                          <  1:1508966/93‑1 (MQ=255)

                                               ccccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgc                                        <  2:1511473/93‑1 (MQ=255)

                                                 ccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgcgc                                      >  1:82712/1‑93 (MQ=255)

                                                  cAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCa                                     <  1:3142558/93‑1 (MQ=255)

                                                   aGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCAt                                    <  2:2446887/93‑1 (MQ=255)

                                                     aaTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATaa                                  >  1:3638159/1‑93 (MQ=255)

                                                      aTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAAt                                 >  1:82465/1‑93 (MQ=255)

                                                      aTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAAt                                 <  2:1457787/93‑1 (MQ=255)

                                                       tttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                >  1:2899180/1‑93 (MQ=255)

                                                        ttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCt                               <  1:2078366/93‑1 (MQ=255)

                                                        ttCTCCGCAGCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCt                               <  2:1599321/93‑1 (MQ=255)

                                                         tctcCGCATAATATAGCGGCTGTTTGGCTAGGCTAACAAAATAGAGTCTAAAGCCACGAGCGCGGGCTACTTTATTCAACGGCGCATAATCTg                              <  2:905166/93‑1 (MQ=255)

                                                          ctcCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGc                             <  1:4404578/93‑1 (MQ=255)

                                                           tcCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGcc                            <  1:119629/93‑1 (MQ=255)

                                                            ccGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCa                           <  1:3459935/93‑1 (MQ=255)

                                                              gCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGc                         >  2:693139/1‑93 (MQ=255)

                                                              gCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGc                         <  1:4644954/93‑1 (MQ=255)

                                                               cATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCa                        >  2:2231723/1‑93 (MQ=255)

                                                               cATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCa                        <  2:1232579/93‑1 (MQ=255)

                                                               cATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCa                        <  2:1916976/93‑1 (MQ=255)

                                                                aTCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAt                       >  2:498622/1‑93 (MQ=255)

                                                                 tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      >  2:4064444/1‑93 (MQ=255)

                                                                 tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      >  1:3438556/1‑93 (MQ=255)

                                                                 acAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      <  1:3045208/92‑1 (MQ=255)

                                                                  cAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCt                     <  1:1094402/93‑1 (MQ=255)

                                                                  cAAAGACCGGCTTTTTTTCCAGGCGACCAGAATAGAGGCCCGGGCGACGCGCGCTCTGGACGTTATTCAACGGCGCATAATCTGCCAACATCt                     <  2:2921760/93‑1 (MQ=255)

                                                                   aaaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTg                    >  2:372133/1‑93 (MQ=255)

                                                                   aaaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTg                    >  2:3800787/1‑93 (MQ=255)

                                                                    aaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGc                   <  1:4696025/93‑1 (MQ=255)

                                                                    aaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGc                   >  2:3047854/1‑93 (MQ=255)

                                                                    aaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGc                   <  2:3205711/93‑1 (MQ=255)

                                                                        ccGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCgg               >  1:2908711/1‑93 (MQ=255)

                                                                         cGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgc              <  2:259310/93‑1 (MQ=255)

                                                                         cGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgc              >  2:1619677/1‑93 (MQ=255)

                                                                         cGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgc              <  1:3829956/93‑1 (MQ=255)

                                                                          ggCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgcg             >  2:3146190/1‑93 (MQ=255)

                                                                           gCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgcgc            <  2:3818902/93‑1 (MQ=255)

                                                                           gCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgcgc            <  2:611066/93‑1 (MQ=255)

                                                                             tGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTg          >  1:2719717/1‑93 (MQ=255)

                                                                                ttttCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTgcgc       >  2:925879/1‑93 (MQ=255)

                                                                                 tttCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCa      <  2:103835/93‑1 (MQ=255)

                                                                                 tttCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCa      >  2:3127481/1‑93 (MQ=255)

                                                                                   tCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCTCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcg    >  2:3259397/1‑93 (MQ=255)

                                                                                   tCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcg    >  2:1709226/1‑93 (MQ=255)

                                                                                    ccAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgc   <  2:1303344/93‑1 (MQ=255)

                                                                                    ccAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgc   <  2:673508/93‑1 (MQ=255)

                                                                                    ccAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgc   <  2:3692082/93‑1 (MQ=255)

                                                                                    ccAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgc   >  2:2298918/1‑93 (MQ=255)

ACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCC  >  NC_000913/4604744‑4604921

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: `‑`