BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	568,823	A→G	G233G (GGA→GGG)	ybcK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	568,823	0	A	G	99.5%	584.7 / ‑3.8	183	coding (699/1527 nt)	ybcK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (92/90); minor base . (1/0); total (93/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 185 total aligned reads displayed.

CAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGA  >  NC_000913/568733‑568912

cAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGcc                                                                                         <  2:1971558/93‑1 (MQ=255)

 aGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCt                                                                                        >  1:1047672/1‑93 (MQ=255)

 aGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCt                                                                                        <  1:243347/93‑1 (MQ=255)

  ggATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCtt                                                                                       >  2:2102029/1‑93 (MQ=255)

   gATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTc                                                                                      <  1:1107642/93‑1 (MQ=255)

     tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTg                                                                                    >  1:1629600/1‑93 (MQ=255)

      ggAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGt                                                                                   >  1:1698597/1‑93 (MQ=255)

      ggAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGt                                                                                   <  2:4551311/93‑1 (MQ=255)

      ggAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGt                                                                                   <  1:4436913/93‑1 (MQ=255)

       gAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTa                                                                                  <  1:1471098/93‑1 (MQ=255)

       gAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTa                                                                                  >  2:953021/1‑93 (MQ=255)

       gAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTa                                                                                  <  1:1519117/93‑1 (MQ=255)

        aaaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTaa                                                                                 <  1:279097/93‑1 (MQ=255)

        aaaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTaa                                                                                 <  2:3879561/93‑1 (MQ=255)

         aaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAt                                                                                >  2:3652613/1‑93 (MQ=255)

          aGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAtt                                                                               <  1:3060300/93‑1 (MQ=255)

          aGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAtt                                                                               <  1:3471334/93‑1 (MQ=255)

            gcgcTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGa                                                                             >  1:26047/1‑93 (MQ=255)

              gcTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaa                                                                           >  1:4258190/1‑93 (MQ=255)

               cTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaa                                                                          >  2:1749386/1‑93 (MQ=255)

               cTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaa                                                                          <  1:3762996/93‑1 (MQ=255)

               cTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaa                                                                          <  1:1154502/93‑1 (MQ=255)

                tCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaaa                                                                         >  1:899933/1‑93 (MQ=255)

                tCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaaa                                                                         >  2:3178460/1‑93 (MQ=255)

                 cATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAat                                                                        <  1:3438631/93‑1 (MQ=255)

                  aTTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatt                                                                       >  2:4489622/1‑93 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      <  2:3264278/93‑1 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      <  1:4616022/93‑1 (MQ=255)

                    tGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAattat                                                                     <  1:3997298/93‑1 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   <  2:1166105/93‑1 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   <  1:347719/93‑1 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   <  1:2001681/93‑1 (MQ=255)

                       aTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAg                                                                  <  1:2833595/93‑1 (MQ=255)

                       aTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAg                                                                  >  1:2135044/1‑93 (MQ=255)

                         gCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCg                                                                >  1:2553224/1‑93 (MQ=255)

                         gCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCg                                                                >  2:110984/1‑93 (MQ=255)

                         gCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCg                                                                >  2:1045946/1‑93 (MQ=255)

                          cAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGa                                                               >  1:1884110/1‑93 (MQ=255)

                            aTAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAAt                                                             <  2:2567890/93‑1 (MQ=255)

                            aTAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAAt                                                             >  1:2116113/1‑93 (MQ=255)

                               gCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATaaa                                                          >  1:3561304/1‑93 (MQ=255)

                               gCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATaaa                                                          <  2:334666/93‑1 (MQ=255)

                               gCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATaaa                                                          >  1:3995475/1‑93 (MQ=255)

                                ctaaGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGG‑CCTTCTGTAATTGAAAAATTATTAGCGAATAAAGc                                                        >  1:1485187/3‑93 (MQ=255)

                                 cAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGc                                                        >  2:873193/1‑93 (MQ=255)

                                  aaGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGct                                                       >  2:536071/1‑93 (MQ=255)

                                   aGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctc                                                      <  2:3976462/93‑1 (MQ=255)

                                   aGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctc                                                      >  1:1054572/1‑93 (MQ=255)

                                   aGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctc                                                      <  1:2213322/93‑1 (MQ=255)

                                     tATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTg                                                    <  2:1278927/93‑1 (MQ=255)

                                      ttttAAATCATCATGCTTTACAGAATTTCTCAGGAACACAAAGTGCATGGGGGCCTTCCGTAATTGAAAAATTATTAGCGAATAAAGCTCCGa                                                   <  2:3130014/92‑1 (MQ=255)

                                      aTTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGa                                                   <  1:3352595/93‑1 (MQ=255)

                                      aTTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGa                                                   >  2:4466790/1‑93 (MQ=255)

                                       tttAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGAt                                                  >  1:2293342/1‑93 (MQ=255)

                                        ttAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATa                                                 <  1:1731021/93‑1 (MQ=255)

                                         tAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAg                                                <  1:1351253/93‑1 (MQ=255)

                                         tAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAg                                                >  1:885834/1‑93 (MQ=255)

                                          aaaTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAgg                                               <  1:3014444/93‑1 (MQ=255)

                                            aTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTa                                             <  1:567893/93‑1 (MQ=255)

                                            aTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTa                                             >  1:3605757/1‑93 (MQ=255)

                                            aTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTa                                             <  1:4206810/93‑1 (MQ=255)

                                             tGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAt                                            >  1:4542904/1‑93 (MQ=255)

                                              gATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAtt                                           <  1:3198719/93‑1 (MQ=255)

                                               aTCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAttt                                          <  1:2928260/93‑1 (MQ=255)

                                                tCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTg                                         >  2:2799155/1‑93 (MQ=255)

                                                  aTGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCg                                       >  2:1922948/1‑93 (MQ=255)

                                                  aTGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCg                                       <  1:1194889/93‑1 (MQ=255)

                                                   tGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGt                                      >  1:4437004/1‑93 (MQ=255)

                                                     cTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTAc                                    >  1:141663/1‑93 (MQ=255)

                                                       gTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCt                                  <  1:75048/93‑1 (MQ=255)

                                                       gTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCt                                  >  2:3453698/1‑93 (MQ=255)

                                                        tAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCtt                                 <  1:4161603/93‑1 (MQ=255)

                                                         aaaGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTc                                >  1:2309115/1‑93 (MQ=255)

                                                         aaaGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTc                                <  2:1794112/93‑1 (MQ=255)

                                                          aaGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCa                               <  1:4666654/93‑1 (MQ=255)

                                                           aGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCat                              <  1:3722514/93‑1 (MQ=255)

                                                             aaTTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCatat                            <  1:1326727/93‑1 (MQ=255)

                                                               tttCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGGATTTGCGTACCTTCATATCg                          >  2:3650028/1‑93 (MQ=255)

                                                                ttCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGt                         <  2:656015/93‑1 (MQ=255)

                                                                ttCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGt                         <  1:3835660/93‑1 (MQ=255)

                                                                  ctcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGc                       <  2:3566114/93‑1 (MQ=255)

                                                                   tcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCa                      >  2:145470/1‑93 (MQ=255)

                                                                     aGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAag                    >  1:1855099/1‑93 (MQ=255)

                                                                      ggAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAaga                   <  1:1689436/93‑1 (MQ=255)

                                                                      ggAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAaga                   >  2:2789065/1‑93 (MQ=255)

                                                                        aaaaGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAgg                 >  2:2471421/1‑93 (MQ=255)

                                                                          aagaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGa               >  2:1377654/1‑93 (MQ=255)

                                                                          aagaaagTGCATGGGGGCCATCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGa               >  2:4578366/1‑93 (MQ=255)

                                                                           agaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGaa              <  1:2697425/93‑1 (MQ=255)

                                                                               agTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAggg          >  2:772859/1‑93 (MQ=255)

                                                                                gTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGa         <  2:2847514/93‑1 (MQ=255)

                                                                                 tGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGAt        >  1:4358924/1‑93 (MQ=255)

                                                                                 tGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCGTCATATCGTGCAAGAGGGAAAGGGAt        >  2:3795685/1‑93 (MQ=255)

                                                                                   cATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATaa      <  1:4070089/93‑1 (MQ=255)

                                                                                    aTGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAg     >  2:1391721/1‑93 (MQ=255)

                                                                                    aTGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAg     <  2:1954486/93‑1 (MQ=255)

                                                                                     tGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGt    >  1:3968555/1‑93 (MQ=255)

                                                                                      ggggGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTg   <  1:2143884/93‑1 (MQ=255)

                                                                                      ggggGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTg   <  1:126824/93‑1 (MQ=255)

                                                                                       ggggCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGa  <  1:2499659/93‑1 (MQ=255)

CAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGA  >  NC_000913/568733‑568912

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