BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,216,301	A→G	S139S (AGT→AGC)	aer ←	Aer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,216,301	0	A	G	99.1%	1103.4 / ‑1.6	350	S139S (AGT→AGC)	aer	Aer
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (185/162); total (188/162)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

Only 100 of 366 total aligned reads displayed.

CATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCA  >  NC_000913/3216211‑3216391

catcaCTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCt                                                                                          >  1:6536004/1‑93 (MQ=255)

 atcaCTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTg                                                                                         >  1:8798725/1‑93 (MQ=255)

  tcaCTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTgg                                                                                        >  1:5205845/1‑93 (MQ=255)

    ccTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTa                                                                                      <  2:3105349/92‑1 (MQ=255)

    ccTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTa                                                                                      <  2:9215993/92‑1 (MQ=255)

    aCTCCACGCGCCCGCTAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTa                                                                                      >  2:1923251/1‑93 (MQ=255)

     cTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTAc                                                                                     >  2:154381/1‑93 (MQ=255)

     cTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTAc                                                                                     >  2:534254/1‑93 (MQ=255)

     cTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTAc                                                                                     >  1:2585485/1‑93 (MQ=255)

      tCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACg                                                                                    >  1:1409843/1‑93 (MQ=255)

        cACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACTAGGCCTTTATGAATACGCTTGCTGGTACGTc                                                                                  <  1:6043301/93‑1 (MQ=255)

        cACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTc                                                                                  <  2:2625189/93‑1 (MQ=255)

        cACGCCCCCCCCAGCGAAGCGGTAATGAAGGCATTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTc                                                                                  <  1:659188/93‑1 (MQ=255)

         ccgcgCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTcc                                                                                 <  1:1604335/92‑1 (MQ=255)

         aCGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTcc                                                                                 >  2:7897952/1‑93 (MQ=255)

          cgcgCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCg                                                                                <  1:5386259/93‑1 (MQ=255)

           gcgcCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCgg                                                                               <  2:4608471/93‑1 (MQ=255)

                cgccAGCGAAGCGGTAATGAAGGCAGTTTACCCATCCAGTCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTc                                                                          >  1:1346586/1‑93 (MQ=255)

                cgccAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTc                                                                          <  2:2938829/93‑1 (MQ=255)

                 gccAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCa                                                                         <  2:3891802/93‑1 (MQ=255)

                  ccAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCaa                                                                        <  1:656837/93‑1 (MQ=255)

                   cAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAAc                                                                       <  1:8436253/93‑1 (MQ=255)

                       gAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCtt                                                                   >  1:3072052/1‑93 (MQ=255)

                       gAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCtt                                                                   >  2:3650834/1‑93 (MQ=255)

                       gAAGCGGCAACGAAAGCAGTTTCCCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGGACGTCCGGCGTTCAACGCtt                                                                   <  1:2059305/93‑1 (MQ=255)

                         aGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTg                                                                 <  2:3401230/93‑1 (MQ=255)

                         aGCAGTAATGAAGGCGGTTTACCCAGACAGCCTTTACGCACCACCAGGCTTTTATGAAGACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTg                                                                 <  1:5708897/93‑1 (MQ=255)

                            gggaaTGAAGTCAGTTGACCCGGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTAc                                                              <  1:4905973/90‑1 (MQ=255)

                            ggTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTAc                                                              <  1:8287617/93‑1 (MQ=255)

                            ggTAATAAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCCTGCTGGTACGTCCGGCGTTTAACGCTTTggaa                                                              >  2:4741006/1‑90 (MQ=255)

                            cgTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTAc                                                              <  1:5231065/92‑1 (MQ=255)

                              tAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAg                                                            <  1:4702606/93‑1 (MQ=255)

                              tAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAg                                                            >  1:4849579/1‑93 (MQ=255)

                              tAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAg                                                            >  1:5086049/1‑93 (MQ=255)

                              tAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAg                                                            >  2:205573/1‑93 (MQ=255)

                               aaTGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGc                                                           >  1:2982723/1‑93 (MQ=255)

                               aaTGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGc                                                           <  2:7617974/93‑1 (MQ=255)

                                 tGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCgg                                                         >  1:9560481/1‑93 (MQ=255)

                                   aaGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCt                                                       >  2:2464945/1‑93 (MQ=255)

                                    aGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTc                                                      >  1:3683745/1‑93 (MQ=255)

                                    aGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTc                                                      <  2:6642510/93‑1 (MQ=255)

                                         gTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGGGGCTCCAccg                                                 >  1:8545524/1‑93 (MQ=255)

                                            tACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCAccgccg                                              >  1:8320228/1‑93 (MQ=255)

                                              cccAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCcgcg                                            >  1:7559142/1‑93 (MQ=255)

                                               ccagccagCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGa                                           >  1:6964764/1‑93 (MQ=255)

                                               ccagccagCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGa                                           <  2:3699420/93‑1 (MQ=255)

                                               ccagccagCCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGa                                           <  2:7494870/93‑1 (MQ=255)

                                                cagccagcCTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGAt                                          >  2:3729650/1‑93 (MQ=255)

                                                 agccagccTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGAtc                                         >  2:448769/1‑93 (MQ=255)

                                                  gccagccTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGAtct                                        >  1:2910120/1‑93 (MQ=255)

                                                  gccagccTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGAtct                                        <  1:8106844/93‑1 (MQ=255)

                                                  gccagccTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGAtct                                        >  1:5508867/1‑93 (MQ=255)

                                                      gccTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTc                                    <  2:1893868/93‑1 (MQ=255)

                                                       ccTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCa                                   <  1:631680/93‑1 (MQ=255)

                                                       ccTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCa                                   >  1:1855062/1‑93 (MQ=255)

                                                        cTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCAt                                  >  1:5918178/1‑93 (MQ=255)

                                                        cTTTACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCAt                                  >  1:3538784/1‑93 (MQ=255)

                                                         tttACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATc                                 <  1:4285514/93‑1 (MQ=255)

                                                         tttACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATc                                 >  2:145491/1‑93 (MQ=255)

                                                         tttACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATc                                 >  1:1451651/1‑93 (MQ=255)

                                                         tttACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATc                                 >  2:5642262/1‑93 (MQ=255)

                                                         tttACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATc                                 >  2:5808444/1‑93 (MQ=255)

                                                          ttACGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATcc                                <  2:5398857/93‑1 (MQ=255)

                                                             cGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTc                             <  2:6143927/93‑1 (MQ=255)

                                                             cGCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTc                             <  1:7052702/93‑1 (MQ=255)

                                                              gCACCACCAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCg                            >  2:716805/1‑93 (MQ=255)

                                                               caccacCAGGCCTTTATGGATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGc                           >  1:2750726/1‑93 (MQ=255)

                                                                accaccAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGcc                          <  2:6075287/93‑1 (MQ=255)

                                                                accaccAGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCAGTCGcc                          >  2:1652016/1‑93 (MQ=255)

                                                                 ccaccaGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGccc                         <  1:9861416/93‑1 (MQ=255)

                                                                  caccaGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCg                        <  2:6294914/93‑1 (MQ=255)

                                                                   cccaGTCGTTCATGAAGACGCTGGCTTGTACGTCCGGCGATCGACGCTTTGTACAGCTGCTCAGCCGCGGCGATCTCCGCTTCCGTCTCCCgg                       <  2:5893482/92‑1 (MQ=255)

                                                                   accaGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCgg                       >  2:9080428/1‑93 (MQ=255)

                                                                     caGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGt                     >  2:5254364/1‑93 (MQ=255)

                                                                      aGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCTGGtc                    >  1:2270231/1‑92 (MQ=255)

                                                                      aGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTa                    >  2:9716922/1‑93 (MQ=255)

                                                                      aGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTa                    <  1:9247920/93‑1 (MQ=255)

                                                                      aGGCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTa                    >  1:7823069/1‑93 (MQ=255)

                                                                        gCCTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACg                  <  1:7249803/93‑1 (MQ=255)

                                                                         ccTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGa                 >  2:1030722/1‑93 (MQ=255)

                                                                         ccTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGa                 <  1:9331012/93‑1 (MQ=255)

                                                                          cTTTATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGaa                >  1:9043860/1‑93 (MQ=255)

                                                                           tttATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAAt               <  1:2518625/93‑1 (MQ=255)

                                                                           tttATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAAt               >  2:8449135/1‑93 (MQ=255)

                                                                           tttATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAAt               <  2:2443425/93‑1 (MQ=255)

                                                                           tttATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAAt               <  1:3984645/93‑1 (MQ=255)

                                                                           tttATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAAt               >  2:6873923/1‑93 (MQ=255)

                                                                            ttATGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATc              <  1:676743/93‑1 (MQ=255)

                                                                             tATGAATACGCTTGCTGGTACGTCCAGCCTTCAACGGTTTGGCCAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCg             <  2:8493382/93‑1 (MQ=255)

                                                                               tGAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGAc           >  2:4984394/1‑93 (MQ=255)

                                                                                gAATACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACa          <  1:3782807/93‑1 (MQ=255)

                                                                                  aTACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACata        <  2:6361103/93‑1 (MQ=255)

                                                                                  aTACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACata        <  2:4868778/93‑1 (MQ=255)

                                                                                  aTACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACata        >  1:1056196/1‑93 (MQ=255)

                                                                                   tACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACatat       >  2:3826900/1‑93 (MQ=255)

                                                                                   tACGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACatat       <  2:9385534/93‑1 (MQ=255)

                                                                                    aCGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACAtata      <  1:5580975/93‑1 (MQ=255)

                                                                                    aCGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACAtata      <  1:5335936/93‑1 (MQ=255)

                                                                                    aCGCTTGCTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACAtata      >  2:4085985/1‑93 (MQ=255)

                                                                                        ttACTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTAGTAATCGACATATAGCCa  >  1:5823561/1‑93 (MQ=255)

CATCACTCCACGCGCCCGCCAGCGAAGCGGTAATGAAGGCAGTTTACCCAGCCAGCCTTTACGCACCACCAGGCCTTTATGAATACGCTTACTGGTACGTCCGGCGTTCAACGCTTTGTACAGCGGCTCCACCGCCGCGATCTCTTCATCCGTCGCCCGGGTACGAATCGACATATAGCCA  >  NC_000913/3216211‑3216391

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