BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	620,296	A→G	V39A (GTA→GCA)	fepG ←	FepG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	620,296	0	A	G	97.7%	859.2 / ‑5.3	263	V39A (GTA→GCA)	fepG	FepG
Reads supporting (aligned to +/- strand): ref base A (0/2); major base G (133/122); minor base T (0/6); total (133/130)
Fisher's exact test for biased strand distribution p-value = 1.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 267 total aligned reads displayed.

ATCAACAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATACCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGTGATG  >  NC_000913/620206‑620385

aTCAACAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGcc                                                                                         >  1:208030/1‑93 (MQ=255)

   aaCAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGc                                                                                      <  2:4877766/93‑1 (MQ=255)

    aCAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCg                                                                                     >  2:5555521/1‑93 (MQ=255)

      aGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAg                                                                                   <  2:3494432/93‑1 (MQ=255)

       gcgcCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAgg                                                                                  >  1:1003948/1‑93 (MQ=255)

       gcgcCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAgg                                                                                  >  2:5245175/1‑93 (MQ=255)

       gcgcCAGCAGCGCGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAgg                                                                                  <  2:2775058/93‑1 (MQ=255)

         gcCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAgggt                                                                                >  1:5028664/1‑91 (MQ=255)

         gcCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGtt                                                                                >  1:3579793/1‑93 (MQ=255)

         gcCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGtt                                                                                >  2:6316316/1‑93 (MQ=255)

         gcCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGtt                                                                                >  1:3775159/1‑93 (MQ=255)

          ccATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGttt                                                                               >  1:4647104/1‑93 (MQ=255)

            aTCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTcc                                                                             <  1:1253730/93‑1 (MQ=255)

            aTCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTcc                                                                             <  2:2543256/93‑1 (MQ=255)

            aTCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTcc                                                                             <  2:6110349/93‑1 (MQ=255)

             tCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCa                                                                            <  2:1438735/93‑1 (MQ=255)

              cAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTTCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAg                                                                           <  1:212112/93‑1 (MQ=255)

                gCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCg                                                                         <  2:4241478/93‑1 (MQ=255)

                 cACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGt                                                                        <  1:3180092/93‑1 (MQ=255)

                    gcgTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGAc                                                                     <  1:2818024/93‑1 (MQ=255)

                    gcgTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGAc                                                                     >  1:5422857/1‑93 (MQ=255)

                      gTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACgg                                                                   >  2:2822578/1‑93 (MQ=255)

                      gTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACgg                                                                   <  1:878017/93‑1 (MQ=255)

                      gTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACgg                                                                   <  1:3221711/93‑1 (MQ=255)

                      gTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACgg                                                                   <  1:1806593/93‑1 (MQ=255)

                       tGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGc                                                                  <  1:3651308/93‑1 (MQ=255)

                       tGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGc                                                                  <  1:887740/93‑1 (MQ=255)

                        ggTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCa                                                                 >  1:1185506/1‑93 (MQ=255)

                           aaaCGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCg                                                              >  1:1476946/1‑93 (MQ=255)

                            aaCGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGc                                                             <  2:3811755/93‑1 (MQ=255)

                            aaCGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGc                                                             <  1:545668/93‑1 (MQ=255)

                            aaCGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGAAGGCACCGc                                                             >  1:5279040/1‑93 (MQ=255)

                              cGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTg                                                           >  1:486292/1‑93 (MQ=255)

                               gCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGc                                                          >  2:3973136/1‑93 (MQ=255)

                               gCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGc                                                          >  1:450640/1‑93 (MQ=255)

                               cccATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATTCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGc                                                          <  1:2146558/92‑1 (MQ=255)

                                ccATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCg                                                         <  2:3768616/93‑1 (MQ=255)

                                ccATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGAAGGCACCGCTGCg                                                         >  2:4833029/1‑93 (MQ=255)

                                  aTTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTa                                                       >  2:2101876/1‑93 (MQ=255)

                                   ttCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTaa                                                      <  2:6253523/93‑1 (MQ=255)

                                   ttCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTaa                                                      >  2:5040183/1‑93 (MQ=255)

                                     cGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATc                                                    >  1:6204586/1‑93 (MQ=255)

                                     cGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATc                                                    >  2:5742572/1‑93 (MQ=255)

                                     cGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATc                                                    <  2:4617480/93‑1 (MQ=255)

                                     cGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATc                                                    <  1:3074577/93‑1 (MQ=255)

                                      ggTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATcc                                                   >  2:3535505/1‑93 (MQ=255)

                                      ggTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATcc                                                   >  2:845934/1‑93 (MQ=255)

                                       gTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATccc                                                  >  1:3920841/1‑93 (MQ=255)

                                          accaccATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAg                                               >  1:6379191/1‑93 (MQ=255)

                                            caccaTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGAt                                             >  2:5387564/1‑93 (MQ=255)

                                              ccaTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATAc                                           >  2:688588/1‑93 (MQ=255)

                                              ccaTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATAc                                           <  1:4454378/93‑1 (MQ=255)

                                               caTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATAcc                                          <  2:5419608/93‑1 (MQ=255)

                                               caTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATAcc                                          >  1:6292523/1‑93 (MQ=255)

                                                aTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCt                                         <  1:2202235/93‑1 (MQ=255)

                                                aTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCt                                         >  2:5932003/1‑93 (MQ=255)

                                                aTCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCt                                         >  2:1956283/1‑93 (MQ=255)

                                                 tCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTg                                        <  1:5645345/93‑1 (MQ=255)

                                                 tCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTg                                        <  2:3018886/93‑1 (MQ=255)

                                                   gTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCa                                      >  2:4756634/1‑93 (MQ=255)

                                                   gTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCa                                      <  1:4618219/93‑1 (MQ=255)

                                                    tCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCaa                                     >  2:4095716/1‑93 (MQ=255)

                                                    tCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCaa                                     >  1:1902861/1‑93 (MQ=255)

                                                     cATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAAc                                    <  1:5165776/93‑1 (MQ=255)

                                                      aTACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAAcc                                   >  2:431787/1‑93 (MQ=255)

                                                      aTACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAAcc                                   <  1:2748228/93‑1 (MQ=255)

                                                      aTACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAAcc                                   <  2:5634418/93‑1 (MQ=255)

                                                       tACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACca                                  >  1:1713679/1‑93 (MQ=255)

                                                        aCTGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACcac                                 <  2:1389865/93‑1 (MQ=255)

                                                          tGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACcacac                               >  1:2543022/1‑93 (MQ=255)

                                                          tGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACcacac                               >  2:3168214/1‑93 (MQ=255)

                                                          tGCGCGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACcacac                               >  2:6422041/1‑93 (MQ=255)

                                                           gcgcGGCGCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACcacaca                              >  1:5598243/1‑93 (MQ=255)

                                                           gcgcGGCGCATCGCCCATCAGCGCGGCGAATACCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACcacaca                              <  2:2883808/93‑1 (MQ=255)

                                                             gcggcgCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGGTGCGTAATCCCCAGATACCTGCAACCACACAgg                            <  2:1807992/93‑1 (MQ=255)

                                                              cggcgCATCGCCCATCAGCGCGGCGAATTCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGc                           <  2:542812/93‑1 (MQ=255)

                                                              cggcgCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGc                           <  1:152735/93‑1 (MQ=255)

                                                              cggcgCATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGc                           >  1:3041894/1‑93 (MQ=255)

                                                                gcgcATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCgg                         <  1:6083935/93‑1 (MQ=255)

                                                                  gcATCGCACATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATCCCTGCAACCCCACAGGCGGaa                       <  1:5828395/93‑1 (MQ=255)

                                                                   cATCGCCCATCAGCGCGGCGAATTCCTGCTAGGTTTCCAGCGTGGCGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGaaa                      <  1:3329912/93‑1 (MQ=255)

                                                                   cATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCCGATACCTGCAACCACACAGGCGGaaa                      >  2:4666871/1‑93 (MQ=255)

                                                                   cATCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGaaa                      <  2:2966827/93‑1 (MQ=255)

                                                                    aTCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAAc                     <  1:5345392/93‑1 (MQ=255)

                                                                     tCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAAcc                    >  1:3430957/1‑93 (MQ=255)

                                                                     tCGCCCATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAAcc                    <  1:2922663/93‑1 (MQ=255)

                                                                       ccccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcag                  <  1:6589300/92‑1 (MQ=255)

                                                                        cccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagc                 >  2:4398431/1‑93 (MQ=255)

                                                                         ccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagca                >  2:3045727/1‑93 (MQ=255)

                                                                         ccATCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACcagca                >  1:2671083/1‑93 (MQ=255)

                                                                           aTCAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCagcagc              >  1:4181750/1‑93 (MQ=255)

                                                                             cAGCGCGGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCaa            <  1:2367355/93‑1 (MQ=255)

                                                                                cgcgGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACa         >  1:4251001/1‑93 (MQ=255)

                                                                                cgcgGCGAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACa         >  1:6509623/1‑93 (MQ=255)

                                                                                 gcggcgAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAg        >  2:3739302/1‑93 (MQ=255)

                                                                                 gcggcgAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAg        <  2:1055294/93‑1 (MQ=255)

                                                                                  cggcgAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCCGCAACCACCCAGGCGGAAACCAGCAGCAAACAgg       >  2:219134/1‑93 (MQ=255)

                                                                                      gAATGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGtgat   >  2:2677081/1‑93 (MQ=255)

                                                                                       aaTGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGtgatg  <  1:6653725/93‑1 (MQ=255)

                                                                                       aaTGCCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGGAATCCCCAGATACCTGCAACCACACAGGCGGAAAACAGCAGCAAACAGGtgatg  >  2:5056073/1‑93 (MQ=255)

ATCAACAGCGCCATCAGCACGCGTGGTAAACGCCATTCGGTGACCACCATCGTCATACTGCGCGGCGCATCGCCCATCAGCGCGGCGAATACCTGCGAGGTTTCCAGCGTGACGGCACCGCTGCGTAATCCCCAGATACCTGCAACCACACAGGCGGAAACCAGCAGCAAACAGGTGATG  >  NC_000913/620206‑620385

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