BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	645,034	A→G	V257A (GTG→GCG)	ybdS ←	YbdS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	645,034	0	A	G	99.2%	405.1 / NA	127	V257A (GTG→GCG)	ybdS	YbdS
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (59/67); total (59/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 131 total aligned reads displayed.

TCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATACC  >  NC_000913/644944‑645122

tctctGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGcgcg                                                                                        <  2:4412469/93‑1 (MQ=255)

tctctGCGTGTCAGCGCACCCATGGTTTTAAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCAGGTTTGTACAGcgcg                                                                                        >  2:1115547/1‑93 (MQ=255)

  tctGCGGGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTa                                                                                      <  1:1029432/93‑1 (MQ=255)

   ctGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAg                                                                                     <  1:308692/93‑1 (MQ=255)

    tGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAgg                                                                                    >  2:2949369/1‑93 (MQ=255)

       gtgtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGaaa                                                                                 >  2:960935/1‑93 (MQ=255)

       gtgtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGaaa                                                                                 >  1:4419823/1‑93 (MQ=255)

        tgtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAg                                                                                <  2:2122751/93‑1 (MQ=255)

         gtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGc                                                                               <  1:97451/93‑1 (MQ=255)

           cAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCa                                                                             <  1:2376042/93‑1 (MQ=255)

           cAGCACACCCATGGTTTTCAGTTCATCACCCGCCCAGGATGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCa                                                                             <  2:1614574/93‑1 (MQ=255)

            aGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCAc                                                                            >  2:3566263/1‑93 (MQ=255)

            aGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCAc                                                                            <  1:3980461/93‑1 (MQ=255)

             gcgcACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACg                                                                           <  1:2713563/93‑1 (MQ=255)

             gcgcACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACg                                                                           >  1:3116475/1‑93 (MQ=255)

              cgcACCCATGGTTTTCTGTTCATCCCCCGCCCCGGTTGCCACTTATTCCCTGTGTGTGCTTTCCGGTTAATCCCGCGCGTAGGAAAACCACgg                                                                          <  2:3617193/93‑1 (MQ=255)

              cgcACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACgg                                                                          >  2:3425531/1‑93 (MQ=255)

               gcACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGgc                                                                         >  2:112952/1‑93 (MQ=255)

                   ccATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCaa                                                                     >  1:3141913/1‑93 (MQ=255)

                   ccATGGTTTTCAGTTCATCACCCGCCCAGGGTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGACACGGCGCaa                                                                     >  2:918606/1‑93 (MQ=255)

                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                    >  1:229010/1‑93 (MQ=255)

                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                    >  2:385373/1‑93 (MQ=255)

                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                    >  2:4359651/1‑93 (MQ=255)

                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                    >  1:1686997/1‑93 (MQ=255)

                     aTGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATg                                                                   >  1:4432629/1‑93 (MQ=255)

                      tGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGa                                                                  <  1:1845654/93‑1 (MQ=255)

                       ggTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGAt                                                                 >  1:4117234/1‑93 (MQ=255)

                        gTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATa                                                                <  1:3557188/93‑1 (MQ=255)

                         ttttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATaa                                                               >  2:4520375/1‑93 (MQ=255)

                         ttttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATaa                                                               >  1:903444/1‑93 (MQ=255)

                         ttttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATaa                                                               <  1:2185823/93‑1 (MQ=255)

                          tttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAg                                                              <  1:428711/93‑1 (MQ=255)

                          tttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAg                                                              >  2:4332089/1‑93 (MQ=255)

                            tCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCa                                                            <  2:772688/93‑1 (MQ=255)

                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                         <  2:2002229/93‑1 (MQ=255)

                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                         <  2:3070328/93‑1 (MQ=255)

                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                         <  1:2661472/93‑1 (MQ=255)

                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                         <  2:3534623/93‑1 (MQ=255)

                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                         <  2:3607229/93‑1 (MQ=255)

                                ttCATCACCCGCCCATGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                        >  2:331274/1‑93 (MQ=255)

                                ttCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                        <  2:2717210/93‑1 (MQ=255)

                                ttCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                        <  1:899672/93‑1 (MQ=255)

                                tgcatcaCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                        <  1:899671/91‑1 (MQ=255)

                                 tcatcaCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATa                                                       <  1:2606289/93‑1 (MQ=255)

                                  catcaCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATaa                                                      <  1:414870/93‑1 (MQ=255)

                                  catcaCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATaa                                                      <  2:3781832/93‑1 (MQ=255)

                                     caCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAAccc                                                   <  2:1292348/93‑1 (MQ=255)

                                        ccGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCaa                                                >  1:1782329/1‑93 (MQ=255)

                                         cGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAAc                                               <  2:3470052/93‑1 (MQ=255)

                                         cGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAAc                                               <  1:4063590/93‑1 (MQ=255)

                                         agCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAAc                                               >  2:4370933/2‑93 (MQ=255)

                                          gCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAAcc                                              >  1:1335374/1‑93 (MQ=255)

                                          gCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAAcc                                              >  1:3428761/1‑93 (MQ=255)

                                           cccAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCg                                             <  1:2523122/93‑1 (MQ=255)

                                             cAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGc                                           >  1:1912702/1‑93 (MQ=255)

                                              aGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCa                                          <  1:3016482/93‑1 (MQ=255)

                                               ggTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAg                                         <  1:2879262/93‑1 (MQ=255)

                                               ggTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAg                                         <  2:217711/93‑1 (MQ=255)

                                               ggTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCACAACCGGCAg                                         >  2:4651998/1‑93 (MQ=255)

                                                  tGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGaa                                      <  2:48708/93‑1 (MQ=255)

                                                  tGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGaa                                      <  1:4240351/93‑1 (MQ=255)

                                                  tGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGaa                                      <  1:1548664/93‑1 (MQ=255)

                                                   gCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGaac                                     >  1:2742452/1‑92 (MQ=255)

                                                   gCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAg                                     >  2:1218436/1‑93 (MQ=255)

                                                   gCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAg                                     >  1:1458133/1‑93 (MQ=255)

                                                     cACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCa                                   <  2:3690436/93‑1 (MQ=255)

                                                      aCTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCag                                  >  1:3220319/1‑93 (MQ=255)

                                                       cttcttCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCaga                                 <  2:355847/93‑1 (MQ=255)

                                                       cttcttCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCaga                                 >  2:716267/1‑93 (MQ=255)

                                                        ttcttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCagag                                >  2:1261716/1‑93 (MQ=255)

                                                         tcttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAgg                               <  1:1866645/93‑1 (MQ=255)

                                                         tcttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAgg                               <  1:1176975/93‑1 (MQ=255)

                                                         tcttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAgg                               <  2:1742767/93‑1 (MQ=255)

                                                          cttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGa                              >  1:3915036/1‑93 (MQ=255)

                                                           ttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGaa                             >  2:1723577/1‑93 (MQ=255)

                                                           ttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGaa                             <  1:586163/93‑1 (MQ=255)

                                                           ttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGaa                             >  1:1966619/1‑93 (MQ=255)

                                                              aCTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCa                          >  1:385590/1‑93 (MQ=255)

                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                        >  1:3862454/1‑93 (MQ=255)

                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                        <  1:4691617/93‑1 (MQ=255)

                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                        <  1:3303356/93‑1 (MQ=255)

                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGACAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                        >  1:3862453/1‑93 (MQ=255)

                                                                 gtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTg                       >  2:1391516/1‑93 (MQ=255)

                                                                    tgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCaa                    >  1:3992091/1‑93 (MQ=255)

                                                                     gtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAAc                   <  2:4232831/93‑1 (MQ=255)

                                                                        aTTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAg                <  1:2876767/93‑1 (MQ=255)

                                                                          ttCCGGTTTGTACAGCGGGGAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCt              <  1:530509/93‑1 (MQ=255)

                                                                            ccGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTaa            >  2:696629/1‑93 (MQ=255)

                                                                            ccGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTaa            <  2:3180152/93‑1 (MQ=255)

                                                                             cGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAAt           >  1:3278587/1‑93 (MQ=255)

                                                                              ggTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTaact          >  1:3595577/1‑91 (MQ=255)

                                                                              ggTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATc          >  2:2635765/1‑93 (MQ=255)

                                                                              ggTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATc          <  2:3996094/93‑1 (MQ=255)

                                                                               gTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCt         <  2:557125/93‑1 (MQ=255)

                                                                               gTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCt         >  2:3024665/1‑93 (MQ=255)

                                                                                tttGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTg        <  2:3030736/93‑1 (MQ=255)

                                                                                 ttGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTgg       <  1:1343876/93‑1 (MQ=255)

                                                                                  tGTACAGCGCGTAGGAAAGCCACAGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGa      >  1:2409461/1‑93 (MQ=255)

                                                                                     aCAGCGCGTAGGAAAACCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATAc   >  2:4515028/1‑93 (MQ=255)

                                                                                      cAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATAcc  <  2:919485/93‑1 (MQ=255)

TCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATACC  >  NC_000913/644944‑645122

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